Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 287,925 A→C 5.0% Q379H (CAA→CAC)  yagH → CP4‑6 prophage; putative xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913287,9250AC 5.0% 2253.4 / 21.1 753Q379H (CAA→CAC) yagHCP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (373/322);  new base C (2/55);  total (375/378)
Fisher's exact test for biased strand distribution p-value = 3.19e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 766 total aligned reads displayed.

AACTGGCGGGACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGG  >  NC_000913/287776‑288073
                                                                                                                                                     |                                                                                                                                                    
aaCTGGCGGGACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCaa                                                                                                                                                      <  2:8699262/150‑1 (MQ=255)
 aCTGGCGGGACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAAt                                                                                                                                                     <  1:8747780/150‑1 (MQ=255)
         gACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTg                                                                                                                                             >  2:1478124/1‑150 (MQ=255)
         gACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTg                                                                                                                                             <  2:8140380/150‑1 (MQ=255)
          aCGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCAGCGCAATCCGATCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTACTTACTGCTCTATGGCAACGACTCGCTCAATACGACCTTCACCCAATCGACCGTgg                                                                                                                                            >  2:7363085/1‑150 (MQ=255)
            gATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGgcg                                                                                                                                          >  1:4724467/1‑150 (MQ=255)
                tCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGcgcc                                                                                                                                      <  2:8422642/150‑1 (MQ=255)
                     gCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTgg                                                                                                                                 <  2:5665667/150‑1 (MQ=255)
                     gCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTgg                                                                                                                                 >  1:502929/1‑150 (MQ=255)
                     gCCAGTTCGCTGGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTgg                                                                                                                                 <  2:7857794/150‑1 (MQ=255)
                       cAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgca                                                                                                                               <  2:2546455/150‑1 (MQ=255)
                        aGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgcag                                                                                                                              >  2:9089757/1‑150 (MQ=255)
                          ttCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGgcagca                                                                                                                            >  1:4318904/1‑150 (MQ=255)
                              cTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTc                                                                                                                        >  1:4146480/1‑150 (MQ=255)
                               ttGACCCGGAGCTGCAGACCCTGCGCATTCGGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTAGGGCTATTTGGCAACGACTCGCTCAATTCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCg                                                                                                                       <  1:6771249/150‑1 (MQ=255)
                               ttGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCg                                                                                                                       >  2:6533489/1‑150 (MQ=255)
                               ttGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCg                                                                                                                       <  2:163688/150‑1 (MQ=255)
                               ttGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCg                                                                                                                       <  2:1070710/150‑1 (MQ=255)
                                 gACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGcc                                                                                                                     >  2:8200226/1‑150 (MQ=255)
                                    ccGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTc                                                                                                                  >  1:7963445/1‑150 (MQ=255)
                                    ccGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTc                                                                                                                  >  2:3948201/1‑150 (MQ=255)
                                    ccGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTc                                                                                                                  >  2:3358782/1‑150 (MQ=255)
                                       gAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCgg                                                                                                               >  2:1239487/1‑150 (MQ=255)
                                        aGCTGCAGACCCTGCGAATTCGTTTCGACGACACCCTCGGCTCGCACACCGCGCGCCCGGGCTTCTTCCGGCTCTAGGGAAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCggg                                                                                                              <  1:3822242/150‑1 (MQ=255)
                                              aGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTAGGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGaaa                                                                                                        <  2:6453185/150‑1 (MQ=255)
                                                   cTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCgg                                                                                                   >  2:313622/1‑150 (MQ=255)
                                                   cTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCgg                                                                                                   >  1:5375205/1‑150 (MQ=255)
                                                     gcgcATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGAt                                                                                                 >  1:1076388/1‑150 (MQ=255)
                                                         aTTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTCCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAg                                                                                             <  2:279476/150‑1 (MQ=255)
                                                          ttCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGt                                                                                            >  2:7683338/1‑150 (MQ=255)
                                                                tCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGc                                                                                      <  1:6671299/150‑1 (MQ=255)
                                                                  gacgacACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCg                                                                                    <  1:2493625/150‑1 (MQ=255)
                                                                  gacgacACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCg                                                                                    <  1:703998/150‑1 (MQ=255)
                                                                     gacACCCTCGGCGCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCCATTCTACCTTCACCCAATCGACAGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTGCTCGCCGGTg                                                                                 >  1:9121537/1‑150 (MQ=255)
                                                                       cacCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCa                                                                               <  1:1227126/150‑1 (MQ=255)
                                                                        acCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCAc                                                                              >  1:5824489/1‑150 (MQ=255)
                                                                          ccTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACtt                                                                            >  1:7456098/1‑150 (MQ=255)
                                                                          ccTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCAAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACtt                                                                            >  1:6295457/1‑150 (MQ=255)
                                                                               gctcgctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCcagc                                                                       >  2:9213673/1‑150 (MQ=255)
                                                                                ctcgctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCcagca                                                                      <  1:3313947/150‑1 (MQ=255)
                                                                                 tcgctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCcagcag                                                                     >  1:4619958/1‑150 (MQ=255)
                                                                                   gctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCagag                                                                   >  1:6058016/1‑150 (MQ=255)
                                                                                       aCCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGcgcg                                                               <  1:3711749/150‑1 (MQ=255)
                                                                                       aCCGCGCGCCCGGGCTACTTACGGCACAATGGCAACGACTCGCTCAATTCGAACTACAACCAATCGAACTTGGCGCGCAGCTGGCAGAACTTCGCCTACCGGGCAGAAACGCTGATGGAGTACTCGCACGTGCACTTCCAGCAGAGcgcg                                                               >  1:3159964/1‑150 (MQ=255)
                                                                                        ccGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgg                                                              <  1:8921922/150‑1 (MQ=255)
                                                                                        ccGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgg                                                              <  1:5165247/150‑1 (MQ=255)
                                                                                          gcgcgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgggg                                                            <  2:3726171/150‑1 (MQ=255)
                                                                                           cgcgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGc                                                           >  2:6968376/1‑150 (MQ=255)
                                                                                              gcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                        <  2:481915/150‑1 (MQ=255)
                                                                                                ccGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAcc                                                      >  1:2230878/1‑150 (MQ=255)
                                                                                                 cGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                     <  2:6016745/150‑1 (MQ=255)
                                                                                                  gggCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTg                                                    >  2:4308274/1‑150 (MQ=255)
                                                                                                   ggCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGc                                                   <  1:732106/150‑1 (MQ=255)
                                                                                                     cttcttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGcta                                                 <  1:1734175/150‑1 (MQ=255)
                                                                                                     cttcttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGcta                                                 <  1:1000402/150‑1 (MQ=255)
                                                                                                     cttcttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGcta                                                 <  1:7218364/150‑1 (MQ=255)
                                                                                                        cttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctacta                                              >  1:7824733/1‑150 (MQ=255)
                                                                                                        cttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctacta                                              >  1:957378/1‑150 (MQ=255)
                                                                                                        cttACGGCTCTATGGCAACGACTCGCTCAAATCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctacta                                              >  2:176094/1‑150 (MQ=255)
                                                                                                            cGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaac                                          <  2:887340/150‑1 (MQ=255)
                                                                                                            cGGCTCTATGGCAACGACTCGCTCAAGTCGACCTTCACCCCCTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctacttgcgg                                          <  2:6224265/150‑6 (MQ=255)
                                                                                                             ggCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaaca                                         <  1:1497969/150‑1 (MQ=255)
                                                                                                             ggCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaaca                                         <  2:6802850/150‑1 (MQ=255)
                                                                                                             ggCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTacaaca                                         >  2:7146135/1‑150 (MQ=255)
                                                                                                               ctctATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGc                                       >  1:159436/1‑150 (MQ=255)
                                                                                                               ctctATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGc                                       >  2:2995684/1‑150 (MQ=255)
                                                                                                               ctctATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGc                                       >  1:8524591/1‑150 (MQ=255)
                                                                                                               ctctATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGc                                       >  1:1383573/1‑150 (MQ=255)
                                                                                                                tctATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCa                                      >  2:7504738/1‑150 (MQ=255)
                                                                                                                tctATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCa                                      >  1:2898528/1‑150 (MQ=255)
                                                                                                                  tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCaaa                                    <  1:1894204/150‑1 (MQ=255)
                                                                                                                       cAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTg                               >  2:7953336/1‑150 (MQ=255)
                                                                                                                       cAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTg                               <  1:3365189/150‑1 (MQ=255)
                                                                                                                       cAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTg                               <  2:4624504/150‑1 (MQ=255)
                                                                                                                        aaCGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTgg                              <  1:6460463/150‑1 (MQ=255)
                                                                                                                        aaCGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTgg                              <  2:2293929/150‑1 (MQ=255)
                                                                                                                         aCGACTCGCTCAATTCGACCTTCACCCATTAGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGa                             <  2:3543568/150‑1 (MQ=255)
                                                                                                                          cGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAg                            >  1:8073560/1‑150 (MQ=255)
                                                                                                                           gACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAACACTGGAGc                           >  2:5165028/1‑150 (MQ=255)
                                                                                                                           gACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGc                           >  1:4941147/1‑150 (MQ=255)
                                                                                                                           gACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGc                           >  2:4979807/1‑150 (MQ=255)
                                                                                                                                gCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTg                      >  2:4764646/1‑150 (MQ=255)
                                                                                                                                 cTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGc                     >  1:3515615/1‑150 (MQ=255)
                                                                                                                                      ttCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtgt                <  2:5081335/150‑1 (MQ=255)
                                                                                                                                      ttCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtgt                <  1:6146463/150‑1 (MQ=255)
                                                                                                                                      ttCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCGTtgt                >  2:4131015/1‑150 (MQ=255)
                                                                                                                                       tCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtgtg               <  1:3017994/150‑1 (MQ=255)
                                                                                                                                       tCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtgtg               <  1:9185468/150‑1 (MQ=255)
                                                                                                                                        cGACCTTCACCCACTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTgg              <  2:3888929/150‑1 (MQ=255)
                                                                                                                                            cTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTa          >  2:4426638/1‑150 (MQ=255)
                                                                                                                                            cTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTa          >  1:395224/1‑150 (MQ=255)
                                                                                                                                            cTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTa          >  1:7656289/1‑150 (MQ=255)
                                                                                                                                              tCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACg        >  2:5856420/1‑150 (MQ=255)
                                                                                                                                               cACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACga       >  1:2637131/1‑150 (MQ=255)
                                                                                                                                                 cccAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgagg     >  1:5732246/1‑150 (MQ=255)
                                                                                                                                                   cAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgaggag   >  2:5382227/1‑150 (MQ=255)
                                                                                                                                                   cAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgaggag   <  2:2370682/150‑1 (MQ=255)
                                                                                                                                                    caTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGaggagg  <  2:2019894/149‑1 (MQ=255)
                                                                                                                                                    aaTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGaggagg  <  1:8259314/150‑1 (MQ=255)
                                                                                                                                                    aaTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCGTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGaggagg  <  1:5571040/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACTGGCGGGACGATTTCGACGCCAGTTCGCTTGACCCGGAGCTGCAGACCCTGCGCATTCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGG  >  NC_000913/287776‑288073

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: ‑