Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,504 T→A 5.8% intergenic (+21/‑111) tnaB → / → mdtL tryptophan transporter of low affinity/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,5040TA 5.8% 409.7 / 9.2 135intergenic (+21/‑111)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (65/62);  new base A (4/4);  total (69/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 137 total aligned reads displayed.

AGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGT  >  NC_000913/3891356‑3891649
                                                                                                                                                    |                                                                                                                                                 
agTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCAtt                                                                                                                                                  >  2:1180163/1‑150 (MQ=255)
    tCCCAATGAGATGGTGTCGGGCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTTTTACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTATCGTGTTTCCTAATTTTGGCTAAATCCTTCTAGAAGCCAGCCATTCGCt                                                                                                                                              <  2:1055913/150‑1 (MQ=255)
     cccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAATAAGCCAGCCATTCGCTg                                                                                                                                             >  2:1524577/1‑150 (MQ=255)
      ccAATCAGCTGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                                                            <  2:534437/150‑1 (MQ=255)
      ccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                                                            <  2:478759/150‑1 (MQ=255)
      ccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTCTCATTCTCGTAGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTTCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                                                            >  1:1159476/1‑150 (MQ=255)
       cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                           <  1:1612912/150‑1 (MQ=255)
       cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                           >  2:683572/1‑150 (MQ=255)
        aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATggctggct                                                                                                                                          <  1:1745722/150‑1 (MQ=255)
            cgATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttctt                                                                                                                                      <  2:873603/149‑1 (MQ=255)
             gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                                                                     <  2:874964/150‑1 (MQ=255)
             gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGATTCTTg                                                                                                                                     >  1:36650/1‑150 (MQ=255)
             gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTGCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTg                                                                                                                                     <  1:1062153/150‑1 (MQ=255)
              aTGTTCACGGTCGGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                                                                    <  2:78476/150‑1 (MQ=255)
               tGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGcc                                                                                                                                   <  2:1800203/150‑1 (MQ=255)
                gTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGAATGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaag                                                                                                                                  <  1:1631526/150‑4 (MQ=255)
                gTTCACGGGCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTGCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                                                                  <  2:1106975/150‑1 (MQ=255)
                  tCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTCCCTAAATTTGGCTAAATCCTACAAGAAGACAGCCATTCGCTGGCTTCTTGCctcg                                                                                                                                >  2:885234/1‑149 (MQ=255)
                  tCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctct                                                                                                                                >  1:199620/1‑150 (MQ=255)
                  tCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTAAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctct                                                                                                                                >  1:14164/1‑150 (MQ=255)
                      ggTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                                                                            <  1:1226364/150‑1 (MQ=255)
                      ggTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                                                                            <  2:827761/150‑1 (MQ=255)
                        tCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                                                                          >  1:1656508/1‑150 (MQ=255)
                           ggggCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                                                                       <  1:828246/150‑1 (MQ=255)
                               cggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACtt                                                                                                                   >  1:852177/1‑150 (MQ=255)
                                ggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTa                                                                                                                  >  1:1662778/1‑150 (MQ=255)
                                ggcAATCTTATTCCGGCGATTCTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGCTTCGGCAACTTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCAAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTa                                                                                                                  >  2:535515/1‑150 (MQ=255)
                                 gcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTAAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                                                                 >  1:843420/1‑150 (MQ=255)
                                 gcAATCTTAATCCGGCGATTGTCATTCTCTTTGGTATAACCGTGAATTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                                                                 >  1:1766039/1‑150 (MQ=255)
                                  cAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                                                                <  2:576555/150‑1 (MQ=255)
                                    aTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTc                                                                                                              <  2:604153/150‑1 (MQ=255)
                                      cTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCa                                                                                                            <  1:860793/150‑1 (MQ=255)
                                         aTTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAAt                                                                                                         >  2:430756/1‑150 (MQ=255)
                                             cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCa                                                                                                     >  1:1523626/1‑150 (MQ=255)
                                                cGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACt                                                                                                  >  1:673929/1‑150 (MQ=255)
                                                 gATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                                                                 <  1:965563/150‑1 (MQ=255)
                                                  aTTGTCATTCTCTTTGGTATACCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCg                                                                                                >  1:300398/1‑150 (MQ=255)
                                                    tGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAGATCACTTATGTCCAAATGGCAACTCGcc                                                                                              <  2:377570/150‑1 (MQ=255)
                                                    tGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                                                                              >  1:1841278/1‑150 (MQ=255)
                                                     gTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATATGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTGCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                                                                             >  2:1343655/1‑150 (MQ=255)
                                                       cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGC‑‑CTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtc                                                                                         <  2:293649/150‑1 (MQ=255)
                                                        aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                          <  1:466284/150‑1 (MQ=255)
                                                         ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtc                                                                                         <  2:624973/150‑1 (MQ=255)
                                                          tctcTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                                                                        <  1:1090365/150‑1 (MQ=255)
                                                          tctcTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                                                                        >  1:139826/1‑150 (MQ=255)
                                                          tctcTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                                                                        >  2:1027613/1‑150 (MQ=255)
                                                           ctctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcat                                                                                       >  1:1819810/1‑148 (MQ=255)
                                                                 ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                 >  1:1270340/1‑150 (MQ=255)
                                                                 ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                 >  1:1006910/1‑150 (MQ=255)
                                                                 ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGCGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                 <  2:1608707/150‑1 (MQ=255)
                                                                  gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcac                                                                                <  1:683572/150‑1 (MQ=255)
                                                                    ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACATAAATTTGGATAAATCCTTCAAGAAGCCAGACATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATACTCCTTcacca                                                                              >  1:135335/1‑150 (MQ=255)
                                                                      aaCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACg                                                                            >  1:754680/1‑150 (MQ=255)
                                                                        ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                                                                          >  1:573747/1‑150 (MQ=255)
                                                                         cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTAt                                                                         >  2:1046078/1‑150 (MQ=255)
                                                                             aTTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCtt                                                                     >  2:979377/1‑150 (MQ=255)
                                                                              ttttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGTCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCttt                                                                    <  2:1727347/150‑1 (MQ=255)
                                                                                 tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                                                                 <  1:982286/150‑1 (MQ=255)
                                                                                 tgtgCGGGGTCGGCTACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTTGCCTGATCCTCCTTCACCACGTATGCTTTGCg                                                                 <  1:1014049/150‑1 (MQ=255)
                                                                                   tgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                                                               <  2:1860247/150‑1 (MQ=255)
                                                                                   ggCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                                                               <  2:961792/149‑1 (MQ=255)
                                                                                    gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                                                              <  1:1831627/150‑1 (MQ=255)
                                                                                    gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                                                              <  1:1583292/150‑1 (MQ=255)
                                                                                        gTTCGGCAACGTCTTTAACGTGTTTCCTAAATTTGGCTAAATCCTTCAAGAAGGCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCtt                                                          <  2:25442/150‑1 (MQ=255)
                                                                                           cGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACt                                                       <  2:373583/150‑1 (MQ=255)
                                                                                           cGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGAGTCACCTTACt                                                       >  2:299678/1‑150 (MQ=255)
                                                                                             gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTAt                                                     <  2:1624913/150‑1 (MQ=255)
                                                                                                aagTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAg                                                  <  2:1674383/148‑1 (MQ=255)
                                                                                                     tttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCCCCACGTATGCTTTGCGTCCCCTTACTATCAGGACGc                                             >  1:438505/1‑150 (MQ=255)
                                                                                                       tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCtt                                           <  1:1088300/150‑1 (MQ=255)
                                                                                                        aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTAACCTTACTATCAGGACGCtta                                          >  1:716727/1‑149 (MQ=255)
                                                                                                        aaCGTGTTACCTAAATGTGGCTAAATCCTTCAAGAAGCCTGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCttt                                          <  1:953480/150‑1 (MQ=255)
                                                                                                                aCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATg                                  >  2:693754/1‑150 (MQ=255)
                                                                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGt                                 >  1:517882/1‑150 (MQ=255)
                                                                                                                    aaaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAATTCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTccc                              <  1:55624/150‑1 (MQ=255)
                                                                                                                    aaaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTccc                              <  2:284563/150‑1 (MQ=255)
                                                                                                                        ttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCtt                          <  1:955713/150‑1 (MQ=255)
                                                                                                                          ggCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGAtttt                        >  2:62007/1‑150 (MQ=255)
                                                                                                                           gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCttttt                       <  1:1017527/150‑1 (MQ=255)
                                                                                                                           gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCGCTCAGGAAATCACTTAAGTCCAAATGGCAACTCGCATGATCCTCGTTCACGACCTATGCTTTGAGTCACAGTACTATCAGGACGCTTGAGCCCTTGTCCCGCttttt                       >  2:206610/1‑150 (MQ=255)
                                                                                                                            cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTg                      <  1:625025/150‑1 (MQ=255)
                                                                                                                            cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCAAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTg                      >  1:876362/1‑150 (MQ=255)
                                                                                                                             tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCATGACGCTTTAGCCCATGTCCCGCTTTTTGa                     >  1:1208036/1‑150 (MQ=255)
                                                                                                                                  tcTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGGCTCTCAGGAAATAACTTATGTCCAAATGGCAACTCGCCTGATACTCCTTCACCACGTTTGCTTTTCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGt                <  1:1524577/149‑1 (MQ=255)
                                                                                                                                    ggcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTAGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAg              >  2:1001118/3‑150 (MQ=255)
                                                                                                                                       aagaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGTCGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGatt           >  2:1119466/1‑150 (MQ=255)
                                                                                                                                       aagaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTGTTTGATTTGTAGttt           >  2:151015/1‑150 (MQ=255)
                                                                                                                                         gaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTg         >  1:1414343/1‑150 (MQ=255)
                                                                                                                                         gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGttttt         >  2:1505745/1‑149 (MQ=255)
                                                                                                                                          aagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCCCCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGc        >  1:1733567/1‑150 (MQ=255)
                                                                                                                                           agccagccATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGcc       <  1:479579/150‑1 (MQ=255)
                                                                                                                                            gccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGccc      <  2:1656508/150‑1 (MQ=255)
                                                                                                                                            gccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCGTTTTGATTTGTAGTTTTGccc      >  1:1653192/1‑150 (MQ=255)
                                                                                                                                             ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCt     <  1:805803/150‑1 (MQ=255)
                                                                                                                                              cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTg    <  1:1345763/150‑1 (MQ=255)
                                                                                                                                              cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTg    >  1:557662/1‑150 (MQ=255)
                                                                                                                                              cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTg    >  1:1100874/1‑150 (MQ=255)
                                                                                                                                              cagccaTTCGCTGGCTTCTTGACTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCATACTTCACCAAGTATGCTTTGAGTCACCTTAATATCAGGACGCTTTAGCACATGTCCAGCTTTTTGATTTGTAGTTTTGCCCTg    >  1:1150981/1‑150 (MQ=255)
                                                                                                                                               agccaTTCGCTGGCTTCTTGCCTCTCCGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCCCCACGTATGCGTTGCGTCACCTTACTATCACGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTgg   >  2:520179/1‑150 (MQ=255)
                                                                                                                                                gccaTTCGCTGGCTTCTTGCCTCTAAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGt  <  2:14164/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGT  >  NC_000913/3891356‑3891649

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: