Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,270,142 G→A 8.4% intergenic (+15/‑96) aphA → / → yjbQ acid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,270,1420GA 8.4% 359.4 / 13.7 119intergenic (+15/‑96)aphA/yjbQacid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (54/55);  new base A (6/4);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 121 total aligned reads displayed.

ATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTG  >  NC_000913/4269994‑4270291
                                                                                                                                                    |                                                                                                                                                     
aTAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                                                                                                      <  1:1366513/150‑1 (MQ=255)
 tAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTg                                                                                                                                                     >  2:39840/1‑150 (MQ=255)
 tAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTg                                                                                                                                                     >  2:224964/1‑150 (MQ=255)
  aaTGATATTACCCCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGa                                                                                                                                                    <  2:1108778/150‑1 (MQ=255)
   aTGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt                                                                                                                                                   >  2:191169/1‑150 (MQ=255)
     gATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGc                                                                                                                                                 <  1:149904/150‑1 (MQ=255)
     gATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGGTTGGTGAAGCGGTGCTCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGc                                                                                                                                                 >  1:256802/1‑150 (MQ=255)
       tatTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGGCAATTCAGACTACTGACAGAGCGGGAGCGCGTGATGctc                                                                                                                                               >  2:856494/1‑150 (MQ=255)
           aCCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTccc                                                                                                                                           <  2:566783/150‑1 (MQ=255)
            ccgccgCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCg                                                                                                                                          <  1:1423718/150‑1 (MQ=255)
            ccgccgCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                                                                          <  2:707644/150‑1 (MQ=255)
            ccgccgCACGCGATGTCGGCGCTCGTGGAATCCCCAATCTCCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGGGCGTATGGTGAAGAGGTGATCGTCAATTCAGAATACTTACAGAGCGGGAGAGCGTGATGCTCTccct                                                                                                                                          >  2:1396042/1‑149 (MQ=255)
              gccgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGaa                                                                                                                                        >  2:341768/1‑150 (MQ=255)
                cgcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATTCTGACAGAGCGGGAGAGCATCACGCTCTCCCgctct                                                                                                                                      <  2:428294/150‑5 (MQ=255)
                 gcACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAAATCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGc                                                                                                                                     >  2:1612392/1‑150 (MQ=255)
                   aCGCGCTGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTg                                                                                                                                   <  2:1470475/150‑1 (MQ=255)
                    cgcgATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCTGAATACTGACAGAGCGGGATAGCGTGATGCTCTCCCGAATGCTGt                                                                                                                                  <  1:1393583/150‑1 (MQ=255)
                      cgATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttt                                                                                                                                >  2:1457026/1‑150 (MQ=255)
                       gATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttt                                                                                                                               >  1:1463602/1‑150 (MQ=255)
                            cGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTCCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                                                                          <  1:191113/150‑1 (MQ=255)
                              gcgcgcGTGGTATCCGCATTCTGCGCGCCTCCAACTCAACCTACAAACCCTTGCCACAAGCGGGGGCGTTTGGTGAAGAGGTGATCGTCAAATCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATc                                                                                                                        <  1:867313/150‑1 (MQ=255)
                              gcgcTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTCTTTTTAATc                                                                                                                        <  2:1314978/150‑1 (MQ=255)
                                gcTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCGTTGCCACAAGCGGGTGCGTTTGGGGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcac                                                                                                                      <  2:911780/150‑1 (MQ=255)
                                gcTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCATTTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcac                                                                                                                      >  2:859046/1‑150 (MQ=255)
                                 cTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcacc                                                                                                                     >  1:1493621/1‑149 (MQ=255)
                                     tGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCtt                                                                                                                 >  2:704400/1‑150 (MQ=255)
                                     tGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCtt                                                                                                                 >  2:1645827/1‑150 (MQ=255)
                                         aTCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                                                                                                             >  2:1237815/1‑150 (MQ=255)
                                          tCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGCTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCAGTTTTTTTAATCACACCTTTATcc                                                                                                            >  1:1645275/1‑150 (MQ=255)
                                           ccGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCGCACCTTTATCCt                                                                                                           >  2:729726/1‑150 (MQ=255)
                                           ccGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCt                                                                                                           >  1:1372946/1‑150 (MQ=255)
                                           ccGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCt                                                                                                           <  1:1454187/150‑1 (MQ=255)
                                           ccGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCt                                                                                                           >  2:912555/1‑150 (MQ=255)
                                              cATTCTGCGCGCCTCCAACTCTCACTACAACCCCTTGCCACAAGCAGGGGCGTTTGGGGCAGAGCTGCTCGTCTATTCAGAATACTGACAGAGCGGGAGAGCGGGATGCTCTCCCGAATGCTGTCTTTTTAATCACACCTTTAACCTTTc                                                                                                        <  2:1093382/150‑1 (MQ=255)
                                                     ggcgcCTCCAACTCTAACTACAACCCCTTGCCACACGCGGGTGCGGTTGGTGAAGTGGGGAGCGCCAATCCAGAATACTGACAGCGCGGGAGAGCGTGTGGCTCTCCGGAATGCTTTTTTTTTAATCACACCTTTATCCTTGCGCTGTCt                                                                                                 <  2:1150920/149‑1 (MQ=255)
                                                     cgcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCt                                                                                                 <  1:489861/150‑1 (MQ=255)
                                                     cgcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCt                                                                                                 >  1:1028164/1‑150 (MQ=255)
                                                         ccTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgct                                                                                             >  2:686695/1‑150 (MQ=255)
                                                          cTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCATAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctg                                                                                            <  1:1270145/150‑1 (MQ=255)
                                                               aCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGTTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAAc                                                                                       <  2:1277919/150‑1 (MQ=255)
                                                                ctctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACt                                                                                      >  2:992778/1‑150 (MQ=255)
                                                                 tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTg                                                                                     >  2:1117350/1‑150 (MQ=255)
                                                                    acctacAAACCCTTGCCACCAGCGGGTGCGTTTGGGGAAGAGGTGATCGTCAATTCAGAATACTCACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAtt                                                                                  <  1:306739/150‑1 (MQ=255)
                                                                    acctacAAAACCTTGCCTCAAGCGGGTGCGTTGGGTGAAGAGGGGATCGTCAATTCAGAATACTGACAGAGCGGGCGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGAtt                                                                                  <  2:62432/150‑1 (MQ=255)
                                                                     cctccAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTa                                                                                 <  1:1577243/150‑1 (MQ=255)
                                                                      ctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTaa                                                                                <  1:420806/150‑1 (MQ=255)
                                                                       tacAAACCGTTGCCACAAGCGGGTGGGTTTGGTGATGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                                                                               <  2:1460997/150‑1 (MQ=255)
                                                                       tacAAACCCTTGCCACAAGCGGGGGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                                                                               <  1:39840/150‑1 (MQ=255)
                                                                        acAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAaga                                                                              <  1:615245/150‑1 (MQ=255)
                                                                         cAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagag                                                                             <  1:192586/150‑1 (MQ=255)
                                                                          aaaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                                                                            <  1:96794/150‑1 (MQ=255)
                                                                           aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagagag                                                                           >  2:928544/1‑150 (MQ=255)
                                                                              ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAACCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGttt                                                                        <  2:878901/150‑1 (MQ=255)
                                                                              ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGGGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGttt                                                                        <  2:1516196/150‑1 (MQ=255)
                                                                               cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGtttt                                                                       >  1:1557077/1‑150 (MQ=255)
                                                                                ttGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTa                                                                      <  2:1649660/150‑1 (MQ=255)
                                                                                    cacaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATcga                                                                  >  2:1540762/1‑148 (MQ=255)
                                                                                    cacaAGCGGGTGCGTATGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTAATCATTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                                                                  >  1:68066/1‑150 (MQ=255)
                                                                                      caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgg                                                                <  1:1262858/150‑1 (MQ=255)
                                                                                       aaGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGa                                                               <  2:712564/150‑1 (MQ=255)
                                                                                            ggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagcagc                                                          >  1:825079/1‑150 (MQ=255)
                                                                                                 gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAtg                                                     >  1:450007/1‑150 (MQ=255)
                                                                                                 gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAtg                                                     >  1:629301/1‑150 (MQ=255)
                                                                                                  tttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAtgt                                                    <  1:217622/150‑1 (MQ=255)
                                                                                                   ttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAtgtg                                                   >  1:1504209/1‑150 (MQ=255)
                                                                                                    tGGTGCAGAGGGGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTgg                                                  <  1:1563128/150‑1 (MQ=255)
                                                                                                    tGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTgg                                                  <  1:1662574/150‑1 (MQ=255)
                                                                                                    tGGGGAAGAGGTGAGCGTCAATTCTGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTgg                                                  <  1:1570176/150‑1 (MQ=255)
                                                                                                        gAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATc                                              >  2:805792/1‑150 (MQ=255)
                                                                                                        gAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATACTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCCCATGTGGTATc                                              >  1:324102/1‑150 (MQ=255)
                                                                                                        gAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGAGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATc                                              >  1:1214044/1‑150 (MQ=255)
                                                                                                         aaGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTCTCCTTTCGCTGTCTGGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGGGGTATCa                                             <  1:284902/150‑1 (MQ=255)
                                                                                                            agGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCaaaa                                          <  1:341768/150‑1 (MQ=255)
                                                                                                                 aTCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGc                                     <  2:62684/150‑1 (MQ=255)
                                                                                                                  tCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCt                                    >  1:103654/1‑150 (MQ=255)
                                                                                                                    gTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCa                                  >  2:520497/1‑150 (MQ=255)
                                                                                                                     tCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCAc                                 <  1:984289/150‑1 (MQ=255)
                                                                                                                          tCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACCCCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTa                            >  1:1534486/1‑150 (MQ=255)
                                                                                                                            aGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAgc                          >  2:1201422/1‑150 (MQ=255)
                                                                                                                             gAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAgcg                         <  1:1645827/150‑1 (MQ=255)
                                                                                                                                  cTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGAAGATCACGCTTAGCGCCaaa                    >  2:1661246/1‑150 (MQ=255)
                                                                                                                                    cagcattcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtc                  >  1:666909/8‑150 (MQ=255)
                                                                                                                                      cAGATCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTCTCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtctc                <  2:276944/150‑1 (MQ=255)
                                                                                                                                      cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCAGCAAACTGATTAAGAGAGTTTTATAAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAAtctc                >  2:1540631/1‑150 (MQ=255)
                                                                                                                                       cattcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGTAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCg               >  1:640695/5‑150 (MQ=255)
                                                                                                                                       agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCg               <  2:1604591/150‑1 (MQ=255)
                                                                                                                                       agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCg               <  2:646296/150‑1 (MQ=255)
                                                                                                                                       agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCg               <  1:1371318/150‑1 (MQ=255)
                                                                                                                                        attcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGt              >  2:1505398/4‑150 (MQ=255)
                                                                                                                                        attcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGt              >  1:432359/4‑150 (MQ=255)
                                                                                                                                        attcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGt              >  1:1625586/4‑150 (MQ=255)
                                                                                                                                         agCGGGAGAGCGTGATGCTCGCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTTGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTg             <  2:845830/150‑1 (MQ=255)
                                                                                                                                           cGGGTGAGCGTGATGCTCTACCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGATAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTggg           <  1:224964/150‑1 (MQ=255)
                                                                                                                                             ggAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTATTCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTTTCAAAAGACGCTCCCGCTTTGCGCAAAATCTCGTGGGtt         >  2:340653/1‑150 (MQ=255)
                                                                                                                                             ggAGAGCGCGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGtt         <  2:1214044/150‑1 (MQ=255)
                                                                                                                                             ggAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGtt         >  2:1578004/1‑150 (MQ=255)
                                                                                                                                                  gCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATc    >  1:442019/1‑150 (MQ=255)
                                                                                                                                                  gCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATc    <  1:744261/150‑1 (MQ=255)
                                                                                                                                                    gTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTg  <  2:243131/150‑1 (MQ=255)
                                                                                                                                                    gTGATGCTCTCCCGAATGCTGAATTTTTAATCACACCATTATCCTTTCGCAGTCTTGCTGCAAACTGATTAAGAGACTTTTATAAAGGAAAAGCACATTTAGTATCAAAAGACACTCACTCTAAGCGACAAATCTCGTGCGATTCATCTg  >  2:273624/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATAATGATATTACCGCCGCACGCGATGTCGGCGCTCGTGGTATCCGCATTCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGTATCAAAAGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTG  >  NC_000913/4269994‑4270291

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