Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,990,922 C→T 8.6% intergenic (+30/+32) 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,9220CT 8.6% 274.1 / 9.8 94intergenic (+30/+32)tyrP/yecAtyrosine transporter/UPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base C (37/47);  new base T (4/4);  total (42/52)
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

AGGCGGGTTACCGGGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCC  >  NC_000913/1990774‑1991071
                                                                                                                                                    |                                                                                                                                                     
aGGCGGGGTACCGGGTCAAAGGGGGTCGTCCGGCGCTGGTGGTGGGGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGcc                                                                                                                                                      <  1:581499/150‑1 (MQ=255)
 ggCGGGTTACCGGGTCAAAGGGGGGCGGCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGGTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCa                                                                                                                                                     <  1:1019433/150‑1 (MQ=255)
  gCGGGTTACCGGGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAg                                                                                                                                                    <  2:15689/150‑1 (MQ=255)
          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:1965718/150‑1 (MQ=255)
                                  gcTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTaa                                                                                                                    <  1:1076629/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                                                                                       <  2:1704549/150‑1 (MQ=255)
                                                               cTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTAcc                                                                                       >  1:616744/1‑150 (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)
                                                                                                                                                    |                                                                                                                                                     
AGGCGGGTTACCGGGTCAAAGGTGGTCGTCCGGCGCTGGTGGTGGTGTTTCTCTGTGGTATTGCTGTGATTGGCGTGCAATTTTTGATTGCGGCAGGGTTGTTACCAGAAGTGGGGTGATCAGATAGCCTCAAATTCCTTATTGGGTGCCAGAATTAACGCTGACACCCAATTTGGCCTCTTAATGCAGGCAGCACTGCTTAAATTTCTTACCACTACCGCACGGGCAAGGATCGTTACGCCCCGGTTTCTCTTCTGCTTTGATCGGTTGCTGAACAGCTTTTTCCTGCGGATGCGCC  >  NC_000913/1990774‑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: