Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,944 G→T 8.3% noncoding (123/140 nt) ryjA ← novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,277,9440GT 8.3% 400.6 / 17.5 147noncoding (123/140 nt)ryjAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base G (67/67);  new base T (3/10);  total (70/77)
Fisher's exact test for biased strand distribution p-value = 8.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 150 total aligned reads displayed.

AGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTTC  >  NC_000913/4277795‑4278093
                                                                                                                                                     |                                                                                                                                                     
   ggCGCTGCGTGACGAACTGGACGGATGTATTGGGTGTGGCTGCCTTTCGCGCAGTGATTTCCCGTTGCGTAACCCTGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgc                                                                                                                                                    <  2:1063295/150‑1 (MQ=255)
       cTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCtt                                                                                                                                                <  2:510083/150‑1 (MQ=255)
        tGCGTGACGAACTGGACGGATGTATTGGGTGTGGCTGCCTTTCGCGCAGTGATTGCCCGGTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCttt                                                                                                                                               <  2:1338667/150‑1 (MQ=255)
          ggTGACTAACTGGACGGATGTATTGGGTGTGTCTGCCTTTCGCGCAGTGATTGCCCGTGGCGTCACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAg                                                                                                                                             <  1:611018/149‑1 (MQ=255)
          cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAg                                                                                                                                             <  1:407023/150‑1 (MQ=255)
          cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGCGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAg                                                                                                                                             <  2:650965/150‑1 (MQ=255)
           gTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAgt                                                                                                                                            <  2:417649/150‑1 (MQ=255)
            tGACTAACTGTACGGATGTATTGGTGGTGGCTGCCTTACGAGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                           <  2:879032/150‑1 (MQ=255)
            tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                           <  2:328268/150‑1 (MQ=255)
            tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                           >  2:1280892/1‑150 (MQ=255)
            tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                           >  2:1165038/1‑150 (MQ=255)
              ccGAACTGGAAGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGACCGTTGCGTAACCCGGGCGACCGCTTAGGAGAATAAGGTACCGGCGCACGCTTGCTGGAAGATGAACTAAACTAAAGCGCCCTTGTGGCGCTTTAgtttt                                                                                                                                         <  2:308816/149‑2 (MQ=255)
              aCGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAgtttt                                                                                                                                         <  1:1183975/150‑2 (MQ=255)
                  aCTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGtttt                                                                                                                                     >  1:502104/1‑150 (MQ=255)
                   caggACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTATTc                                                                                                                                    >  2:1683347/3‑150 (MQ=255)
                     ggtcGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCg                                                                                                                                  <  1:522647/147‑1 (MQ=255)
                     ggACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCg                                                                                                                                  >  2:944286/1‑150 (MQ=255)
                       ccGGAAGTATTGGTTGTGGCTCCCTTTCGCGCAGGGATTGCCCGGTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTAAGTTTGTTTTCCGgt                                                                                                                                <  1:1528979/149‑1 (MQ=255)
                            tGTATTGGTGGTGGCTGCCTTTCGCGCAGGGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttg                                                                                                                           <  1:134626/150‑1 (MQ=255)
                            ggTATTGGTTGTGGCTGCCTATGGCGCAGTGATTGCCCGTTGCGTAAACCGGGCGTCCGCTTAGTAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttg                                                                                                                           <  1:129359/149‑1 (MQ=255)
                             gTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCGTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgt                                                                                                                          >  2:778025/1‑150 (MQ=255)
                              tATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtc                                                                                                                         <  2:1022934/150‑1 (MQ=255)
                              tATTGGTTGTGGCTGCCTTTCGCGCAGTGATAGCCCTTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAATGTACAGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCAGgtctttgtc                                                                                                                         >  1:238103/1‑150 (MQ=255)
                              tATTGGATGTGGCTGCCTTTCGAGCAGTGATTGCCCGTTGCGTAACCCGGGAGACCGATTAGGAGAAGAAGGTACCGCCGCACGATTGCTGGAAGATCAACAAAAAAAAAGCCCCACATGGCCAGTTTATTTTGTTTTCGGgaatttgac                                                                                                                         >  1:580994/1‑148 (MQ=255)
                               aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtct                                                                                                                        >  1:1621211/1‑150 (MQ=255)
                                 tGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAAAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtcttt                                                                                                                      >  2:1236750/1‑150 (MQ=255)
                                 tGGTTGTGGCTGCCTTTCGCGCAGGGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtcttt                                                                                                                      <  1:1165021/150‑1 (MQ=255)
                                  ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtc                                                                                                                     >  2:1434971/1‑150 (MQ=255)
                                    ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTGGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctc                                                                                                                   <  2:1496270/150‑1 (MQ=255)
                                      gtgGCCGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGTAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctctc                                                                                                                 <  1:336642/150‑1 (MQ=255)
                                        ggCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTa                                                                                                               >  2:227766/1‑150 (MQ=255)
                                         gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTAt                                                                                                              >  2:1092511/1‑150 (MQ=255)
                                            gCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGCCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCATTAGTTAGTTTTCCGGTCTTTGTCTTTCTCTCTATccc                                                                                                           >  2:1666186/1‑150 (MQ=255)
                                             ccTTTCGCGCAGTGATTGCCCGTGGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCg                                                                                                          <  1:1224226/150‑1 (MQ=255)
                                              cTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTGGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGc                                                                                                         <  2:52704/150‑1 (MQ=255)
                                                    cgcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTac                                                                                                   >  1:144220/1‑150 (MQ=255)
                                                    cgcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTac                                                                                                   >  1:139057/1‑150 (MQ=255)
                                                    cgcAGTGATTGCCCGTTGCGGAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTac                                                                                                   <  1:1366150/150‑1 (MQ=255)
                                                       aGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTAcaca                                                                                                <  1:186446/150‑1 (MQ=255)
                                                          gATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACagga                                                                                             >  1:788327/1‑150 (MQ=255)
                                                          gATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTGTCCGGTCTTTGTCTGTCTCTCTATCCCGCTGGTACACagga                                                                                             >  1:841948/1‑150 (MQ=255)
                                                             tGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAggg                                                                                          <  2:1267271/150‑1 (MQ=255)
                                                                 cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTc                                                                                      <  2:504957/150‑1 (MQ=255)
                                                                 cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTc                                                                                      >  2:127633/1‑150 (MQ=255)
                                                                   ttGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccc                                                                                    <  1:775639/150‑1 (MQ=255)
                                                                    tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTcccc                                                                                   <  2:1485812/150‑1 (MQ=255)
                                                                    tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTcccc                                                                                   >  1:65552/1‑150 (MQ=255)
                                                                     gCGTAACCCGGGAGACCGCTTAGGAGAAGAAGGTACCTGCGCACGCTTGATGGAAGATGAACAAAACTAAAGCGCCACAATGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccccc                                                                                  <  2:350871/150‑1 (MQ=255)
                                                                      cGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCg                                                                                 <  2:709459/150‑1 (MQ=255)
                                                                       gTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGa                                                                                >  2:954048/1‑150 (MQ=255)
                                                                          aCCCGGTCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTAGCTGGAAGAAGAACAAAACTAAAGCGCAACAAGGGCGCGTTAGTTTGTTTTCCGGTCTTTCTCTTTCTCTCTATCCCGCTTGTACACAGGACGGTTTCCCCCGACGt                                                                             >  1:182344/1‑150 (MQ=255)
                                                                             cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                          >  2:1368338/1‑150 (MQ=255)
                                                                             cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                          >  1:1493650/1‑150 (MQ=255)
                                                                               ggCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAaca                                                                        <  1:954048/150‑1 (MQ=255)
                                                                                   aCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACAcc                                                                    >  1:737484/1‑150 (MQ=255)
                                                                                      gCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCa                                                                 >  1:1668693/1‑150 (MQ=255)
                                                                                       cTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCAt                                                                >  1:1513005/1‑150 (MQ=255)
                                                                                           ggAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAAAACACCTCATTCGa                                                            >  1:374920/1‑150 (MQ=255)
                                                                                            gagaAGAAGGTACCGGCGCATGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACATCATTagag                                                           >  1:691801/1‑150 (MQ=255)
                                                                                            gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                           <  1:208664/150‑1 (MQ=255)
                                                                                            gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                           >  1:221590/1‑150 (MQ=255)
                                                                                            gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                           <  2:139057/150‑1 (MQ=255)
                                                                                              gaagaaGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGCCGTCAACACACCTCATTCGAGCa                                                         >  2:35488/1‑150 (MQ=255)
                                                                                              gaagaaGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCa                                                         >  1:1378194/1‑150 (MQ=255)
                                                                                              gaagaaGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCa                                                         >  1:348179/1‑150 (MQ=255)
                                                                                                   agGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgtgg                                                    <  1:1165038/150‑1 (MQ=255)
                                                                                                        ccGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                               <  1:944286/150‑1 (MQ=255)
                                                                                                        ccGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                               >  1:657732/1‑150 (MQ=255)
                                                                                                        ccGGCGCACGCGTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                               <  2:114070/150‑1 (MQ=255)
                                                                                                         cGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAgg                                              >  1:880973/1‑150 (MQ=255)
                                                                                                          ggCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGt                                             >  1:739702/1‑150 (MQ=255)
                                                                                                             gcACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGATGGTTTCCCCCGCCGTCAACACCCCTCATTCGAGCACGGGGTGGAggatac                                          >  2:266344/1‑146 (MQ=255)
                                                                                                                  cTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTg                                     <  2:673951/150‑1 (MQ=255)
                                                                                                                      cgggAAGATGCACAAAACTCAAGCGCCACCAGGGCGCTTTTGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGGTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                 <  1:1368785/148‑1 (MQ=255)
                                                                                                                      cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCACCACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                 >  1:1452674/1‑150 (MQ=255)
                                                                                                                      cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                 <  1:1100048/150‑1 (MQ=255)
                                                                                                                      cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGATGgtgt                                 >  2:191109/1‑150 (MQ=255)
                                                                                                                      cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCGTTAGTTTGTTTTCCGGGCTTTGGCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                 <  2:1401834/150‑1 (MQ=255)
                                                                                                                        ggAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                               <  2:221590/150‑1 (MQ=255)
                                                                                                                        ggAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                               >  2:1423395/1‑150 (MQ=255)
                                                                                                                        ggAAGAAGAACAGAACTACATCGCCACAAGGGCGCTTTAGTTTGTTGTCCGGTCTTTGTCTTTCTCTCTATCCCGCTTGTACAGCGTTGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                               <  2:1514812/150‑1 (MQ=255)
                                                                                                                          aatatcaaCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACCCAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGAt                             <  1:887341/144‑1 (MQ=255)
                                                                                                                          aaGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTAGAt                             >  2:462728/1‑150 (MQ=255)
                                                                                                                              tGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCtt                         >  1:236018/1‑150 (MQ=255)
                                                                                                                               taaCAAAACTAAATCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGTGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCttt                        <  2:886300/149‑1 (MQ=255)
                                                                                                                               gAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGATGTTCCGGTTGGTGTTGATGCttt                        <  2:1585337/150‑1 (MQ=255)
                                                                                                                                aaCAAAACTAAAGATCAACAAGGGCGATTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTGTATCCCTCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCTAGAACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTa                       <  2:182344/150‑1 (MQ=255)
                                                                                                                                   aaaaCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAtt                    >  2:1013347/1‑150 (MQ=255)
                                                                                                                                     aaCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGt                  >  1:934283/1‑150 (MQ=255)
                                                                                                                                      aCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa                 >  1:553450/1‑150 (MQ=255)
                                                                                                                                      aCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa                 <  1:854196/150‑1 (MQ=255)
                                                                                                                                      aCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGAGGTTACGGGTGGTGGTGATTCTTTAATTGTAt                >  2:1289061/1‑150 (MQ=255)
                                                                                                                                        ttaaccgcgACTAGGGTGGTTTAGTTTGTTTTCCGGTATTTGTCTTTCTCTGTATCGGTTTTGTAAACAGTAGGGTTTCCCCCGACGTCTACACACCTCATTCGAGAACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTATTTGTATg               <  2:580994/145‑1 (MQ=255)
                                                                                                                                        aaaaGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTAGACAGTATGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTTGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATg               <  1:958796/149‑1 (MQ=255)
                                                                                                                                           aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa            <  1:1423395/150‑1 (MQ=255)
                                                                                                                                           aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa            >  1:1663267/1‑150 (MQ=255)
                                                                                                                                           aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa            <  1:1248553/150‑1 (MQ=255)
                                                                                                                                               ccACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCgc        >  1:1133482/1‑149 (MQ=255)
                                                                                                                                                    aGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGttt   >  2:1100992/1‑150 (MQ=255)
                                                                                                                                                     gggCGCTTTATTTTGTTTTCCGGTCTTTCTCTTACACTCTATCCCGCTGGTACACCCCAGGGTTTCCCCCGCCGTCAACACACCTCATGCGCGCACGGGGTGGAGGGTCCTGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTttc  <  2:293013/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTTC  >  NC_000913/4277795‑4278093

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: