Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,941 C→T 7.6% noncoding (126/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,9410CT 7.6% 431.5 / 17.6 144noncoding (126/140 nt)ryjAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base C (66/67);  new base T (3/8);  total (69/75)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

Only 100 of 149 total aligned reads displayed.

CTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACG  >  NC_000913/4277792‑4278089
                                                                                                                                                     |                                                                                                                                                    
cTTCGTGGCGCTGCGTGACGAACTGGACGGCTGTATTGGGGGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCcac                                                                                                                                                      <  1:1525720/150‑1 (MQ=255)
   aGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCCGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGTAGATGAACAAAACTAAAGCGCCACAAg                                                                                                                                                   >  2:336642/1‑150 (MQ=255)
      ggCGCTGCGTGACGAACTGGACGGATGTATTGGGTGTGGCTGCCTTTCGCGCAGTGATTTCCCGTTGCGTAACCCTGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgc                                                                                                                                                <  2:1063295/150‑1 (MQ=255)
           tGCGTGACGAACTGGACGGATGTATTGGGTGTGGCTGCCTTTCGCGCAGTGATTGCCCGGTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCttt                                                                                                                                           <  2:1338667/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:1280892/1‑150 (MQ=255)
                 ccGAACTGGAAGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGACCGTTGCGTAACCCGGGCGACCGCTTAGGAGAATAAGGTACCGGCGCACGCTTGCTGGAAGATGAACTAAACTAAAGCGCCCTTGTGGCGCTTTAgtttt                                                                                                                                     <  2:308816/149‑2 (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)
                          ccGGAAGTATTGGTTGTGGCTCCCTTTCGCGCAGGGATTGCCCGGTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTAAGTTTGTTTTCCGgt                                                                                                                            <  1:1528979/149‑1 (MQ=255)
                             gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCcaca                                                                                                                                                     <  1:1586466/122‑1 (MQ=255)
                             gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCcaca                                                                                                                                                     >  2:1586466/1‑122 (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)
                                    tGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAAAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtcttt                                                                                                                  >  2:1236750/1‑150 (MQ=255)
                                     ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGCCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCTCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtc                                                                                                                 >  2:1368785/1‑150 (MQ=255)
                                         gtgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctctc                                                                                                             <  1:1319690/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)
                                                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)
                                                            tGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACagg                                                                                          <  1:1637731/150‑1 (MQ=255)
                                                             gATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTGTCCGGTCTTTGTCTGTCTCTCTATCCCGCTGGTACACagga                                                                                         >  1:841948/1‑150 (MQ=255)
                                                              tttGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggag                                                                                        <  1:814997/149‑1 (MQ=255)
                                                              aTTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggag                                                                                        >  2:512692/1‑150 (MQ=255)
                                                              aTTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggag                                                                                        <  1:59848/150‑1 (MQ=255)
                                                              aTTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCGTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggag                                                                                        >  1:1730503/1‑150 (MQ=255)
                                                                tGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAggg                                                                                      <  2:1267271/150‑1 (MQ=255)
                                                                    cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTc                                                                                  <  2:504957/150‑1 (MQ=255)
                                                                      ttGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccc                                                                                <  1:775639/150‑1 (MQ=255)
                                                                       tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTcccc                                                                               <  2:1485812/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)
                                                                            aaCCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACg                                                                          >  2:596334/1‑150 (MQ=255)
                                                                             aCCCGGTCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTAGCTGGAAGAAGAACAAAACTAAAGCGCAACAAGGGCGCGTTAGTTTGTTTTCCGGTCTTTCTCTTTCTCTCTATCCCGCTTGTACACAGGACGGTTTCCCCCGACGt                                                                         >  1:182344/1‑150 (MQ=255)
                                                                              cccGGGCGCCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCAGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTGAGTGTGTGTTCCGGTCTTTGGCTGTCGCTCTATCCCGCTGGGACACAGGAGGGGTTCCCCCGACGta                                                                        >  1:204655/1‑149 (MQ=255)
                                                                                  ggCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAaca                                                                    <  1:954048/150‑1 (MQ=255)
                                                                                      aCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACAcc                                                                >  1:657602/1‑150 (MQ=255)
                                                                                      aCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACAcc                                                                >  1:737484/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:139057/150‑1 (MQ=255)
                                                                                               gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                       <  1:208664/150‑1 (MQ=255)
                                                                                               gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                       <  2:144220/150‑1 (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)
                                                                                                           ccGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                           <  1:944286/150‑1 (MQ=255)
                                                                                                           ccGGCGCACGCGTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAg                                           <  2:114070/150‑1 (MQ=255)
                                                                                                            cGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAgg                                          >  1:880973/1‑150 (MQ=255)
                                                                                                                gcACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTcc                                      >  1:1224590/1‑150 (MQ=255)
                                                                                                                         cTGGAAGATGAACAAAACTAAAGCGGCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                             <  1:227766/150‑1 (MQ=255)
                                                                                                                         cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCACCACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                             >  1:1452674/1‑150 (MQ=255)
                                                                                                                         cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                             <  1:1100048/150‑1 (MQ=255)
                                                                                                                         cTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGATGgtgt                             >  2:191109/1‑150 (MQ=255)
                                                                                                                           ggAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                           <  2:394281/150‑1 (MQ=255)
                                                                                                                           ggAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                           >  2:1423395/1‑150 (MQ=255)
                                                                                                                           ggAAGAAGAACAGAACTACATCGCCACAAGGGCGCTTTAGTTTGTTGTCCGGTCTTTGTCTTTCTCTCTATCCCGCTTGTACAGCGTTGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                           <  2:1514812/150‑1 (MQ=255)
                                                                                                                              agagGAACAACACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGACGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGGTCCGGTTGGTGTTGATg                        <  1:116534/150‑1 (MQ=255)
                                                                                                                              aGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATg                        <  1:143203/150‑1 (MQ=255)
                                                                                                                              aGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATg                        >  2:1697388/1‑150 (MQ=255)
                                                                                                                               gATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGc                       <  1:1368338/150‑1 (MQ=255)
                                                                                                                                 tGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCtt                     >  1:236018/1‑150 (MQ=255)
                                                                                                                                  taaCAAAACTAAATCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGTGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCttt                    <  2:886300/149‑1 (MQ=255)
                                                                                                                                    aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                  >  1:1136958/1‑150 (MQ=255)
                                                                                                                                    aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                  >  2:1497036/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)
                                                                                                                                           aaaaGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTAGACAGTATGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTTGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATg           <  1:958796/149‑1 (MQ=255)
                                                                                                                                              aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa        >  1:1663267/1‑150 (MQ=255)
                                                                                                                                              aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa        <  1:1423395/150‑1 (MQ=255)
                                                                                                                                                  ccACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCgc    >  1:1133482/1‑149 (MQ=255)
                                                                                                                                                  ccACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGa    <  2:502104/150‑1 (MQ=255)
                                                                                                                                                   cacaAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGAc   >  2:1405618/1‑150 (MQ=255)
                                                                                                                                                    acaATGTCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTGTCTATCACGCTGGTACACAGGAGGGTTTCCACCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACg  <  2:1148677/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACG  >  NC_000913/4277792‑4278089

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: