Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,388,295 A→T 7.2% M338L (ATG→TTG)  pyrD → dihydroorotate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,388,2950AT 7.2% 345.7 / 14.6 111M338L (ATG→TTG) pyrDdihydroorotate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (47/55);  new base T (4/4);  total (52/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

Only 100 of 113 total aligned reads displayed.

TTGCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAG  >  NC_002947/2388148‑2388440
                                                                                                                                                   |                                                                                                                                                 
ttGCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATg                                                                                                                                                 >  1:2145208/1‑150 (MQ=255)
     gATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCcgcg                                                                                                                                            >  2:1680504/1‑150 (MQ=255)
     gATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTACCGAGAAGATCGCGTCCGGGGCGAGCTTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGATGGATGCTATCACGTCACTGCcgcg                                                                                                                                            >  1:2299154/1‑150 (MQ=255)
       tCATTGCCGCCGGGGGCATTAGCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGt                                                                                                                                          <  2:1239825/150‑1 (MQ=255)
       tCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGt                                                                                                                                          <  2:729029/150‑1 (MQ=255)
          ttGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTACCGCGGCATTGCcgcgatagc                                                                                                                                       <  2:395290/150‑6 (MQ=255)
                ccgGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAgccc                                                                                                                                 <  2:1431420/150‑1 (MQ=255)
                    gggCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAgcccgccc                                                                                                                             <  1:240441/150‑1 (MQ=255)
                    gggCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAgcccgccc                                                                                                                             >  1:2104582/1‑150 (MQ=255)
                        aTTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCtgt                                                                                                                         <  1:1380273/150‑1 (MQ=255)
                          tACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTgg                                                                                                                       >  2:1874444/1‑150 (MQ=255)
                          tACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTgg                                                                                                                       <  1:1956865/150‑1 (MQ=255)
                           aCCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggg                                                                                                                      <  2:370450/150‑1 (MQ=255)
                            ccGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGGAGATATACTAGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCTGTGACCAGCCAGCCCCCTTgggg                                                                                                                     >  2:506855/1‑150 (MQ=255)
                            ccGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGACCGCCCCTGTgggg                                                                                                                     <  2:2299154/150‑1 (MQ=255)
                             cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggggg                                                                                                                    >  2:692194/1‑150 (MQ=255)
                             cGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCACACCTGTgtggg                                                                                                                    >  2:172997/1‑150 (MQ=255)
                                    cgccATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggt                                                                                                             <  1:1077701/150‑1 (MQ=255)
                                     gccATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtg                                                                                                            >  1:402449/1‑150 (MQ=255)
                                      ccATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtgg                                                                                                           >  2:769568/1‑150 (MQ=255)
                                       cATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtggg                                                                                                          >  2:1608103/1‑150 (MQ=255)
                                        aTGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtgggt                                                                                                         <  2:283825/150‑1 (MQ=255)
                                        aTGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCgggtgggt                                                                                                         >  1:283825/1‑150 (MQ=255)
                                               cGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGAccg                                                                                                                                     >  1:1631390/1‑114 (MQ=255)
                                               cGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGAccg                                                                                                                                     <  2:1631390/115‑2 (MQ=255)
                                                gagaAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtga                                                                                                 >  1:2226962/1‑150 (MQ=255)
                                                  gaagaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCAGCCCCGGCGGGGGGACGGGGGGGGGCGGCtgatg                                                                                               >  2:2323285/1‑150 (MQ=255)
                                                   aagaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtgatga                                                                                              <  1:769568/150‑1 (MQ=255)
                                                   aagaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtgatga                                                                                              <  2:402449/150‑1 (MQ=255)
                                                    agaTCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAg                                                                                             <  2:1660336/150‑1 (MQ=255)
                                                       tCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggg                                                                                                                      <  2:1646316/122‑1 (MQ=255)
                                                       tCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTggg                                                                                                                      >  1:1646316/1‑122 (MQ=255)
                                                        cgcgGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGttt                                                                                         >  2:162114/1‑150 (MQ=255)
                                                           ggCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGcc                                                                                      >  2:1537627/1‑150 (MQ=255)
                                                             ccGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATACGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCaa                                                                                    <  1:25147/150‑1 (MQ=255)
                                                              cGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAg                                                                                   <  1:692194/150‑1 (MQ=255)
                                                                  gCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCGCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCAt                                                                               >  1:321425/1‑150 (MQ=255)
                                                                         tGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCTACCGGTGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAg                                                                        <  1:172997/150‑1 (MQ=255)
                                                                         tGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAg                                                                        >  2:1499761/1‑150 (MQ=255)
                                                                         tGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGGGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAg                                                                        <  1:399952/150‑1 (MQ=255)
                                                                           gTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCtga                                                                                                 >  1:1838060/1‑123 (MQ=255)
                                                                                     aCTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAgg                                                                       <  1:2210252/139‑1 (MQ=255)
                                                                                     aCTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAgg                                                                       >  2:2210252/1‑139 (MQ=255)
                                                                                      cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAg                                                                        >  1:2401258/1‑137 (MQ=255)
                                                                                      cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAg                                                                        <  2:2401258/137‑1 (MQ=255)
                                                                                      cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCGCCGTCAAgg                                                           >  1:715233/1‑150 (MQ=255)
                                                                                      cTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAgg                                                           <  2:2226962/150‑1 (MQ=255)
                                                                                       tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGAc                                                                                                                                       <  1:2094778/73‑1 (MQ=255)
                                                                                       tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGAc                                                                                                                                       >  2:2094778/1‑73 (MQ=255)
                                                                                       tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                          >  1:422528/1‑150 (MQ=255)
                                                                                       tCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAggg                                                          >  1:1359493/1‑150 (MQ=255)
                                                                                         gggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGc                                                        >  1:982175/1‑150 (MQ=255)
                                                                                          ggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGcc                                                                                <  2:1655164/125‑2 (MQ=255)
                                                                                          ggCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCa                                                                                >  1:1655164/1‑125 (MQ=255)
                                                                                           gCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTGAAGGGGccc                                                      <  1:1680504/150‑1 (MQ=255)
                                                                                            cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAg                                                            >  1:1624388/1‑143 (MQ=255)
                                                                                            cTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGcccc                                                     >  2:1666338/1‑150 (MQ=255)
                                                                                               cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGCTGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                  <  2:2145208/150‑1 (MQ=255)
                                                                                               cATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                  <  1:162114/150‑1 (MQ=255)
                                                                                                 tCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAgg                                                                       >  2:886712/1‑127 (MQ=255)
                                                                                                 tCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAgg                                                                       <  1:886712/127‑1 (MQ=255)
                                                                                                 tCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGc                                                >  1:2152107/1‑150 (MQ=255)
                                                                                                  cTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGcc                                               <  2:2237312/150‑1 (MQ=255)
                                                                                                   tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                              >  1:1922842/1‑150 (MQ=255)
                                                                                                   tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                              >  1:51958/1‑150 (MQ=255)
                                                                                                   tACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCACCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCt                                              <  2:1799697/150‑1 (MQ=255)
                                                                                                     cAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGTCCCTGGGCCTTc                                            <  1:1537627/150‑1 (MQ=255)
                                                                                                      aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCcgcg                                                                                                                                            <  1:2214994/53‑1 (MQ=255)
                                                                                                      aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCcgcg                                                                                                                                            >  2:2214994/1‑53 (MQ=255)
                                                                                                      aaaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcg                                           >  2:1073357/1‑150 (MQ=255)
                                                                                                       aaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgc                                          <  1:1311353/150‑1 (MQ=255)
                                                                                                       aaGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTcgc                                          <  2:1140955/150‑1 (MQ=255)
                                                                                                                cgcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAATGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGccc                                 <  1:2290997/150‑1 (MQ=255)
                                                                                                                 gcTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCg                                >  1:2126103/1‑150 (MQ=255)
                                                                                                                    gATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGAt                             <  2:2368780/150‑1 (MQ=255)
                                                                                                                       tCGCGAAGCGGTGGATGCTATCGCGGCAAGCCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATggg                          <  1:506855/150‑1 (MQ=255)
                                                                                                                         gcgAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGc                        >  1:277001/1‑150 (MQ=255)
                                                                                                                            aaGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTg                     <  2:1922842/150‑1 (MQ=255)
                                                                                                                              gCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCa                   <  2:935101/150‑1 (MQ=255)
                                                                                                                               cGGTGGATGCTATCGCGGAAATGCAGAGGTGACCAGCCCGCCAATGTGGGGGGGCGAGTGGGTGAAGCTGATGAAGTTTGAAAAAAAAAAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCAt                  <  2:982175/150‑1 (MQ=255)
                                                                                                                                ggTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATg                 <  2:715233/150‑1 (MQ=255)
                                                                                                                                 gTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGcc                                                                    >  1:1213613/1‑98 (MQ=255)
                                                                                                                                 gTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGcc                                                                    <  2:1213613/98‑1 (MQ=255)
                                                                                                                                 gTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATgg                >  1:1628375/1‑150 (MQ=255)
                                                                                                                                       gtcaccgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                  <  1:1157301/105‑1 (MQ=255)
                                                                                                                                       gtcaccgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTgg                                                  >  2:1157301/6‑110 (MQ=255)
                                                                                                                                        tcaccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa            >  2:406461/5‑147 (MQ=255)
                                                                                                                                        tcaccgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTaaa            <  1:406461/143‑1 (MQ=255)
                                                                                                                                         tATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGa        >  1:1311625/1‑150 (MQ=255)
                                                                                                                                          accgcgGCATTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa       >  1:1878107/3‑150 (MQ=255)
                                                                                                                                          accgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATg                 >  1:1956792/3‑140 (MQ=255)
                                                                                                                                          accgcgGCATTGCCGCGATGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATg                 <  2:1956792/138‑1 (MQ=255)
                                                                                                                                          aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa       <  1:253891/150‑1 (MQ=255)
                                                                                                                                          aTCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGaa       <  2:2129621/150‑1 (MQ=255)
                                                                                                                                           tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAAt      >  2:5773/1‑150 (MQ=255)
                                                                                                                                           tCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAAt      >  2:28010/1‑150 (MQ=255)
                                                                                                                                              cgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa   <  1:28010/150‑1 (MQ=255)
                                                                                                                                              cgGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCa   <  2:2152107/150‑1 (MQ=255)
                                                                                                                                               ggCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAg  >  2:1423286/1‑150 (MQ=255)
                                                                                                                                                  aaTGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTa              >  1:588063/1‑135 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTGCCGATCATTGCCGCCGGGGGCATTACCGAAGGCCGCCATGCTGCCGAGAAGATCGCGGCCGGGGCGAGCCTGGTGCAGATATACTCGGGCTTCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGCCCCTGTGGGGGGGCGGGTGGGTGCAGCTGATGAGGTTTGCCAAGGCATAAAAAAGGGCCCCTTCAAGGGGCCCCTGGGCCTTCGCCCGCCGCCCGGATGGGGCGTGCATGGTAAATCGAATTCAG  >  NC_002947/2388148‑2388440

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