Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,388,379 C→G 9.2% intergenic (+73/+71) pyrD → / ← rmf dihydroorotate dehydrogenase/ribosome modulation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,388,3790CG 9.2% 223.2 / 9.9 75intergenic (+73/+71)pyrD/rmfdihydroorotate dehydrogenase/ribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base G (5/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAC  >  NC_002947/2388232‑2388526
                                                                                                                                                   |                                                                                                                                                   
taCTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCaa                                                                                                                                                   <  2:456205/150‑1 (MQ=255)
  cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAgg                                                                                                                                                 <  1:556408/150‑1 (MQ=255)
   tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                                                                                                                >  1:1265065/1‑150 (MQ=255)
      ggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGcc                                                                                                                                             >  1:1508599/1‑150 (MQ=255)
        cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccct                                                                                                                                           <  2:619099/150‑2 (MQ=255)
         ttCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCAGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCCTAAAAAAGGGCCCATTGAAGGGGCCCCt                                                                                                                                          >  1:1744393/1‑150 (MQ=255)
           cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  1:1084020/150‑1 (MQ=255)
           cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  2:598663/150‑1 (MQ=255)
            aTCTACAAAGGGACGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTggg                                                                                                                                       >  2:1553282/1‑150 (MQ=255)
                 cAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTc                                                                                                                                  <  2:507419/150‑1 (MQ=255)
                  aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcg                                                                                                                                 <  1:22715/150‑1 (MQ=255)
                    aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgcc                                                                                                                               <  1:531907/150‑1 (MQ=255)
                    aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGAAAGAGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgcc                                                                                                                               <  2:1808580/150‑1 (MQ=255)
                      ggCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAg                                                                                                                                                  <  1:806358/129‑1 (MQ=255)
                      ggCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAg                                                                                                                                                  >  2:806358/1‑129 (MQ=255)
                         cGGCGCTGATTCGCGAAGCGGTGGATGCTATGGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCACGGGGCCCCTGGGCCGTCGCccgccg                                                                                                                          >  1:166665/1‑150 (MQ=255)
                           gcgcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTACAGCTGATGAGGTTTGCCAAGTCATAAAAAAGGGCCCCTTCAAGGGGCCCCTTGGCCTTCGCCCgccgcc                                                                                                                        >  2:214440/1‑150 (MQ=255)
                               tGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGa                                                                                                                    >  1:270299/1‑150 (MQ=255)
                                gATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGAt                                                                                                                   >  1:1675187/1‑150 (MQ=255)
                                 aTTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGGGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATg                                                                                                                  <  1:870412/150‑1 (MQ=255)
                                 aTTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATg                                                                                                                  <  1:1293749/150‑1 (MQ=255)
                                  ttCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATgg                                                                                                                 <  2:1675187/150‑1 (MQ=255)
                                       gAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCtt                                                                                                                                   >  2:692087/1‑127 (MQ=255)
                                       gAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCtt                                                                                                                                   <  1:692087/127‑1 (MQ=255)
                                                      accgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                             >  2:985969/3‑150 (MQ=255)
                                                      aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa                                                                                             <  2:507352/150‑1 (MQ=255)
                                                       ccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        <  1:236152/105‑1 (MQ=255)
                                                       ccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                                                                                                        >  2:236152/2‑106 (MQ=255)
                                                          cgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa                                                                                         >  1:710710/1‑150 (MQ=255)
                                                           ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAg                                                                                        <  1:1355595/150‑1 (MQ=255)
                                                           ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGATAACCCAAAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAg                                                                                        <  1:1553282/150‑1 (MQ=255)
                                                            gCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccctttt                                                                                                                                        >  1:587466/1‑97 (MQ=255)
                                                            gCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccctttt                                                                                                                                        <  2:587466/101‑5 (MQ=255)
                                                            gCAATGCCGCGGTGACCAGCCCGCCCCTGTCGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGTCCCTGGGCCTTCGCCCGCCTCCCGGATGGGGCGTGCATCGCAATTGGTATTCAgt                                                                                       <  2:1019615/150‑1 (MQ=255)
                                                             cAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa                                                                                                  <  1:401449/138‑1 (MQ=255)
                                                             cAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa                                                                                                  >  2:401449/1‑138 (MQ=255)
                                                                  ccGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTc                                                                                 <  1:985969/150‑1 (MQ=255)
                                                                   cgcgGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGc                                                                                                          <  2:942249/124‑1 (MQ=255)
                                                                   cgcgGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGc                                                                                                          >  1:942249/1‑124 (MQ=255)
                                                                                              gggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAATAAAGGGCGCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGa                                                     <  1:1451342/150‑1 (MQ=255)
                                                                                              gggCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGa                                                     <  1:1895254/150‑1 (MQ=255)
                                                                                                    gtgggtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTAt                                               <  2:270299/150‑1 (MQ=255)
                                                                                                      gggtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATgg                                             <  2:1739318/150‑1 (MQ=255)
                                                                                                       ggtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCCGCCAACGGCGTGATTTTCGTTGAGTCTATGGa                                            >  1:852504/1‑150 (MQ=255)
                                                                                                        gtgCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGAt                                           <  2:1265065/150‑1 (MQ=255)
                                                                                                             gCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATgg                                                                                                                 >  1:1461199/1‑75 (MQ=255)
                                                                                                             gCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGACGCCCGGATgg                                                                                                                 <  2:1461199/75‑1 (MQ=255)
                                                                                                                     ggTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCa                              <  2:1815804/150‑1 (MQ=255)
                                                                                                                     ggTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGACTACGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCa                              <  1:1550920/150‑1 (MQ=255)
                                                                                                                      gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGtt                                                       <  1:1627517/124‑1 (MQ=255)
                                                                                                                      gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGtt                                                       >  2:1627517/1‑124 (MQ=255)
                                                                                                                      gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCa                              >  1:1097958/1‑149 (MQ=255)
                                                                                                                      gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCa                              <  2:1097958/149‑1 (MQ=255)
                                                                                                                      gTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAg                             >  1:5638/1‑150 (MQ=255)
                                                                                                                           ccAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCAtt                        >  1:771011/1‑150 (MQ=255)
                                                                                                                            cAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTc                       <  1:1675494/150‑1 (MQ=255)
                                                                                                                            cAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTAGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTc                       >  1:539751/1‑150 (MQ=255)
                                                                                                                              aGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCg                     >  2:1461942/1‑150 (MQ=255)
                                                                                                                               ggCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGt                    <  1:368207/150‑1 (MQ=255)
                                                                                                                                gCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGCGCCAGTCATTCCGTc                   <  2:1419817/150‑1 (MQ=255)
                                                                                                                                  aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTccc                 >  2:1860931/1‑150 (MQ=255)
                                                                                                                                  aTAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTccc                 >  1:892740/1‑150 (MQ=255)
                                                                                                                                  aTAAAAAAGGGACCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTccc                 >  2:1860664/1‑150 (MQ=255)
                                                                                                                                      ccaggggCCCCTTGAAGGGGCCCTTGGGCATTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTGCAGTGTCTGGGTGCCAGCGGTGCCAGTCATTCCGTCGCAGtt             >  2:722142/5‑150 (MQ=255)
                                                                                                                                          gggCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTTATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCa         >  2:611531/1‑150 (MQ=255)
                                                                                                                                            gCCCCTTGAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCAcc       >  1:884006/1‑150 (MQ=255)
                                                                                                                                            gCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCAcc       >  1:923929/1‑150 (MQ=255)
                                                                                                                                             ccccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcg      <  2:710710/150‑1 (MQ=255)
                                                                                                                                             ccccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcg      >  1:650597/1‑150 (MQ=255)
                                                                                                                                              cccTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACcgc     >  2:191829/1‑150 (MQ=255)
                                                                                                                                                cTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCACACCCGTGATTTTCGTTGAGTCTATTGAGGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGa   >  1:697659/1‑150 (MQ=255)
                                                                                                                                                cTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGa   >  1:1920646/1‑150 (MQ=255)
                                                                                                                                                 ttCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAc  >  2:1602148/1‑150 (MQ=255)
                                                                                                                                                 ttCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCATTCAGCCAACGGCGTGATTTTCGTTGAGt                                                   >  2:1863256/1‑101 (MQ=255)
                                                                                                                                                 ttCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCATTCAGCCAACGGCGTGATTTTCGTTGAGt                                                   <  1:1863256/101‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAGTGTCCTCGTTCAGCCAACGGCGTGATTTTCGTTGAGTCTATGGATGCCAGCGGTGCCAGTCATTCCGTCCCAGTTATCACCGCGAC  >  NC_002947/2388232‑2388526

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