Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,654,142 T→A 5.6% noncoding (288/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,1420TA 5.6% 321.5 / 6.9 117noncoding (288/304 nt)ryfAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (61/48);  new base A (3/4);  total (64/53)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 117 total aligned reads displayed.

GCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTT  >  NC_000913/2653996‑2654290
                                                                                                                                                  |                                                                                                                                                     
gCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTg                                                                                                                                                    >  2:1555374/1‑150 (MQ=255)
    aaTGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                                                                <  1:907414/150‑1 (MQ=255)
    aaTGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                                                                <  2:1067851/150‑1 (MQ=255)
    aaTGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGggcg                                                                                                                                                <  2:607103/150‑1 (MQ=255)
      tGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGggcggc                                                                                                                                              <  1:118645/150‑1 (MQ=255)
       gAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCa                                                                                                                                             <  2:622009/150‑1 (MQ=255)
          gATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCatat                                                                                                                                         >  1:1630281/1‑150 (MQ=255)
            tACGGGGCTTATGTCTCGTACTTATTGTTGCTGGTGCGCTGTTCACCGCGGTAAATAATCACCGGAGTCTCTCCGGCGACAATTTACTCGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcgg                                                                                                                                        <  1:1667752/150‑5 (MQ=255)
             aCGGTGCTTTTGTTTCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCAAAGGGCGGCattcggg                                                                                                                                       <  2:1250032/150‑6 (MQ=255)
             aCGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATAtttt                                                                                                                                      >  2:951091/1‑150 (MQ=255)
              cGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                     >  1:246995/1‑150 (MQ=255)
              cGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTa                                                                                                                                     >  1:1761265/1‑150 (MQ=255)
               ggTGCTTTTGTTTCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAg                                                                                                                                    <  2:1655011/150‑1 (MQ=255)
               ggTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGTTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAg                                                                                                                                    <  2:60102/150‑1 (MQ=255)
                 tGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGAt                                                                                                                                  <  1:869147/150‑1 (MQ=255)
                 tGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGAt                                                                                                                                  <  1:370724/150‑1 (MQ=255)
                  gccgttGTATCGTACTTCTTGTTGCTGGGGCGCTGTTAACCGAGGTACAAAATAACCGGAGTCTCTCCGGCGACAATATACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAAGGCCGCCCCTTGGGCGGCATATTTTAGAtt                                                                                                                                 <  1:284439/146‑1 (MQ=255)
                   cTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTa                                                                                                                                >  2:948457/1‑150 (MQ=255)
                      ttGTCTCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATcc                                                                                                                             <  1:1361523/150‑1 (MQ=255)
                       tGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCg                                                                                                                            >  2:186/1‑150 (MQ=255)
                                   aTTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTaaa                                                                                                                <  2:508715/150‑1 (MQ=255)
                                    ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAg                                                                                                               <  2:1795662/150‑1 (MQ=255)
                                    ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAg                                                                                                               <  1:1656091/150‑1 (MQ=255)
                                    ttGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGGGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAg                                                                                                               <  2:336770/150‑1 (MQ=255)
                                        ttCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGCGCGGCATATTTTACATTATCCGATTCTGTTTAAAGTCAc                                                                                                           >  2:1622679/1‑150 (MQ=255)
                                          cTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGc                                                                                                         <  2:35462/150‑1 (MQ=255)
                                           tGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACACCCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCa                                                                                                        >  2:230713/1‑150 (MQ=255)
                                            ggTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaa                                                                                                       >  1:1159798/1‑150 (MQ=255)
                                              tGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCaaaa                                                                                                     >  1:1818859/1‑150 (MQ=255)
                                               gcgcTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGTTGGTTAACAACCTTCAGAGCCGCAAGTAAGCCCGAATGCCGCCCTTTTTGCGGCATATTTTAGATTATAAGATTCTGTTTAAAGTCACGCaaaaa                                                                                                    >  1:1645683/1‑150 (MQ=255)
                                                 gcTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTGTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAac                                                                                                  >  2:759737/1‑150 (MQ=255)
                                                     ttAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAaccacc                                                                                              <  1:1485793/150‑1 (MQ=255)
                                                      tAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCAccc                                                                                             >  1:60627/1‑150 (MQ=255)
                                                       aaCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCAcccc                                                                                            >  2:384535/1‑150 (MQ=255)
                                                       aaCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCAcccc                                                                                            >  2:1101168/1‑150 (MQ=255)
                                                       aaCCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCAcacc                                                                                            >  2:416326/1‑150 (MQ=255)
                                                         ccGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAg                                                                                          >  1:1055007/1‑150 (MQ=255)
                                                             ggTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCTGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGAc                                                                                      >  1:947641/1‑150 (MQ=255)
                                                             ggTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGAAAACAACCTTCAGAGCAGCAAGTAAGCCCGAAAGCCGCCCTTTGGGCGGCATATTTTAGATTAACCCATTCAGTTTAAAGTCACGCACACAACCAcccccgagcc                                                                                      >  1:713466/1‑144 (MQ=255)
                                                               tAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGt                                                                                    >  2:983104/1‑150 (MQ=255)
                                                                aaaTAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGtt                                                                                   <  2:1822640/150‑1 (MQ=255)
                                                                aaaTAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCACGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGAATCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGtt                                                                                   >  2:586702/1‑150 (MQ=255)
                                                                  ataataACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCa                                                                                 <  1:587536/150‑1 (MQ=255)
                                                                  ataataACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGAATATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCa                                                                                 <  1:1611402/150‑1 (MQ=255)
                                                                  ataataACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCAGCCCTTTGGGCGGCAAATTTTACCTTAACCTAATCTGTTTAAATTCACGCAAAATAAGACCCCAGAGACTTTCa                                                                                 >  1:1151694/1‑150 (MQ=255)
                                                                    aataaCCGGAGTCTCTCCTGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATa                                                                               <  1:759737/150‑1 (MQ=255)
                                                                    aataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATa                                                                               <  2:243228/150‑1 (MQ=255)
                                                                     ataaCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACGTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAg                                                                              <  2:961513/150‑1 (MQ=255)
                                                                        aCCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAAt                                                                           >  2:1348244/1‑150 (MQ=255)
                                                                          cGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTATAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATgg                                                                         >  2:693304/1‑150 (MQ=255)
                                                                              gTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCa                                                                     >  1:478073/1‑150 (MQ=255)
                                                                                ctctcCTGCGAAAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCAGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCTAAAAACCACCCCAGCGACGTTCATAGAATTGCGCAat                                                                   <  1:1282714/150‑1 (MQ=255)
                                                                                ctctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAat                                                                   >  1:427500/1‑150 (MQ=255)
                                                                                 tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTTTACGATTCTGTTTAAAGTAAAGCAAAAAACCACCCCAGAGACGTTAATAGAATGGCGCAata                                                                  >  1:1312882/1‑150 (MQ=255)
                                                                                 tctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAata                                                                  >  1:1275333/1‑150 (MQ=255)
                                                                                  ctcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAatat                                                                 >  1:94161/1‑150 (MQ=255)
                                                                                   tcCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATg                                                                <  1:402707/150‑1 (MQ=255)
                                                                                    ccGGCGACAATTTACTGGTGGTTAACATCCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTGTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGGCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATTTGt                                                               >  1:1119650/1‑150 (MQ=255)
                                                                                      ggCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTc                                                             <  2:1761878/150‑1 (MQ=255)
                                                                                       gCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCg                                                            <  1:948457/150‑1 (MQ=255)
                                                                                             aTTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGtaata                                                      <  2:246995/150‑1 (MQ=255)
                                                                                             aTTTACTGGTGGTTAACAACCTTAAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTATATTATCAGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGAGCAATATGTTCGGtaata                                                      >  1:1500289/1‑150 (MQ=255)
                                                                                                  cTGGTGGTTGACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAAACACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTg                                                 >  2:1201212/1‑150 (MQ=255)
                                                                                                    ggtggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGc                                               <  1:1555374/150‑1 (MQ=255)
                                                                                                      gggtTAACAACCTTCAGAGCATGAAGTAAGCCCGAAAGCCGCCCTTTTGGCGGCATATTTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCt                                              <  2:504290/149‑1 (MQ=255)
                                                                                                       ggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTaa                                            >  2:1480272/1‑150 (MQ=255)
                                                                                                       ggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTaa                                            >  2:1286472/1‑150 (MQ=255)
                                                                                                       ggtTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTaa                                            >  2:1280041/1‑150 (MQ=255)
                                                                                                                ccTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCttttt                                   >  1:812459/1‑150 (MQ=255)
                                                                                                                 cTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTc                                  >  1:1105004/1‑150 (MQ=255)
                                                                                                                 cTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCATTGGGCGGCATATTTTAGATTATACGATTCTGTTTAAATTCACGCAAAAAACCCCCCCCGCCACGTTCAAATAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCAttttta                                  >  1:1499060/1‑149 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAATAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTcc                                 <  1:762766/150‑1 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACTTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTcc                                 <  2:1609275/150‑1 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTcc                                 <  2:1630281/150‑1 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACAACCCCAGCGACGTTCATATAATGGCGCAATATTTTCGGTAATAAAGTGGCTAATTCCTTTTTcc                                 >  1:1644842/1‑150 (MQ=255)
                                                                                                                  ttCAGAGCAGCAAGTAAGCCAGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATTGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTcc                                 <  2:733356/150‑1 (MQ=255)
                                                                                                                   tCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTccc                                >  1:1542345/1‑150 (MQ=255)
                                                                                                                    cAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTcccc                               <  1:62190/150‑1 (MQ=255)
                                                                                                                        gcagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                           <  2:980638/150‑1 (MQ=255)
                                                                                                                        gcagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                           >  2:1672137/1‑150 (MQ=255)
                                                                                                                        gcagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttt                           >  1:1210644/1‑150 (MQ=255)
                                                                                                                         cagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCttttt                          <  2:594444/150‑1 (MQ=255)
                                                                                                                         cagcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCttttt                          <  2:1215388/150‑1 (MQ=255)
                                                                                                                          agcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCtttttt                         >  2:636496/1‑150 (MQ=255)
                                                                                                                           gcaAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTc                        >  1:1819215/1‑150 (MQ=255)
                                                                                                                           gcaAGTAGCCCGAATGCCCCCCTTTGGGCGGCATATTTTAGATTAAACGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTAATAGAATGGCGCAATATGTTAGGTAATAAAGTGGCTAATTCCTTTTTACCCatttttaa                       >  2:1753299/1‑148 (MQ=255)
                                                                                                                            caAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGCATGGCGCAATATGTTCGGTAATAAAGTGGCTAGTTGCTTTTTCCCCTTTTTTCa                       >  2:1103588/1‑150 (MQ=255)
                                                                                                                            caAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCa                       <  2:1055007/150‑1 (MQ=255)
                                                                                                                             aaGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCAc                      >  2:672110/1‑150 (MQ=255)
                                                                                                                                tAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTg                   >  1:1692743/1‑150 (MQ=255)
                                                                                                                                   taaaataTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACa                >  2:1399222/7‑150 (MQ=255)
                                                                                                                                    cccGAATTCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGATTGGCGCAATATGTTCGGTAATATAGTGGCTAATTCCTTTTTCCCCTTTTTTCACATGACaa               >  2:708077/1‑150 (MQ=255)
                                                                                                                                         aTGCCGCCCAAAGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTATTAAAGTGGCTAATTCCTTTTTCCCCGTTTTTCACCTGACAAATATc          >  1:173002/1‑150 (MQ=255)
                                                                                                                                           gccgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGGTCATAGAATGGCGCAATATGTTAGGTAATAAAGTGGCTAATTAATTTTTCCCATTTTTTGACCTGACAAATATAGa        >  2:751157/1‑150 (MQ=255)
                                                                                                                                           gccgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTAACGCAAAAAAACACCCCCGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTATTTTCACCTGACAAATATCGa        >  1:1401136/1‑150 (MQ=255)
                                                                                                                                            ccgccCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCCCGCAAAAAACAACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTACCCTTTTTTCACATGACAAATATCGAt       >  1:1472616/1‑150 (MQ=255)
                                                                                                                                               cccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTgg    <  2:1146171/150‑1 (MQ=255)
                                                                                                                                                ccTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGt   >  2:1446908/1‑150 (MQ=255)
                                                                                                                                                 cTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGtt  <  1:384535/150‑1 (MQ=255)
                                                                                                                                                 cTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTACCCTTTTTTCACCTGACAAATATCGATTGGtt  >  2:1243668/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCACAATGAAGATACGGTGCTTTTGTATCGTACTTATTGTTTCTGGTGCGCTGTTAACCGAGGTAAATAATAACCGGAGTCTCTCCGGCGACAATTTACTGGTGGTTAACAACCTTCAGAGCAGCAAGTAAGCCCGAATGCCGCCCTTTGGGCGGCATATTTTAGATTATCCGATTCTGTTTAAAGTCACGCAAAAAACCACCCCAGCGACGTTCATAGAATGGCGCAATATGTTCGGTAATAAAGTGGCTAATTCCTTTTTCCCCTTTTTTCACCTGACAAATATCGATTGGTT  >  NC_000913/2653996‑2654290

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: