Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,943 A→G 7.6% noncoding (124/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,9430AG 7.6% 420.7 / 19.0 146noncoding (124/140 nt)ryjAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base A (66/68);  new base G (3/8);  total (70/76)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

Only 100 of 148 total aligned reads displayed.

AGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTT  >  NC_000913/4277795‑4278092
                                                                                                                                                    |                                                                                                                                                     
   ggCGCTGCGTGACGAACTGGACGGATGTATTGGGTGTGGCTGCCTTTCGCGCAGTGATTTCCCGTTGCGTAACCCTGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgc                                                                                                                                                   <  2:1063295/150‑1 (MQ=255)
       cTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCtt                                                                                                                                               <  2:510083/150‑1 (MQ=255)
        tGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGAGAATGCCCGTGGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCttt                                                                                                                                              <  1:1528473/150‑1 (MQ=255)
          ggTGACTAACTGGACGGATGTATTGGGTGTGTCTGCCTTTCGCGCAGTGATTGCCCGTGGCGTCACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAg                                                                                                                                            <  1:611018/149‑1 (MQ=255)
          cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGCGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAg                                                                                                                                            <  2:650965/150‑1 (MQ=255)
           gTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAgt                                                                                                                                           <  2:417649/150‑1 (MQ=255)
           gTGAAGAACTGGACTGATGTATTGGGTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgt                                                                                                                                           <  1:1489975/150‑1 (MQ=255)
            tGACTAACTGTACGGATGTATTGGTGGTGGCTGCCTTACGAGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                          <  2:879032/150‑1 (MQ=255)
            tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCTCCACAAGGGCGCTTTAgtt                                                                                                                                          >  1:254476/1‑150 (MQ=255)
            tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                          <  2:328268/150‑1 (MQ=255)
            tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                          >  2:1165038/1‑150 (MQ=255)
               cGAACTGGACGGATGTATTGGTGGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtttgt                                                                                                                                       <  1:1380490/150‑1 (MQ=255)
                gAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAgttttgt                                                                                                                                      <  1:916444/150‑4 (MQ=255)
                 aaCTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAgttttgtt                                                                                                                                     <  1:485010/150‑5 (MQ=255)
                  aCTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGtttt                                                                                                                                    >  1:502104/1‑150 (MQ=255)
                   caggACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTATTc                                                                                                                                   >  2:1683347/3‑150 (MQ=255)
                     ggACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCg                                                                                                                                 >  2:944286/1‑150 (MQ=255)
                            tGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATGGCCCGTGGCGTAACCCGGGCGACCGCTTAGGAGAAGACGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttg                                                                                                                          <  1:535790/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)
                               aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtct                                                                                                                       >  1:1621211/1‑150 (MQ=255)
                                 tGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAAAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtcttt                                                                                                                     >  2:1236750/1‑150 (MQ=255)
                                  ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGCCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCTCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtc                                                                                                                    >  2:1368785/1‑150 (MQ=255)
                                  ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtc                                                                                                                    >  2:1434971/1‑150 (MQ=255)
                                    ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTGGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctc                                                                                                                  <  2:1496270/150‑1 (MQ=255)
                                      gtgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctctc                                                                                                                <  1:1319690/150‑1 (MQ=255)
                                      gtgGCCGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGTAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctctc                                                                                                                <  1:336642/150‑1 (MQ=255)
                                         gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTAt                                                                                                             >  2:1092511/1‑150 (MQ=255)
                                             ccTTTCGCGCAGTGATTGCCCGTGGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCg                                                                                                         <  1:1224226/150‑1 (MQ=255)
                                              cTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTGGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGc                                                                                                        <  2:52704/150‑1 (MQ=255)
                                              cTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGc                                                                                                        >  2:208664/1‑150 (MQ=255)
                                               tttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCt                                                                                                       >  2:324083/1‑150 (MQ=255)
                                                    cgcAGTGATTGCCCGTTGCGGAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTac                                                                                                  <  1:1366150/150‑1 (MQ=255)
                                                       aGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTAcaca                                                                                               <  1:186446/150‑1 (MQ=255)
                                                         tGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACagg                                                                                             <  1:1637731/150‑1 (MQ=255)
                                                          gATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTGTCCGGTCTTTGTCTGTCTCTCTATCCCGCTGGTACACagga                                                                                            >  1:841948/1‑150 (MQ=255)
                                                           tttGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggag                                                                                           <  1:814997/149‑1 (MQ=255)
                                                           aTTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggag                                                                                           <  1:59848/150‑1 (MQ=255)
                                                           aTTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCGTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggag                                                                                           >  1:1730503/1‑150 (MQ=255)
                                                                 cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTc                                                                                     <  2:504957/150‑1 (MQ=255)
                                                                 cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTc                                                                                     >  2:127633/1‑150 (MQ=255)
                                                                   ttGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccc                                                                                   <  1:775639/150‑1 (MQ=255)
                                                                     gCGTAACCCGGGAGACCGCTTAGGAGAAGAAGGTACCTGCGCACGCTTGATGGAAGATGAACAAAACTAAAGCGCCACAATGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccccc                                                                                 <  2:350871/150‑1 (MQ=255)
                                                                      cGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCg                                                                                <  2:709459/150‑1 (MQ=255)
                                                                       gTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGa                                                                               >  2:954048/1‑150 (MQ=255)
                                                                           cccGGGCGCCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCAGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTGAGTGTGTGTTCCGGTCTTTGGCTGTCGCTCTATCCCGCTGGGACACAGGAGGGGTTCCCCCGACGta                                                                           >  1:204655/1‑149 (MQ=255)
                                                                             cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                         >  2:1368338/1‑150 (MQ=255)
                                                                               ggCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAaca                                                                       <  1:954048/150‑1 (MQ=255)
                                                                                   aCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACAcc                                                                   >  1:657602/1‑150 (MQ=255)
                                                                                   aCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACAcc                                                                   >  1:737484/1‑150 (MQ=255)
                                                                                      gCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCa                                                                >  1:1668693/1‑150 (MQ=255)
                                                                                       cTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCAt                                                               >  1:1513005/1‑150 (MQ=255)
                                                                                          aGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCg                                                            <  2:127832/150‑1 (MQ=255)
                                                                                           ggAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAAAACACCTCATTCGa                                                           >  1:374920/1‑150 (MQ=255)
                                                                                            gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCTCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                          <  2:254476/150‑1 (MQ=255)
                                                                                            gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                          <  2:144220/150‑1 (MQ=255)
                                                                                            gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                          >  1:221590/1‑150 (MQ=255)
                                                                                              gaagaaGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGCCGTCAACACACCTCATTCGAGCa                                                        >  2:35488/1‑150 (MQ=255)
                                                                                              gaagaaGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCa                                                        >  1:1378194/1‑150 (MQ=255)
                                                                                              ccagaagGTACCGGCGCACGCTTGCTGTAAGATGTACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCATTTTGTAAATTTTAGGGTTTCCCCCGAGGTCAACACACCTCATTCGAGCa                                                        <  2:1314175/148‑1 (MQ=255)
                                                                                                   agGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgtgg                                                   <  1:1165038/150‑1 (MQ=255)
                                                                                                    ggTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgtggt                                                  >  1:358840/1‑150 (MQ=255)
                                                                                                        ccGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                              >  1:657732/1‑150 (MQ=255)
                                                                                                        ccGGCGCACGCGTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                              <  2:114070/150‑1 (MQ=255)
                                                                                                         cGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAgg                                             >  1:880973/1‑150 (MQ=255)
                                                                                                             gcACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGATGGTTTCCCCCGCCGTCAACACCCCTCATTCGAGCACGGGGTGGAggatac                                         >  2:266344/1‑146 (MQ=255)
                                                                                                             gcACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTcc                                         >  1:1224590/1‑150 (MQ=255)
                                                                                                                  cTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTg                                    <  2:673951/150‑1 (MQ=255)
                                                                                                                      cTGGAAGATGAACAAAACTAAAGCGGCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                <  1:227766/150‑1 (MQ=255)
                                                                                                                      cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCACCACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                >  1:1452674/1‑150 (MQ=255)
                                                                                                                      cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                <  1:1100048/150‑1 (MQ=255)
                                                                                                                      cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCGTTAGTTTGTTTTCCGGGCTTTGGCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                                <  2:1401834/150‑1 (MQ=255)
                                                                                                                        ggAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                              <  2:221590/150‑1 (MQ=255)
                                                                                                                          aaGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTAGAt                            >  2:462728/1‑150 (MQ=255)
                                                                                                                           agagGAACAACACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGACGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGGTCCGGTTGGTGTTGATg                           <  1:116534/150‑1 (MQ=255)
                                                                                                                           aGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATg                           >  2:1697388/1‑150 (MQ=255)
                                                                                                                           aGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATg                           <  1:143203/150‑1 (MQ=255)
                                                                                                                            gATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGc                          <  1:1368338/150‑1 (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)
                                                                                                                                 aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                     >  2:1497036/1‑150 (MQ=255)
                                                                                                                                 aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                     >  1:1136958/1‑150 (MQ=255)
                                                                                                                                 aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGGCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                     <  1:733520/150‑1 (MQ=255)
                                                                                                                                 aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTGTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                     >  1:74877/1‑150 (MQ=255)
                                                                                                                                   aaaaCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAtt                   >  2:1013347/1‑150 (MQ=255)
                                                                                                                                    aaaCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTa                      >  1:869855/1‑146 (MQ=255)
                                                                                                                                     aaCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGt                 >  1:934283/1‑150 (MQ=255)
                                                                                                                                      aCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa                >  1:553450/1‑150 (MQ=255)
                                                                                                                                      aCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa                <  2:1621211/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)
                                                                                                                                              gcCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGCCGTCAACACCCCTCATTCGCGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGATTTAATTGTATGTCACCg        >  1:950292/1‑150 (MQ=255)
                                                                                                                                               ccACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCgc       >  1:1133482/1‑149 (MQ=255)
                                                                                                                                                cacaAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGAc      >  2:1405618/1‑150 (MQ=255)
                                                                                                                                                    aGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGttt  >  2:1100992/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTT  >  NC_000913/4277795‑4278092

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: