Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,907,574 T→C 8.8% D2G (GAC→GGC)  phoU ← negative regulator of PhoR/PhoB two‑component regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,907,5740TC 8.8% 649.5 / 41.6 216D2G (GAC→GGC) phoUnegative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base T (117/80);  new base C (5/14);  total (122/94)
Fisher's exact test for biased strand distribution p-value = 7.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

Only 100 of 216 total aligned reads displayed.

CGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCG  >  NC_000913/3907426‑3907720
                                                                                                                                                    |                                                                                                                                                  
cGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTc                                                                                                                                                   <  2:543853/150‑1 (MQ=255)
 gCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTcc                                                                                                                                                  <  2:2938714/150‑1 (MQ=255)
  ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCa                                                                                                                                                 >  1:1049349/1‑150 (MQ=255)
  ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCa                                                                                                                                                 >  1:2900814/1‑150 (MQ=255)
   cAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCAt                                                                                                                                                >  1:2090963/1‑150 (MQ=255)
      aTCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAAc                                                                                                                                             <  2:873092/150‑1 (MQ=255)
      aTCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAAc                                                                                                                                             <  2:347174/150‑1 (MQ=255)
      aTCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAAc                                                                                                                                             <  2:2422262/150‑1 (MQ=255)
         gCTGTCATGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCa                                                                                                                                          >  2:1567863/1‑150 (MQ=255)
          cTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCAc                                                                                                                                         >  2:2377333/1‑150 (MQ=255)
           tGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACt                                                                                                                                        >  1:778826/1‑150 (MQ=255)
           tGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACt                                                                                                                                        >  1:2641129/1‑150 (MQ=255)
             tCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTcc                                                                                                                                      <  2:1364257/150‑1 (MQ=255)
             tCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTcc                                                                                                                                      <  2:774851/150‑1 (MQ=255)
               cTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTg                                                                                                                                    >  1:1907926/1‑150 (MQ=255)
                tGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGa                                                                                                                                   <  1:1599005/150‑1 (MQ=255)
                tGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGa                                                                                                                                   >  1:3030190/1‑150 (MQ=255)
                 ggTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGaa                                                                                                                                  <  2:321641/150‑1 (MQ=255)
                      tGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAAc                                                                                                                             <  2:2740948/150‑1 (MQ=255)
                      tGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAAc                                                                                                                             <  2:2189872/150‑1 (MQ=255)
                      tGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAAc                                                                                                                             <  1:2614943/150‑1 (MQ=255)
                       gCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAAcc                                                                                                                            >  2:628234/1‑150 (MQ=255)
                           cgcgGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTaa                                                                                                                        <  1:2535957/150‑1 (MQ=255)
                            gcgGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAAc                                                                                                                       >  2:2874401/1‑150 (MQ=255)
                            gcgGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAAc                                                                                                                       >  2:1192057/1‑150 (MQ=255)
                            gcgGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAAc                                                                                                                       >  2:1618984/1‑150 (MQ=255)
                                  aTTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAATTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCgg                                                                                                                 >  2:1299483/1‑150 (MQ=255)
                                       aTCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATg                                                                                                            >  2:1128604/1‑150 (MQ=255)
                                        tCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGt                                                                                                           >  1:113216/1‑150 (MQ=255)
                                         cAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTa                                                                                                          <  2:439243/150‑1 (MQ=255)
                                         cAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTa                                                                                                          <  2:493088/150‑1 (MQ=255)
                                         cAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTa                                                                                                          <  2:1735797/150‑1 (MQ=255)
                                          aGAACGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGGTGAACTGGCCGGACATATGTTTCGTAAGATTGAGACGGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGGGagggag                                                                                                         >  2:1250227/1‑145 (MQ=255)
                                          aGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAg                                                                                                         >  2:1517760/1‑150 (MQ=255)
                                          aGAAAGCTGCTGCTCCACCATGCCGCACATGGACATCACCTGCGTGCGGAAACTATCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGAACGGTGAGGTat                                                                                                         >  2:1391368/1‑149 (MQ=255)
                                           gAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGt                                                                                                        >  2:2921995/1‑150 (MQ=255)
                                           gAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGt                                                                                                        >  1:1352272/1‑150 (MQ=255)
                                           gAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGt                                                                                                        >  1:1876327/1‑150 (MQ=255)
                                            aaaGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtc                                                                                                       >  1:1849019/1‑150 (MQ=255)
                                                   ctgctCCACCCGGCCGCCCCTGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGACCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtt                                                                                                <  2:2885354/150‑1 (MQ=255)
                                                   ctgctCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtt                                                                                                >  1:1137703/1‑150 (MQ=255)
                                                   ctgctCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtt                                                                                                >  1:2198105/1‑150 (MQ=255)
                                                     gctCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttg                                                                                              >  1:2011135/1‑150 (MQ=255)
                                                       tCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttgtt                                                                                            >  1:1019363/1‑150 (MQ=255)
                                                        cccccaGGCCGCCCATGGTCATCACCTGCGTGCGGATCCTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTAAGAAGATTGAGACGGTCCATAACGCACTCCGGAATCAACCGAAACGACCGGTGATGTAGTCTTCTgtttgttt                                                                                           <  1:2686116/147‑1 (MQ=255)
                                                           ccttGCCGCCATTGGTAATCACCTGCGTGCGGATCCTTTCCAGTTCGGCTTTGAACTGGCCGGAAATTTGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttctt                                                                                        <  1:1253490/147‑1 (MQ=255)
                                                            caTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCCTAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttcttc                                                                                       >  1:2556881/1‑150 (MQ=255)
                                                            caTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttcttc                                                                                       >  1:355175/1‑150 (MQ=255)
                                                            caTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttcttc                                                                                       >  1:1766346/1‑150 (MQ=255)
                                                               gccgccCTTGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCt                                                                                    <  1:2346591/150‑1 (MQ=255)
                                                               gccgccCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCt                                                                                    <  1:2346146/150‑1 (MQ=255)
                                                               gccgccCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCt                                                                                    <  2:2854897/150‑1 (MQ=255)
                                                                ccgccCCTGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTg                                                                                   <  1:1306525/150‑1 (MQ=255)
                                                                ccgccCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTGTATTAAGATTGAGACTGTCCATAACGCACTCCTGACTCACCCGTAACGACAGGGGAGGTCGTCTTCTGTTTGTTTCTTCGCTg                                                                                   >  2:2360779/1‑150 (MQ=255)
                                                                     cATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTg                                                                              <  2:2034325/150‑1 (MQ=255)
                                                                      aTGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTgg                                                                             >  1:1224313/1‑150 (MQ=255)
                                                                          tcatcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGaa                                                                         >  1:2452025/1‑150 (MQ=255)
                                                                             tcaCCTGCGTGCGTATAATTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGATAATGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAg                                                                      <  1:1592600/150‑1 (MQ=255)
                                                                             tcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGAGGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAg                                                                      >  2:76171/1‑150 (MQ=255)
                                                                              caCCTGGGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTAGTAAGTTTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGa                                                                     <  1:2015835/150‑1 (MQ=255)
                                                                              caCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGa                                                                     >  1:531869/1‑150 (MQ=255)
                                                                               aCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGAt                                                                    >  2:2749126/1‑150 (MQ=255)
                                                                               aCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGAt                                                                    >  2:1201707/1‑150 (MQ=255)
                                                                                     gtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCg                                                              >  1:1452274/1‑150 (MQ=255)
                                                                                       gcgGATACTTTCCAGNTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCgtg                                                            >  2:303811/1‑150 (MQ=255)
                                                                                        cgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCgtgt                                                           >  1:1104558/1‑150 (MQ=255)
                                                                                         ggATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGtt                                                          >  2:581363/1‑150 (MQ=255)
                                                                                          gATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATACCGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTg                                                         <  2:2846533/150‑1 (MQ=255)
                                                                                          gATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTg                                                         <  2:2590445/150‑1 (MQ=255)
                                                                                            tACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCt                                                       >  2:2044914/1‑150 (MQ=255)
                                                                                            tACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCt                                                       >  1:589928/1‑150 (MQ=255)
                                                                                            tACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCAGAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCt                                                       <  1:2874401/150‑1 (MQ=255)
                                                                                             aCTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTg                                                      <  2:2064677/150‑1 (MQ=255)
                                                                                             aCTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTg                                                      <  2:2064702/150‑1 (MQ=255)
                                                                                               tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGGGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGaa                                                    <  2:94148/150‑1 (MQ=255)
                                                                                               tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGaa                                                    >  2:725600/1‑150 (MQ=255)
                                                                                                ttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAAc                                                   >  1:406485/1‑150 (MQ=255)
                                                                                                    aGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaa                                               <  1:2025118/150‑1 (MQ=255)
                                                                                                    aGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaa                                               <  2:3068420/150‑1 (MQ=255)
                                                                                                     gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaat                                              <  1:1679351/150‑1 (MQ=255)
                                                                                                     gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaat                                              >  2:1869684/1‑150 (MQ=255)
                                                                                                     gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaat                                              <  1:85544/150‑1 (MQ=255)
                                                                                                         ggCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTcaatcaat                                          <  2:268775/150‑1 (MQ=255)
                                                                                                           cGTTGAACTGGCCGAAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTc                                        <  1:2251636/150‑1 (MQ=255)
                                                                                                             ttGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGc                                      <  2:1126109/150‑1 (MQ=255)
                                                                                                                aaCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCa                                   <  2:113225/150‑1 (MQ=255)
                                                                                                                  cTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAgg                                 >  1:1928798/1‑150 (MQ=255)
                                                                                                                   tGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGt                                >  2:2441786/1‑150 (MQ=255)
                                                                                                                               tGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCg                    <  1:2052101/150‑1 (MQ=255)
                                                                                                                                ttttATTATGATTGAGACGGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTTGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATGGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgt                   <  1:666143/149‑1 (MQ=255)
                                                                                                                                gTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgt                   <  2:1814418/150‑1 (MQ=255)
                                                                                                                                gTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgt                   >  2:715551/1‑150 (MQ=255)
                                                                                                                                gTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgt                   >  1:2774016/1‑150 (MQ=255)
                                                                                                                                 tttATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgtg                  >  1:2835905/1‑150 (MQ=255)
                                                                                                                                        aGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTAGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCgg           >  1:1931140/1‑150 (MQ=255)
                                                                                                                                         gATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGa          >  1:2113479/1‑150 (MQ=255)
                                                                                                                                          aTTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaa         >  2:2534730/1‑150 (MQ=255)
                                                                                                                                            tGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaaca       >  2:459452/1‑150 (MQ=255)
                                                                                                                                               gaCTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAAcg    <  2:1679473/150‑1 (MQ=255)
                                                                                                                                                 cTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAAcgcg  <  2:1359502/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCG  >  NC_000913/3907426‑3907720

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