Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,654,143 T→A 9.5% noncoding (289/304 nt) ryfA → novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,654,1430TA 9.5% 442.4 / 30.4 137noncoding (289/304 nt)ryfAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (67/57);  new base A (5/8);  total (72/65)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

Only 100 of 138 total aligned reads displayed.

GCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGT  >  NC_000913/2653996‑2654289
                                                                                                                                                   |                                                                                                                                                  
gCACAATGAAGATATGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCttt                                                                                                                                                   >  1:794218/1‑149 (MQ=255)
 cacaATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTg                                                                                                                                                  <  2:17018/149‑1 (MQ=255)
    aaTGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggc                                                                                                                                               >  3:1038008/1‑149 (MQ=255)
     aTGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                                                              >  4:502567/1‑149 (MQ=255)
     aTGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                                                              >  2:575889/1‑149 (MQ=255)
       gAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcggc                                                                                                                                            >  2:815460/1‑149 (MQ=255)
         aGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCat                                                                                                                                          >  4:261868/1‑149 (MQ=255)
         aGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCat                                                                                                                                          <  2:1055629/149‑1 (MQ=255)
           aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattc                                                                                                                                        <  4:715889/149‑3 (MQ=255)
           aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattc                                                                                                                                        <  1:391706/149‑3 (MQ=255)
           aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattc                                                                                                                                        <  1:472199/149‑3 (MQ=255)
           aTACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattc                                                                                                                                        <  3:1001031/149‑3 (MQ=255)
             aCGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcgg                                                                                                                                      <  1:946554/149‑5 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                    >  2:1020154/1‑149 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                    <  1:963118/149‑1 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                    <  2:401800/149‑1 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                    >  1:288175/1‑149 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTAACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                    <  4:74498/149‑1 (MQ=255)
                gTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAg                                                                                                                                   <  3:1013173/149‑1 (MQ=255)
                gTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATGTTAg                                                                                                                                   >  4:936904/1‑149 (MQ=255)
                  gCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGCTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGAt                                                                                                                                 <  1:815460/149‑1 (MQ=255)
                   cTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGAtt                                                                                                                                >  1:301390/1‑149 (MQ=255)
                    ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                               >  2:668791/1‑149 (MQ=255)
                    ttttGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                               >  3:576205/1‑149 (MQ=255)
                           tCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcgg                                                                                                                                      <  2:252788/135‑5 (MQ=255)
                           tCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcgg                                                                                                                                      >  1:252788/1‑131 (MQ=255)
                             gTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAttt                                                                                                                                      >  2:657731/1‑133 (MQ=255)
                             gTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAttt                                                                                                                                      <  1:657731/133‑1 (MQ=255)
                                cTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATc                                                                                                                             <  2:1051412/139‑1 (MQ=255)
                                cTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATc                                                                                                                             >  1:1051412/1‑139 (MQ=255)
                                cTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGtt                                                                                                                   >  1:976628/1‑149 (MQ=255)
                                cTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGtt                                                                                                                   <  2:977043/149‑1 (MQ=255)
                                   aTTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaa                                                                                                                >  2:551716/1‑149 (MQ=255)
                                    ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                                                               <  3:824817/149‑1 (MQ=255)
                                    ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                                                               <  3:261868/149‑1 (MQ=255)
                                     tGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGAtt                                                                                                                        >  3:509333/1‑139 (MQ=255)
                                       tttCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTc                                                                                                            <  3:1023238/149‑1 (MQ=255)
                                       tttCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTc                                                                                                            >  3:642576/1‑149 (MQ=255)
                                         tCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAGCCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCAc                                                                                                          <  4:754927/149‑1 (MQ=255)
                                              tGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaa                                                                                                     >  3:561273/1‑149 (MQ=255)
                                               gcgcTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaa                                                                                                    <  2:441139/149‑1 (MQ=255)
                                                 gcTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaaaa                                                                                                  <  3:944745/149‑1 (MQ=255)
                                                  cTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                                                                 <  3:301439/149‑1 (MQ=255)
                                                   tGTTAACCGAGGTAAATAATAACCGGAGGCTGTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGACTGCCGCCCATTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAacc                                                                                                >  4:644595/1‑149 (MQ=255)
                                                       aaCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCAccc                                                                                            >  1:679244/1‑149 (MQ=255)
                                                               tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGAc                                                                                     <  3:962785/148‑1 (MQ=255)
                                                                 aataatAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGtt                                                                                  >  1:900629/1‑149 (MQ=255)
                                                                  ataataACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTc                                                                                 >  3:488393/1‑149 (MQ=255)
                                                                   taataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCa                                                                                >  1:767543/1‑149 (MQ=255)
                                                                   taataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCa                                                                                >  1:872857/1‑149 (MQ=255)
                                                                        aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                                                                           <  4:642576/149‑1 (MQ=255)
                                                                        aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGaa                                                                           <  2:577912/149‑1 (MQ=255)
                                                                           ggAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATgg                                                                        <  3:936904/149‑1 (MQ=255)
                                                                             aGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGgcg                                                                      <  4:1038008/149‑1 (MQ=255)
                                                                               tctctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCa                                                                    >  4:567368/1‑149 (MQ=255)
                                                                               tctctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCa                                                                    >  4:523834/1‑149 (MQ=255)
                                                                                 tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                                                                  <  4:1041958/149‑1 (MQ=255)
                                                                                 tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                                                                  >  3:292967/1‑149 (MQ=255)
                                                                                 tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                                                                  <  3:285591/149‑1 (MQ=255)
                                                                                       gCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                            <  2:722457/149‑1 (MQ=255)
                                                                                          aCAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGt                                                         >  2:272935/1‑149 (MQ=255)
                                                                                            aaTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaa                                                       >  4:297561/1‑149 (MQ=255)
                                                                                             aTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaat                                                      >  4:79966/1‑149 (MQ=255)
                                                                                                 aCTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCCGCATATTTTGGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAg                                                  <  1:935687/149‑1 (MQ=255)
                                                                                                   tggtggTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                            >  3:1009796/1‑137 (MQ=255)
                                                                                                   tggtggTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTg                                                >  1:75847/1‑149 (MQ=255)
                                                                                                     gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGc                                              >  3:150651/1‑149 (MQ=255)
                                                                                                     gtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGc                                              <  3:309302/149‑1 (MQ=255)
                                                                                                      tggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGTCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCt                                             >  2:936482/1‑149 (MQ=255)
                                                                                                         ttAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGTCTAAt                                          >  1:47998/1‑149 (MQ=255)
                                                                                                           aacaacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTc                                        >  2:757237/1‑149 (MQ=255)
                                                                                                              aacCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAaccac                                                                                              >  1:660522/1‑92 (MQ=255)
                                                                                                               acCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCttt                                    >  1:671875/1‑149 (MQ=255)
                                                                                                               acCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCttt                                    <  2:671875/149‑1 (MQ=255)
                                                                                                                 cTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCttttt                                  >  4:288325/1‑149 (MQ=255)
                                                                                                                 cTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCttttt                                  >  1:493867/1‑149 (MQ=255)
                                                                                                                         cagcaAGTAAGCCAGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                          <  3:288325/149‑1 (MQ=255)
                                                                                                                          agcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCttttt                         >  4:333454/1‑149 (MQ=255)
                                                                                                                          agcaAGTAAGCCCGAATGCCGCCCTTTGCGCCGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATCGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCttttt                         >  4:644246/1‑149 (MQ=255)
                                                                                                                           gcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttttt                        >  3:742067/1‑149 (MQ=255)
                                                                                                                              aGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCAc                     <  1:757237/149‑1 (MQ=255)
                                                                                                                                tAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCt                   <  3:661794/149‑1 (MQ=255)
                                                                                                                                  aGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGa                 <  1:919572/149‑1 (MQ=255)
                                                                                                                                  aGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGa                 >  2:226111/1‑149 (MQ=255)
                                                                                                                                   taaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTa                                            >  1:1001671/7‑121 (MQ=255)
                                                                                                                                   taaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTa                                            <  2:1001671/115‑1 (MQ=255)
                                                                                                                                   taaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGAc                >  4:185005/7‑149 (MQ=255)
                                                                                                                                   taaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGAc                >  4:666876/7‑149 (MQ=255)
                                                                                                                                   taaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAAACACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTGTCCCCTTTTTTCACCTGAc                >  3:50369/7‑149 (MQ=255)
                                                                                                                                   gCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGAc                <  2:199948/149‑1 (MQ=255)
                                                                                                                                     ccGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACaa              >  4:22992/1‑149 (MQ=255)
                                                                                                                                     ccGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACaa              >  1:590778/1‑149 (MQ=255)
                                                                                                                                      cGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACaaa             <  3:567368/149‑1 (MQ=255)
                                                                                                                                       gAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATaaa                                                   <  4:52059/110‑1 (MQ=255)
                                                                                                                                          tGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTTGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATc         <  2:911536/149‑1 (MQ=255)
                                                                                                                                             cgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGAt      >  3:837205/1‑149 (MQ=255)
                                                                                                                                               cccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTg    <  1:272935/149‑1 (MQ=255)
                                                                                                                                               cccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTg    <  1:23497/149‑1 (MQ=255)
                                                                                                                                               cccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTg    >  4:784800/1‑149 (MQ=255)
                                                                                                                                                 cTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGt  <  4:300956/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGT  >  NC_000913/2653996‑2654289

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: