Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,279,858 A→C 6.5% intergenic (+30/‑122) ghxP → / → yjcE guanine/hypoxanthine transporter GhxP/putative transporter YjcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,279,8580AC 6.5% 938.0 / 40.4 289intergenic (+30/‑122)ghxP/yjcEguanine/hypoxanthine transporter GhxP/putative transporter YjcE
Reads supporting (aligned to +/- strand):  ref base A (138/131);  new base C (6/13);  total (144/145)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

GCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATA  >  NC_000913/4279710‑4280006
                                                                                                                                                    |                                                                                                                                                    
    cgcgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGc                                                                                                                                                 <  1:951024/150‑1 (MQ=255)
    cgcgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGc                                                                                                                                                 >  1:1824614/1‑150 (MQ=255)
         ctctGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGccc                                                                                                                                            >  1:573650/1‑150 (MQ=255)
         ctctGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGccc                                                                                                                                            >  2:1329961/1‑150 (MQ=255)
               tGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttc                                                                                                                                      <  2:1931070/150‑2 (MQ=255)
               tGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttc                                                                                                                                      <  1:1802050/150‑2 (MQ=255)
                gATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttca                                                                                                                                     <  1:645484/150‑3 (MQ=255)
                 tttGGTCGTGTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGGCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttcag                                                                                                                                    <  2:200285/149‑4 (MQ=255)
                 aTTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttcag                                                                                                                                    <  1:1459678/150‑4 (MQ=255)
                  ttGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAtt                                                                                                                                   <  1:1601594/150‑1 (MQ=255)
                  ttGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttcaga                                                                                                                                   <  1:1927974/150‑5 (MQ=255)
                   tGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAttt                                                                                                                                  >  1:1139602/1‑150 (MQ=255)
                    ggTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAtttt                                                                                                                                 >  1:1166688/1‑150 (MQ=255)
                    ggTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTAtttt                                                                                                                                 >  1:735781/1‑150 (MQ=255)
                    ggTCGTCTGGGTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCGGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGCAGCCATTGCCACCCGttttcagaat                                                                                                                                 <  2:695787/150‑7 (MQ=255)
                        gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTcc                                                                                                                             <  1:1332997/150‑1 (MQ=255)
                        gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTcc                                                                                                                             >  1:888756/1‑150 (MQ=255)
                        gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTcc                                                                                                                             <  1:1125697/150‑1 (MQ=255)
                          cTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGc                                                                                                                           <  1:630384/150‑1 (MQ=255)
                          cTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGc                                                                                                                           <  1:84715/150‑1 (MQ=255)
                                tCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtg                                                                                                                     <  1:508303/150‑1 (MQ=255)
                                tCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtg                                                                                                                     >  2:1071769/1‑150 (MQ=255)
                                 ccGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTtc                                                                                                                                <  2:1555819/138‑1 (MQ=255)
                                 ccGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgt                                                                                                                    >  1:758246/1‑150 (MQ=255)
                                  cGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtg                                                                                                                   >  1:359567/1‑150 (MQ=255)
                                  cGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtg                                                                                                                   >  1:1890300/1‑150 (MQ=255)
                                   ggCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgt                                                                                                                  <  1:276559/150‑1 (MQ=255)
                                    gCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                                                 <  2:958561/150‑1 (MQ=255)
                                    gCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                                                 <  2:819259/150‑1 (MQ=255)
                                       aaTGGCGCAAGTTGAACATAGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGGGGCAAAGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGc                                                                                                              >  2:40817/1‑150 (MQ=255)
                                        aTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCt                                                                                                             >  2:1811464/1‑150 (MQ=255)
                                            cgcAAGTTGAACATCGGTACGGTCGTTATCGCGGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTg                                                                                                         <  2:575751/150‑1 (MQ=255)
                                                 gTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGata                                                                                                    >  2:550352/1‑150 (MQ=255)
                                                    gAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTc                                                                                                 >  1:312692/1‑150 (MQ=255)
                                                    gAACATCGGTACGGGCGTTATCGCCGTGGCGCTGGTGACCTGCTATGCGGGTGGCGGGGCTAGCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAAAGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTc                                                                                                 <  2:1396404/150‑1 (MQ=255)
                                                       cATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGa                                                                                              >  2:775605/1‑150 (MQ=255)
                                                         tCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaaa                                                                                            <  1:1083127/150‑1 (MQ=255)
                                                          cGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaaaa                                                                                           >  2:1004748/1‑150 (MQ=255)
                                                          cGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaaaa                                                                                           >  2:1304590/1‑150 (MQ=255)
                                                          cGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaaaa                                                                                           <  2:540594/150‑1 (MQ=255)
                                                           ggTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGaaaaa                                                                                          >  1:372841/1‑150 (MQ=255)
                                                              aCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAAt                                                                                       <  2:1139602/150‑1 (MQ=255)
                                                              aCGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAAt                                                                                       <  2:1061673/150‑1 (MQ=255)
                                                                ggTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATga                                                                                     <  2:1843438/150‑1 (MQ=255)
                                                                 ggcGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATgag                                                                                    <  2:463167/148‑1 (MQ=255)
                                                                 gTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATgag                                                                                    >  2:1363808/1‑150 (MQ=255)
                                                                   cGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGaa                                                                                  >  2:725506/1‑150 (MQ=255)
                                                                      tATCGCCGTGGGGCGGGGGACCTGCTATGCGGGTGGCTGGGCGATCAAATCATTCATCGGTTTCTGAACCCGGGGGGCTATGGGTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGACAAtgc                                                                               <  2:1627280/150‑3 (MQ=255)
                                                                      tATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTc                                                                               >  2:3341/1‑150 (MQ=255)
                                                                       aTCGCGGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCa                                                                              <  2:878431/150‑1 (MQ=255)
                                                                        tCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAg                                                                             <  1:1708115/150‑1 (MQ=255)
                                                                         cGCCGGGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAAGCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAgg                                                                            <  2:1505717/150‑1 (MQ=255)
                                                                          gCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGc                                                                           >  1:1539828/1‑150 (MQ=255)
                                                                            cGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCAt                                                                         <  2:1774901/150‑1 (MQ=255)
                                                                                       tGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtct                                                              >  2:297334/1‑150 (MQ=255)
                                                                                         aCCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCa                                                            <  1:1041616/150‑1 (MQ=255)
                                                                                         aCCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCa                                                            <  1:1650988/150‑1 (MQ=255)
                                                                                          ccTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAg                                                           <  1:1103143/150‑1 (MQ=255)
                                                                                           cTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGTCTGCCCGTTTTTATTTTCTCGGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                                          <  1:1848372/150‑1 (MQ=255)
                                                                                           cTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                                          >  1:1513436/1‑150 (MQ=255)
                                                                                           cTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                                          >  1:1255630/1‑150 (MQ=255)
                                                                                            ttCTATGCGGTTGGATGGTCTATATAATCTTTCATCTGATTCTGAAAACGGGTTGCAATGGCTGCCAGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCTATATTCTTAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGAt                                                         <  1:810575/150‑1 (MQ=255)
                                                                                             tCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATc                                                        <  1:672809/150‑1 (MQ=255)
                                                                                              cTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCCGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCCCACATTGACTGAGCTTTTGCGATATTCTGAAAAAAAAGAGAATTAAGGCATAACGTCTTTCTAACATCa                                                       >  2:1740902/1‑150 (MQ=255)
                                                                                                aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATa                                                     >  2:912688/1‑150 (MQ=255)
                                                                                                 tGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaa                                                    >  1:1995839/1‑150 (MQ=255)
                                                                                                 tGCGGGTGGCTGGGCTATCTAATCTTTAATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaa                                                    >  1:1310488/1‑150 (MQ=255)
                                                                                                    gggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATaaaaa                                                 <  1:628939/150‑1 (MQ=255)
                                                                                                     ggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATc                                                        >  2:277362/1‑142 (MQ=255)
                                                                                                          cTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCg                                           >  2:1176679/1‑150 (MQ=255)
                                                                                                           tGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGc                                          <  1:550352/150‑1 (MQ=255)
                                                                                                           tGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGc                                          <  1:1920150/150‑1 (MQ=255)
                                                                                                             ggCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCaa                                        <  1:1223408/150‑1 (MQ=255)
                                                                                                             ggCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCaa                                        <  2:1659348/150‑1 (MQ=255)
                                                                                                               cTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAAc                                      <  2:1909368/150‑1 (MQ=255)
                                                                                                                 aTCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAg                                    >  1:612905/1‑150 (MQ=255)
                                                                                                                  tCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAgg                                   <  1:453637/150‑1 (MQ=255)
                                                                                                                   cTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAggg                                  >  1:1001322/1‑150 (MQ=255)
                                                                                                                     aaTCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGaa                                >  2:56751/1‑150 (MQ=255)
                                                                                                                       tCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCATGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACg                              <  2:1123662/150‑1 (MQ=255)
                                                                                                                         tttCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCa                            >  1:1811231/1‑150 (MQ=255)
                                                                                                                            cATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATgg                         <  2:982202/150‑1 (MQ=255)
                                                                                                                             attggATTCTTTAAAAGGGTTGCAATTGCTTCCCGTTTTTAATTTTTCCCGAAATTTTTTTATCTTTTTCGATATTCTTAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGa                        <  1:696350/147‑1 (MQ=255)
                                                                                                                             aTCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGa                        <  1:1065039/150‑1 (MQ=255)
                                                                                                                                ggATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAAc                                      <  1:591351/133‑1 (MQ=255)
                                                                                                                                ggATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAt                     >  1:15750/1‑150 (MQ=255)
                                                                                                                                ggATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAt                     >  2:1081709/1‑150 (MQ=255)
                                                                                                                                 gATTCTGAAAACGGGTGTCAATTGCTGCCCGTTTTTATTTTCTCCGCACATAGTTTGAGCTTTTGCGATATTCTGAAAAAAATGAGATTTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtc                    <  1:1281042/150‑1 (MQ=255)
                                                                                                                                 gATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTATCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtg                    >  2:1707907/1‑149 (MQ=255)
                                                                                                                                  aTTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAATAACATCGCAAACAGGGAACGCATGGAAAtct                   <  2:359567/150‑1 (MQ=255)
                                                                                                                                    tCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttc                 <  2:592157/150‑1 (MQ=255)
                                                                                                                                     cTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttct                >  2:970802/1‑150 (MQ=255)
                                                                                                                                     cTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttct                >  2:514296/1‑150 (MQ=255)
                                                                                                                                      tGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGGGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATcttctt               >  1:978281/1‑150 (MQ=255)
                                                                                                                                       gAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAAtcttct                <  1:1758408/148‑1 (MQ=255)
                                                                                                                                        aaaaCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCa             <  2:584131/150‑1 (MQ=255)
                                                                                                                                          aaCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGc                                          <  2:654430/119‑1 (MQ=255)
                                                                                                                                            cGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCAt         <  2:626728/150‑1 (MQ=255)
                                                                                                                                                  gCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGAt   <  2:1255630/150‑1 (MQ=255)
                                                                                                                                                   cAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATa  >  1:212274/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTTCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCATGGAAATCTTCTTCACCATACTGATA  >  NC_000913/4279710‑4280006

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