Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 431,111 A→T 5.8% intergenic (+159/+18) yajD → / ← tsx HNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913431,1110AT 5.8% 371.3 / 7.7 127intergenic (+159/+18)yajD/tsxHNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K
Reads supporting (aligned to +/- strand):  ref base A (59/60);  new base T (4/4);  total (63/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

AAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAG  >  NC_000913/430962‑431259
                                                                                                                                                     |                                                                                                                                                    
aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATaaa                                                                                                                                                      >  2:1801002/1‑150 (MQ=255)
 aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAAt                                                                                                                                                     >  1:1166336/1‑150 (MQ=255)
   aTTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATcc                                                                                                                                                   >  1:1160618/1‑150 (MQ=255)
      gCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCgg                                                                                                                                                <  2:159903/150‑1 (MQ=255)
       ccTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGc                                                                                                                                               <  1:1236443/150‑1 (MQ=255)
       ccTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGc                                                                                                                                               >  1:740492/1‑150 (MQ=255)
       ccTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGc                                                                                                                                               <  2:1474721/150‑1 (MQ=255)
          gATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtt                                                                                                                                            <  1:306354/150‑1 (MQ=255)
             gcgcTACGCTTATCAGGCCAACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATACATCCCGGCATGTTc                                                                                                                                         >  1:1165920/1‑150 (MQ=255)
               gcaaCGCATATCAGGCCTACGTTATTTCAGCAAAATATTGAAATTACGTGCTTATGTTGGCCGGACAAATCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttctga                                                                                                                                       <  1:332389/147‑5 (MQ=25)
                cTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATa                                                                                                                                      >  2:1403631/1‑150 (MQ=255)
                 tACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATaa                                                                                                                                     <  1:35439/150‑1 (MQ=255)
                    gCTTCTCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGAATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCa                                                                                                                                  <  2:1155077/150‑1 (MQ=255)
                      ttATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGa                                                                                                                                >  2:739967/1‑150 (MQ=255)
                        aTCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAg                                                                                                                              <  2:298378/150‑1 (MQ=255)
                          cAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGtt                                                                                                                            >  2:1424422/1‑150 (MQ=255)
                          cAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGtt                                                                                                                            >  2:417006/1‑150 (MQ=255)
                          cAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGAATATTCATAATCAGAAGtt                                                                                                                            >  2:1503585/1‑150 (MQ=255)
                              ccTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGTAATAAATCCCGGCATTTTCATAATCAGAAGTTGTaa                                                                                                                        <  1:1601168/150‑1 (MQ=255)
                              ccTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTaa                                                                                                                        <  1:1499993/150‑1 (MQ=255)
                              ccTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTaa                                                                                                                        >  2:138793/1‑150 (MQ=255)
                               gtACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGACGCATCCCGCATGAACAAAGCACACGTTGTTTACAATAAGAAATTCCGGGAATAAATCCCGGGAATTTAATAATCAGAAGTTGTAAc                                                                                                                       >  2:1361220/2‑150 (MQ=255)
                                  cGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACcta                                                                                                                    >  2:913611/1‑150 (MQ=255)
                                   gTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACctac                                                                                                                   <  1:1398389/150‑1 (MQ=255)
                                    ttATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACctact                                                                                                                  >  1:1373249/1‑150 (MQ=255)
                                    ttATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACctact                                                                                                                  >  2:752385/1‑150 (MQ=255)
                                        ttCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCa                                                                                                              <  2:582979/150‑1 (MQ=255)
                                         tCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAg                                                                                                             <  1:956270/150‑1 (MQ=255)
                                          cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAgg                                                                                                            >  1:875011/1‑150 (MQ=255)
                                          cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAgg                                                                                                            >  1:1208849/1‑150 (MQ=255)
                                            gCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTa                                                                                                          <  1:1571424/150‑1 (MQ=255)
                                               atctatTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAacc                                                                                                       <  2:1127675/147‑1 (MQ=255)
                                               atatatTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGATGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTACATAATCAGAAGTTGAAACCTACTACCAGGTAaca                                                                                                       >  1:1767688/1‑149 (MQ=255)
                                                       aaTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCAGAATCAGAAGTTGTAACCTACTACCAGGTAGCCACCCAAGc                                                                                               >  1:275292/1‑150 (MQ=255)
                                                          ttACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCgg                                                                                            <  2:1216629/150‑1 (MQ=255)
                                                          ttACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCgg                                                                                            <  1:86778/150‑1 (MQ=255)
                                                            aCGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTa                                                                                          <  1:1751870/150‑1 (MQ=255)
                                                              gTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTaga                                                                                        <  2:539845/150‑1 (MQ=255)
                                                                gCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGc                                                                                      <  2:1265076/150‑1 (MQ=255)
                                                                gCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGc                                                                                      <  2:432398/150‑1 (MQ=255)
                                                                gCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACAACCAGGTAACCACCCCAGCCGGTAGAGc                                                                                      <  1:1029011/150‑1 (MQ=255)
                                                                       ttgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGtt                                                                               >  1:1067377/1‑150 (MQ=255)
                                                                        tgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTg                                                                              >  2:1786797/1‑150 (MQ=255)
                                                                         ggCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGa                                                                             >  2:1582015/1‑150 (MQ=255)
                                                                         ggCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGa                                                                             >  2:24967/1‑150 (MQ=255)
                                                                         ggCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGa                                                                             >  2:1237440/1‑150 (MQ=255)
                                                                         ggCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGa                                                                             <  2:292988/150‑1 (MQ=255)
                                                                          gCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                            >  2:1465892/1‑150 (MQ=255)
                                                                           ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAg                                                                           >  1:1546722/1‑150 (MQ=255)
                                                                              gACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACAAGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCATGTAACCACCACATCAGGTAGAGCGAACGTTGAActtt                                                                        >  2:702565/1‑149 (MQ=255)
                                                                                 aaaGCGTTTACGCCGCATCCGGCATGCACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCACGGCATTTTCATAATCAGAAGTTGTAACCTACTACAAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccg                                                                     >  1:227582/1‑150 (MQ=255)
                                                                                  aaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgt                                                                    >  1:546857/1‑150 (MQ=255)
                                                                                   aGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgtt                                                                   <  1:1080333/150‑1 (MQ=255)
                                                                                      gTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgcc                                                                <  1:1234446/150‑1 (MQ=255)
                                                                                      gTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgcc                                                                >  2:1676382/1‑150 (MQ=255)
                                                                                        ttACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGa                                                              <  1:297542/150‑1 (MQ=255)
                                                                                          aCGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAag                                                            >  1:167783/1‑150 (MQ=255)
                                                                                             ccgcATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGTGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttc                                                         <  2:468130/150‑1 (MQ=255)
                                                                                             ccgcATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttc                                                         <  2:1320442/150‑1 (MQ=255)
                                                                                             ccgcATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttc                                                         <  1:382459/150‑1 (MQ=255)
                                                                                                cATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagt                                                      <  2:631307/150‑1 (MQ=255)
                                                                                                    cGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTTCCGAAGTTCAGTTCTg                                                  >  2:1667079/1‑150 (MQ=255)
                                                                                                    cGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTg                                                  <  2:881987/150‑1 (MQ=255)
                                                                                                    cGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTg                                                  <  2:1524190/150‑1 (MQ=255)
                                                                                                      gCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCa                                                <  2:1050500/150‑1 (MQ=255)
                                                                                                        aTGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtc                                              <  2:1690071/150‑1 (MQ=255)
                                                                                                          gAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtcgt                                            >  2:1665878/1‑150 (MQ=255)
                                                                                                            aCAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtcgtcg                                          <  1:721269/150‑1 (MQ=255)
                                                                                                            aCAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtcgtcg                                          >  1:641471/1‑150 (MQ=255)
                                                                                                              aaaGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGtt                                        <  2:1547744/150‑1 (MQ=255)
                                                                                                               atgCACACGTTGTAAACAATACGAAATGCCGGGAATAAATCCCGGCATTTTCATAAACAGACGTTGTACCCAACTCCCCGTTAACCACCCCACCCGGTTGAGCGACCGTTGAAGTTGCCTTGCCGAAGTTCAGTTCTGCATCGTCGTTcc                                      <  2:631272/148‑1 (MQ=255)
                                                                                                                aGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTcc                                      <  2:1087460/150‑1 (MQ=255)
                                                                                                                   acacGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACt                                   <  2:1474315/150‑1 (MQ=255)
                                                                                                                      cGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGac                                >  1:12643/1‑150 (MQ=255)
                                                                                                                          gttAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGaccacc                            >  1:386764/1‑150 (MQ=255)
                                                                                                                          gttAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGaccacc                            >  2:1320826/1‑150 (MQ=255)
                                                                                                                             aaCAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACcgtc                         >  1:337307/1‑150 (MQ=255)
                                                                                                                              aCAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACcgtcg                        <  1:695233/150‑1 (MQ=255)
                                                                                                                               cAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACcgtcgt                       >  2:1339385/1‑150 (MQ=255)
                                                                                                                                aaTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCCGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTg                      >  2:1055380/1‑150 (MQ=255)
                                                                                                                                aaTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTg                      <  1:835214/150‑1 (MQ=255)
                                                                                                                                 tatgaaaaTGCCGTGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTACCCTACTACCATGTAACCACCACAGCCGGTATAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCTGTTCTGCATCGTCGTTCCACTGCACACCGTCGTGc                     >  2:679809/6‑150 (MQ=255)
                                                                                                                                  tCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGcc                    >  2:1135741/1‑150 (MQ=255)
                                                                                                                                   cAGAAAAGCCGGGAATAAATCCCGGCATATTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTTGAGAGAACGTTGAAGATGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGcct                   >  2:906818/1‑149 (MQ=255)
                                                                                                                                    aGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAg                  >  2:416470/1‑150 (MQ=255)
                                                                                                                                    aGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAg                  >  2:1627511/1‑150 (MQ=255)
                                                                                                                                    aGAAATGCAGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAg                  <  1:1125083/150‑1 (MQ=255)
                                                                                                                                      aaaTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTa                <  1:200975/150‑1 (MQ=255)
                                                                                                                                        aTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAAc              >  1:1544202/1‑150 (MQ=255)
                                                                                                                                        aTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAAc              >  1:400087/1‑150 (MQ=255)
                                                                                                                                        aTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAAc              >  1:1272063/1‑150 (MQ=255)
                                                                                                                                        aTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAAc              >  2:1825867/1‑150 (MQ=255)
                                                                                                                                          gCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGa            >  2:537267/1‑150 (MQ=255)
                                                                                                                                          gCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGa            >  1:241720/1‑150 (MQ=255)
                                                                                                                                            cGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAActgga          >  2:1816285/1‑146 (MQ=255)
                                                                                                                                            cGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGc          <  1:1801002/150‑1 (MQ=255)
                                                                                                                                              ggAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTa        <  1:1668703/150‑1 (MQ=255)
                                                                                                                                               gAATAAATCCCGGCATTTTCTTAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTCAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTAc       <  1:949820/150‑1 (MQ=255)
                                                                                                                                               gAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTAc       >  2:1515523/1‑150 (MQ=255)
                                                                                                                                                    aaTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCTAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACag  <  2:157427/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGACAG  >  NC_000913/430962‑431259

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