Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,378,108 +C coding (47/471 nt) argR → l‑arginine‑responsive arginine metabolism regulon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,378,1071.C96.4% 173.1 / 0.7 56coding (46/471 nt)argRl‑arginine‑responsive arginine metabolism regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (26/27);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

CTCTTTCATAACATTATTTCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑G‑C‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCTGTACGTACACGCAATGCCAAAAT  >  CP009273/3377959‑3378249
                                                                                                                                                      |                                                                                                                                                  
ctctTTCATAACAGTATTTCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGttaaagcatttaaa‑g                                                                                                                                                     <  1:12756/149‑1 (MQ=255)
    ttCATAACATTGTTGCAGCCGTCTGCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCCACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGttaaagcatttaaa‑g‑ccat                                                                                                                                                <  2:507341/149‑4 (MQ=255)
     tCATAACATTATTTCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTCGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGttaaagcatttaca‑g‑ccatt                                                                                                                                               >  1:1111174/1‑145 (MQ=255)
          ccATTCTTTCAGCCTTCTTCAGGGCTGATTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTa                                                                                                                                          <  2:345536/148‑1 (MQ=255)
          aCATTATTTCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTCAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTa                                                                                                                                          >  1:397113/1‑149 (MQ=255)
             ttattTCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑aag                                                                                                                                    >  2:1100389/1‑149 (MQ=255)
             ttattTCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGACTAAGCAAGAAGAACTAGTTAAAGCATTTGAA‑GCC‑ATTACGTA‑‑‑aag                                                                                                                                    >  1:214321/1‑149 (MQ=255)
             ttattTCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGAACAATAATGTTGTATCAACCACCATATCGGGAGACTTATGCGAAGCGCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑GTTACTTA‑‑‑aag                                                                                                                                    >  1:214443/1‑149 (MQ=255)
                 ttCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑Aagaaga                                                                                                                                <  1:575489/149‑1 (MQ=255)
                 ttCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑Aagaaga                                                                                                                                <  2:504192/149‑1 (MQ=255)
                 ttCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑G‑CGATTACTGA‑‑‑Aagaaga                                                                                                                                >  2:308551/1‑149 (MQ=255)
                  tCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAagag                                                                                                                               >  1:356308/1‑149 (MQ=255)
                       cttcttCCGGGCTGTCTGTTTGCCTAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAAtt                                                                                                                          <  1:776387/149‑1 (MQ=255)
                       cttcttCAGGGCGGACGGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAAtt                                                                                                                          <  2:943092/149‑1 (MQ=255)
                       cttcttCAGGGATGAGTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAATCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAAtt                                                                                                                          <  1:911039/149‑1 (MQ=255)
                        ttcttcAGGGCTGGCTGTTTGCATAAAAATTCATCTGTATGCACCATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAGGCTCGGCTAAGCAAGAAGCACTAGTGAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAACAGAGAAActt                                                                                                                         >  1:615681/1‑149 (MQ=255)
                         tcttcAGGGCTGCCTTTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCCCCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTa                                                                                                                        <  2:125511/149‑1 (MQ=255)
                          cttcAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATCCGAAGCTCGGCTACGCCAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATCACTTA‑‑‑AAGAAGAGAAATTTAg                                                                                                                       >  2:813261/1‑149 (MQ=255)
                           ttcAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGc                                                                                                                      >  2:1012567/1‑149 (MQ=255)
                                ggCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCa                                                                                                                 <  1:502502/149‑1 (MQ=255)
                                   tGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAggg                                                                                                              >  2:804392/1‑149 (MQ=255)
                                   tGACTGTTTGCATAAAAATTCATCTGTATGCACAAGAATGTTGTATCAACCACCATATCGGGAGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATAACTTA‑‑‑AAGAAGAGAAATATAGCCCCGAggg                                                                                                              >  1:689475/1‑149 (MQ=255)
                                       tGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGaa                                                                                                          >  2:966857/1‑149 (MQ=255)
                                             cATAAAAATTCATCTTTATGCACAATAATGTTGTATCAACCACGATATCGGGTGAATTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTACA‑GAC‑ATTACCTA‑‑‑AAGAAGAGAACTTTAGCTCCCAGGGCGAAAtcgtc                                                                                                    <  1:562315/149‑1 (MQ=255)
                                               tAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTcgc                                                                                                  <  2:682879/149‑1 (MQ=255)
                                               tAAAAATTCATCCGTACGCACAATAATGCTGTATCAACCACCATAACGGGTGACATATGCGAAGCTCGACTAAGCAAGAAGAACTAGTTAAAGCAATAACA‑GCC‑ACAACACA‑‑‑AAGAAGAGAAATTCAGCTCACAGAGCGAAATCGgcgc                                                                                                  >  1:630125/1‑149 (MQ=255)
                                                          cTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTGAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGCAATTTAGCTCCCAGGGCGAGATCGGCGCCGCGTTGCAga                                                                                       >  1:523899/1‑148 (MQ=255)
                                                            gTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAATAGAAATTTAGATCCCAGGGCGAAATCATCGCCGCGTTGgaggag                                                                                     >  1:1093643/1‑149 (MQ=255)
                                                             tATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGGGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGc                                                                                    >  2:795449/1‑149 (MQ=255)
                                                                gCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAg                                                                                 <  1:62513/149‑1 (MQ=255)
                                                                gCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTAGGCTAAGCAAGACGAACTAGTTAAAGCCTTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTGCCCGCGGTGCCGGAGCAAg                                                                                 >  1:583509/1‑149 (MQ=255)
                                                                  acaATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGCAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGAAGGAGCAAGGc                                                                               >  1:35350/1‑149 (MQ=255)
                                                                   caATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTAAAGAAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAg                                                                                 <  1:153624/149‑1 (MQ=255)
                                                                    aataatGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCtt                                                                             <  2:356308/149‑1 (MQ=255)
                                                                     ataatGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCttt                                                                            <  1:1008767/149‑1 (MQ=255)
                                                                     ataatGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCttt                                                                            <  1:606627/149‑1 (MQ=255)
                                                                           ttgtATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCCGGAGCAAGGCTTTGACact                                                                      >  1:545958/1‑147 (MQ=255)
                                                                            tgtATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAata                                                                     <  2:254927/149‑1 (MQ=255)
                                                                             gtATCACCCACCATATCGGGTGACTTATGCTAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAatat                                                                    <  2:462515/149‑1 (MQ=255)
                                                                               aTCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTa                                                                  <  1:836717/149‑1 (MQ=255)
                                                                                 cAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAAt                                                                <  1:1079646/149‑1 (MQ=255)
                                                                                     caccaTATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGt                                                            <  1:1100389/149‑1 (MQ=255)
                                                                                        caTATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAGGGCTTGGACAATATTAATCAGTCTa                                                         >  1:201608/1‑149 (MQ=255)
                                                                                          tttCGGGTGACTTATGCGAAGCGCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTaaa                                                       <  1:189516/147‑1 (MQ=255)
                                                                                                gTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTcg                                                 >  1:1111720/1‑149 (MQ=255)
                                                                                                  gACTTAGGCGAAGCTCCGCTAAGCAAGAAGAACGATTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGACGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTcgcg                                               <  1:999197/149‑1 (MQ=255)
                                                                                                          cGAAGCGCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTGCTTA‑‑‑TATAATATAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGa                                       <  2:1069902/149‑1 (MQ=255)
                                                                                                           gAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGAc                                      >  2:472864/1‑149 (MQ=255)
                                                                                                                  ggCTAAGCAAGAACAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGCAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGCCTAAAGTCTCGCGGATGTTGACCAAGtgt                               >  1:797479/1‑147 (MQ=255)
                                                                                                                   gCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTg                              >  1:338842/1‑149 (MQ=255)
                                                                                                                   gCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTg                              <  2:233872/149‑1 (MQ=255)
                                                                                                                   gCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAACTCTCGCGGAGGTTGACCAAGTTTg                              >  2:629436/1‑149 (MQ=255)
                                                                                                                    cTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTgg                             <  1:370545/149‑1 (MQ=255)
                                                                                                                        gccagaagaACTAGTTAAAGCATTTAAATC‑C‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGATGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGTATGTTGACCAAGTTTGGTGCt                         >  1:647629/4‑149 (MQ=255)
                                                                                                                        gCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGACATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCt                         <  2:397113/149‑1 (MQ=255)
                                                                                                                        gCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCt                         <  2:31780/149‑1 (MQ=255)
                                                                                                                        gCAAGAAGAACTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCt                         >  1:508614/1‑149 (MQ=255)
                                                                                                                             aagaaCTAGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCTgtacg                    <  1:962232/149‑1 (MQ=255)
                                                                                                                                    aGTTAAAGCATTTAAA‑GCC‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAACCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCTGTACGTACACGc             >  2:322114/1‑149 (MQ=255)
                                                                                                                                       taaagcatttaaa‑gcc‑atTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCTGTACGTACACGCAAt          <  2:496794/149‑1 (MQ=255)
                                                                                                                                            catttaaa‑gcc‑atTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCTGTACGTACACGCAATGCCaa     >  1:326103/1‑149 (MQ=255)
                                                                                                                                               ttaaa‑gcc‑atTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGTGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCTGTACGTACACGCAATGCCAaaat  <  2:214321/149‑1 (MQ=255)
                                                                                                                                               ttaaa‑gcc‑atGACCTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTATAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCTGTACGTACACGCAATGCCAaaat  <  2:214443/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CTCTTTCATAACATTATTTCAGCCTTCTTCAGGGCTGACTGTTTGCATAAAAATTCATCTGTATGCACAATAATGTTGTATCAACCACCATATCGGGTGACTTATGCGAAGCTCGGCTAAGCAAGAAGAACTAGTTAAAGCATTTAAA‑G‑C‑ATTACTTA‑‑‑AAGAAGAGAAATTTAGCTCCCAGGGCGAAATCGTCGCCGCGTTGCAGGAGCAAGGCTTTGACAATATTAATCAGTCTAAAGTCTCGCGGATGTTGACCAAGTTTGGTGCTGTACGTACACGCAATGCCAAAAT  >  CP009273/3377959‑3378249

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑