Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,990,927 A→G 8.6% intergenic (+35/+27) tyrP → / ← yecA tyrosine transporter/UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,990,9270AG 8.6% 263.0 / 13.3 95intergenic (+35/+27)tyrP/yecAtyrosine transporter/UPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base A (39/47);  new base G (4/4);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CCGGGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCC  >  NC_000913/1990784‑1991071
                                                                                                                                               |                                                                                                                                                
ccGGGTCAAAGGTGGTCGTCCGGCGCTGGTGTTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACg                                                                                                                                            >  1:1078567/1‑150 (MQ=255)
 cGGGGCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGGATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGc                                                                                                                                           <  2:899345/150‑1 (MQ=255)
  gggTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCt                                                                                                                                          <  2:249303/150‑1 (MQ=255)
   ggTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGGTTCTCTGTGGTATTGCTGTGATTGGCGTGCAAATTTTGATTGCGGCAGGGTTGTTACCAGATGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTg                                                                                                                                         >  1:1012712/1‑150 (MQ=255)
     gcAAAGGTGGTCGGCCGGCGCGGGTGGTGGTGTTTCTCTGGGGTATTGCGGGGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGCAGTGGGGGGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGac                                                                                                                                       <  2:1755950/149‑1 (MQ=255)
      cAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGaca                                                                                                                                      >  1:1051616/1‑150 (MQ=255)
       aaaGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGAATTTTTTGATTGCGGCAGGGTTTTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGacac                                                                                                                                     <  1:1200560/150‑1 (MQ=255)
         aGGTGGTCGTCCGGCGCTGGTGGTGGTGGTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACAccc                                                                                                                                   <  2:1334183/150‑1 (MQ=255)
              gtcgtcCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAtaa                                                                                                                              <  1:2017171/150‑3 (MQ=255)
               tcgtcCTGCGCTGGATGTTGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTTCAATTTTTGATTGCTTCAGTGTTGTTACCAGAAGTGGAGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAtaag                                                                                                                             <  1:34250/150‑4 (MQ=255)
               tcgtcCGGCGCTGGTGGGGGGGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCTGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAtaag                                                                                                                             <  1:1473123/150‑4 (MQ=255)
                cgtcCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTgg                                                                                                                            <  2:615146/150‑1 (MQ=255)
                cgtcCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCCGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGTCAGCGTTAATTCTGGCACCCAAtaagg                                                                                                                            <  1:205128/150‑5 (MQ=38)
                  tcCGGCGCTGGTGGGTGTGTTTCTCTGTGGTATTGCTGGGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGcc                                                                                                                          <  2:1316411/150‑1 (MQ=255)
                  gccGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGcc                                                                                                                          <  1:15995/149‑1 (MQ=255)
                        gcTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTaa                                                                                                                    <  1:1076629/150‑1 (MQ=255)
                        gcTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTaa                                                                                                                    <  1:1965718/150‑1 (MQ=255)
                         cTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAAt                                                                                                                   <  1:37747/150‑1 (MQ=255)
                          tggtggtggTGTTTCTATGTGGTATTGCTGTGCTGGGCGTTCAATTTTTGATTGCGGGAGGGGTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGTGAGACAGAATTAACGCTGACACCCAATTTGGCCTCTTAATg                                                                                                                  <  1:920130/150‑1 (MQ=255)
                           ggtggtggtGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgc                                                                                                                 >  2:1407563/1‑150 (MQ=255)
                            gtggttgtgTTTCTCTGTGGTATTGCGGTGATTGGCGTGGAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGTTGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgca                                                                                                                <  2:1438734/150‑1 (MQ=255)
                            gtggtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATATGGCCTCTAAATgca                                                                                                                >  2:436062/1‑150 (MQ=255)
                            gtggtggtgTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATATTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATgca                                                                                                                >  1:1528805/1‑150 (MQ=255)
                              ggtggtgTTTCTCTGTGGTATTGCTGTGCTTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGCTATCCTCAAATTCCTTATTGGTTGCCAGAATTAACGATGACACCAATTTTGGCCTCTTAATgcagg                                                                                                              >  2:1083081/1‑150 (MQ=255)
                                    gtTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGCGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAAATTGGCCTCTTAATGCAGGCAGCAc                                                                                                        >  2:239645/1‑150 (MQ=255)
                                      ttCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGAGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTg                                                                                                      <  1:391917/150‑1 (MQ=255)
                                       tctcTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGTGTGATCAGATAGCCTCAAATTCCTTATTTGGTGCCAGAATTAAAGCTGACACACAATTTGGCCTCTTAATGCAGGCAGCACTga                                                                                                     >  1:227168/1‑149 (MQ=255)
                                       tctcTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGc                                                                                                     <  1:1906931/150‑1 (MQ=255)
                                       tctcTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGc                                                                                                     <  1:1573508/150‑1 (MQ=255)
                                        ctctGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCt                                                                                                    >  2:620183/1‑150 (MQ=255)
                                        ctctGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACAAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGACAGATTTAACGCTGACACCAAATTTGGCCTCTTAATGCAGGCAGCAATgat                                                                                                    >  1:1283055/1‑148 (MQ=255)
                                          ctGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTa                                                                                                  <  2:5330/150‑1 (MQ=255)
                                           gggggTATGGCGGTGATTGGCGTGCAATGTTTGATTGCGGCAGGGGGGGTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTaa                                                                                                 <  1:1285102/147‑1 (MQ=255)
                                              ggTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAAtt                                                                                              <  2:405673/150‑1 (MQ=255)
                                               gTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAAttt                                                                                             >  1:1170099/1‑150 (MQ=255)
                                                tctagCTGTGATTGGCGTGCAATTTTTGATAGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATACCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTc                                                                                            <  2:1022368/146‑1 (MQ=255)
                                                tATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCATGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTc                                                                                            <  1:1789674/150‑1 (MQ=255)
                                                gaTTGCTCTGATTGGCGTGCAATGTTTGATTGCGGCAGGGGGGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTc                                                                                            <  1:388602/149‑1 (MQ=255)
                                                 aTTGCTGTGATTGGCGTGCAATTTTTGTTTGCGGCAGGGGTGTTACCAGAAGTGGGGCGAGCCTATACCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCt                                                                                           <  2:284070/150‑1 (MQ=255)
                                                 aTTGCTGGGATTGGCGTGCAATTTTTGATTGCGGCAGGGGTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGGGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCt                                                                                           <  1:1270690/150‑1 (MQ=255)
                                                  ttGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCtt                                                                                          <  2:1051616/150‑1 (MQ=255)
                                                   tGCTGTGATGGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTa                                                                                         <  1:1591970/150‑1 (MQ=255)
                                                    gCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCATGCAGCACTGCTTAAATTTCTTAc                                                                                        >  1:1646140/1‑150 (MQ=255)
                                                     cTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTAcc                                                                                       >  1:616744/1‑150 (MQ=255)
                                                     cTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTAcc                                                                                       <  2:1704549/150‑1 (MQ=255)
                                                        tgATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACt                                                                                    >  1:690460/1‑150 (MQ=255)
                                                           ttGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTAcc                                                                                 <  2:32389/150‑1 (MQ=255)
                                                            tGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCg                                                                                <  1:20755/150‑1 (MQ=255)
                                                                gTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACg                                                                            >  1:579675/1‑150 (MQ=255)
                                                                gTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACg                                                                            >  1:2016886/1‑150 (MQ=255)
                                                                   cAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGc                                                                         >  2:1614463/1‑150 (MQ=255)
                                                                    aaTTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCATTTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCa                                                                        >  2:975296/1‑150 (MQ=255)
                                                                     aTTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCaa                                                                       >  2:184987/1‑150 (MQ=255)
                                                                         ttGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGAt                                                                   >  2:541985/1‑150 (MQ=255)
                                                                                cGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCATGGATCGTTACg                                                            >  1:499519/1‑150 (MQ=255)
                                                                                  gCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGcc                                                          >  1:905956/1‑150 (MQ=255)
                                                                                        ttgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGtt                                                    <  1:279884/150‑1 (MQ=255)
                                                                                        ttgttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGtt                                                    <  1:460117/150‑1 (MQ=255)
                                                                                          gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                  <  2:1709846/150‑1 (MQ=255)
                                                                                          gttACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACACAATTTGGCCTCTTAATGCAGTCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTtc                                                  <  2:812595/150‑1 (MQ=255)
                                                                                             aCCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTctct                                               >  2:723204/1‑150 (MQ=255)
                                                                                              ccAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCtctt                                              >  1:827642/1‑150 (MQ=255)
                                                                                               cAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCtctta                                             >  1:116581/1‑149 (MQ=255)
                                                                                                 gtaGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCAGTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTg                                           <  2:271496/148‑1 (MQ=255)
                                                                                                    gTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGAGCGTTACGCCCCGGTTTCTCTTCTGCtt                                        <  2:692156/150‑1 (MQ=255)
                                                                                                     tGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCttt                                       >  1:779340/1‑150 (MQ=255)
                                                                                                        ggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGAt                                    >  2:1615818/1‑150 (MQ=255)
                                                                                                        ggTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAAAGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTAACACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTATCTTCTGCTTTGAt                                    >  1:1499889/1‑150 (MQ=255)
                                                                                                           gATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCgg                                 >  1:640368/1‑150 (MQ=255)
                                                                                                              cAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTg                              <  1:788877/150‑1 (MQ=255)
                                                                                                              cAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTg                              <  1:1771797/150‑1 (MQ=255)
                                                                                                               aGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGc                             <  1:1963930/150‑1 (MQ=255)
                                                                                                                gATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCt                            >  2:497159/1‑150 (MQ=255)
                                                                                                                 aTAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCTCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg                           <  2:609475/150‑1 (MQ=255)
                                                                                                                 aTAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTg                           <  2:1139098/150‑1 (MQ=255)
                                                                                                                    gCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAAc                        <  2:1078567/150‑1 (MQ=255)
                                                                                                                     ccTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACa                       >  2:578094/1‑150 (MQ=255)
                                                                                                                     ccTCAAATTCCTTATTGGGAGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACa                       <  1:194294/150‑1 (MQ=255)
                                                                                                                       tCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGc                     >  1:456961/1‑150 (MQ=255)
                                                                                                                           aTTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTACTGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCtttt                 <  1:406261/150‑1 (MQ=255)
                                                                                                                            ggccaaaTTGGGTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGTCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTAGATCGGTTGCTTAACAGCttttt                >  1:304034/7‑150 (MQ=37)
                                                                                                                             tCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTc               >  1:1786768/1‑150 (MQ=255)
                                                                                                                                aaaTTGGGTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTg            >  2:459465/3‑150 (MQ=255)
                                                                                                                                 tcttGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGc           <  1:1407563/148‑1 (MQ=255)
                                                                                                                                  aTTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCg          <  1:1836558/150‑1 (MQ=255)
                                                                                                                                   ttGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCgg         <  1:483230/150‑1 (MQ=255)
                                                                                                                                    tGGGTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGa        >  2:1611795/1‑150 (MQ=255)
                                                                                                                                    tGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTTAATTTCTTACCACTTCCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGCACAGCTTTTTCCTGCGGa        >  1:251283/1‑150 (MQ=255)
                                                                                                                                    tGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGAggt        >  2:149886/1‑149 (MQ=255)
                                                                                                                                      ggTGTCAGCGTTAATTCTGGCACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATg      >  2:535784/1‑150 (MQ=255)
                                                                                                                                        tGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgcg    >  1:1975968/1‑150 (MQ=255)
                                                                                                                                         gCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATgcgc   >  2:786019/1‑150 (MQ=255)
                                                                                                                                          ccAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGTGcc  >  2:583783/1‑150 (MQ=255)
                                                                                                                                          ccAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGcc  >  2:1943484/1‑150 (MQ=255)
                                                                                                                                          ccAGAAATAACGCCGACGCCCAATAGGGCCTCCTAATGCACGCAGCACTGCTTAACTTTCATAACACTACCGCACGGGCAAGGAGCGTTACGCCCCGGTTTCTCTTCAGCTTTGATCGGTAGCTGAACAGCTTTTTCCTGCGGATGCGcc  <  1:294665/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CCGGGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCC  >  NC_000913/1990784‑1991071

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: