Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,940 A→C 7.7% noncoding (127/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,9400AC 7.7% 398.1 / 19.0 144noncoding (127/140 nt)ryjAnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base A (66/67);  new base C (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 = 5.45e-01

Only 100 of 152 total aligned reads displayed.

CCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACG  >  NC_000913/4277791‑4278089
                                                                                                                                                     |                                                                                                                                                     
ccTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCca                                                                                                                                                       <  1:589807/150‑1 (MQ=255)
ccTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCca                                                                                                                                                       >  1:67727/1‑150 (MQ=255)
ccTAAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCca                                                                                                                                                       <  1:1064143/150‑1 (MQ=255)
 cTTCGTGGCGCTGCGTGACGAACTGGACGGCTGTATTGGGGGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCcac                                                                                                                                                      <  1:1525720/150‑1 (MQ=255)
       ggCGCTGCGTGACGAACTGGACGGATGTATTGGGTGTGGCTGCCTTTCGCGCAGTGATTTCCCGTTGCGTAACCCTGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgc                                                                                                                                                <  2:1063295/150‑1 (MQ=255)
           cTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCtt                                                                                                                                            <  2:510083/150‑1 (MQ=255)
            tGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGAGAATGCCCGTGGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCttt                                                                                                                                           <  1:1528473/150‑1 (MQ=255)
              cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGCGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAg                                                                                                                                         <  2:650965/150‑1 (MQ=255)
               gTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAgt                                                                                                                                        <  2:417649/150‑1 (MQ=255)
               gTGAAGAACTGGACTGATGTATTGGGTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgt                                                                                                                                        <  1:1489975/150‑1 (MQ=255)
                tGACGAACTTGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGGGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCTGAGCACGCTTGCTGGAATATTAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                       <  2:924047/150‑1 (MQ=255)
                tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCTCCACAAGGGCGCTTTAgtt                                                                                                                                       >  1:254476/1‑150 (MQ=255)
                tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                       >  2:1165038/1‑150 (MQ=255)
                tGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtt                                                                                                                                       <  2:328268/150‑1 (MQ=255)
                  aCGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCTTGTGGCGCTTTAgtttt                                                                                                                                     <  1:1183975/150‑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)
                              gATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCcaca                                                                                                                                                     <  1:1586466/122‑1 (MQ=255)
                                tGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATGGCCCGTGGCGTAACCCGGGCGACCGCTTAGGAGAAGACGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttg                                                                                                                       <  1:535790/150‑1 (MQ=255)
                                tGTATTGGTGGTGGCTGCCTTTCGCGCAGGGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttg                                                                                                                       <  1:134626/150‑1 (MQ=255)
                                 gTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCGTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgt                                                                                                                      >  2:778025/1‑150 (MQ=255)
                                   aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtct                                                                                                                    >  1:1621211/1‑150 (MQ=255)
                                     tGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAAAAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtcttt                                                                                                                  >  2:1236750/1‑150 (MQ=255)
                                     tGGTTGTGGCTGCCTTTCGCGCAGGGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtcttt                                                                                                                  <  1:1165021/150‑1 (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)
                                                gCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGCCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCATTAGTTAGTTTTCCGGTCTTTGTCTTTCTCTCTATccc                                                                                                       >  2:1666186/1‑150 (MQ=255)
                                                  cTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTGGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGc                                                                                                     <  2:52704/150‑1 (MQ=255)
                                                  cTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGc                                                                                                     >  2:208664/1‑150 (MQ=255)
                                                   tttCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCt                                                                                                    >  2:324083/1‑150 (MQ=255)
                                                        cgcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTac                                                                                               >  1:144220/1‑150 (MQ=255)
                                                        cgcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTac                                                                                               >  1:139057/1‑150 (MQ=255)
                                                              gATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACagga                                                                                         >  1:788327/1‑150 (MQ=255)
                                                              gATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTGTCCGGTCTTTGTCTGTCTCTCTATCCCGCTGGTACACagga                                                                                         >  1:841948/1‑150 (MQ=255)
                                                               aTTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACaggag                                                                                        <  1:59848/150‑1 (MQ=255)
                                                                  gCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCTCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGt                                                                                     >  2:1100048/1‑150 (MQ=255)
                                                                     cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTc                                                                                  >  2:127633/1‑150 (MQ=255)
                                                                     cGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTc                                                                                  <  2:504957/150‑1 (MQ=255)
                                                                       ttGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccc                                                                                <  1:775639/150‑1 (MQ=255)
                                                                        tGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTcccc                                                                               >  1:65552/1‑150 (MQ=255)
                                                                         gCGTAACCCGGGAGACCGCTTAGGAGAAGAAGGTACCTGCGCACGCTTGATGGAAGATGAACAAAACTAAAGCGCCACAATGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTccccc                                                                              <  2:350871/150‑1 (MQ=255)
                                                                           gTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGa                                                                            >  2:954048/1‑150 (MQ=255)
                                                                              aCCCGGTCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTAGCTGGAAGAAGAACAAAACTAAAGCGCAACAAGGGCGCGTTAGTTTGTTTTCCGGTCTTTCTCTTTCTCTCTATCCCGCTTGTACACAGGACGGTTTCCCCCGACGt                                                                         >  1:182344/1‑150 (MQ=255)
                                                                               cccGGGCGCCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCAGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTGAGTGTGTGTTCCGGTCTTTGGCTGTCGCTCTATCCCGCTGGGACACAGGAGGGGTTCCCCCGACGta                                                                        >  1:204655/1‑149 (MQ=255)
                                                                                 cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                      >  1:1493650/1‑150 (MQ=255)
                                                                                 cGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCaa                                                                      >  2:1368338/1‑150 (MQ=255)
                                                                                   ggCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAaca                                                                    <  1:954048/150‑1 (MQ=255)
                                                                                       aCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACAcc                                                                >  1:737484/1‑150 (MQ=255)
                                                                                       aCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACAcc                                                                >  1:657602/1‑150 (MQ=255)
                                                                                          gCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCa                                                             >  1:1668693/1‑150 (MQ=255)
                                                                                              aGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCg                                                         <  2:127832/150‑1 (MQ=255)
                                                                                               ggAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAAAACACCTCATTCGa                                                        >  1:374920/1‑150 (MQ=255)
                                                                                                gagaAGAAGGTACCGGCGCATGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACATCATTagag                                                       >  1:691801/1‑150 (MQ=255)
                                                                                                gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCTCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                       <  2:254476/150‑1 (MQ=255)
                                                                                                gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                       <  1:208664/150‑1 (MQ=255)
                                                                                                gagaAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAg                                                       <  2:144220/150‑1 (MQ=255)
                                                                                                  gaagaaGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGCCGTCAACACACCTCATTCGAGCa                                                     >  2:35488/1‑150 (MQ=255)
                                                                                                  ccagaagGTACCGGCGCACGCTTGCTGTAAGATGTACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCTCTATCCATTTTGTAAATTTTAGGGTTTCCCCCGAGGTCAACACACCTCATTCGAGCa                                                     <  2:1314175/148‑1 (MQ=255)
                                                                                                        ggTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACgtggt                                               >  1:358840/1‑150 (MQ=255)
                                                                                                          tACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGtggtgg                                             >  2:1012317/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)
                                                                                                                 gcACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTcc                                      >  1:1224590/1‑150 (MQ=255)
                                                                                                                  cACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGAGCTTTAGTTTGTTTACGGGTATTTGTCTTTCAGTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCCAGCACGTGGGGCAGGTTCCg                                     >  2:105391/1‑150 (MQ=255)
                                                                                                                      cTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTg                                 <  2:673951/150‑1 (MQ=255)
                                                                                                                          cgggAAGATGCACAAAACTCAAGCGCCACCAGGGCGCTTTTGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGGTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                             <  1:1368785/148‑1 (MQ=255)
                                                                                                                          cTGGAAGATGAACAAAACTAAAGCGGCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGgtgt                             <  1:227766/150‑1 (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                        <  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)
                                                                                                                                   gAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGATGTTCCGGTTGGTGTTGATGCttt                    <  2:1585337/150‑1 (MQ=255)
                                                                                                                                    aaCAAAACTAAAGATCAACAAGGGCGATTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTGTATCCCTCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCTAGAACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTa                   <  2:182344/150‑1 (MQ=255)
                                                                                                                                     aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                  >  1:1136958/1‑150 (MQ=255)
                                                                                                                                     aCAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTaa                  >  2:1497036/1‑150 (MQ=255)
                                                                                                                                       aaaaCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAAtt                >  2:1013347/1‑150 (MQ=255)
                                                                                                                                        aaaCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTa                   <  2:869855/146‑1 (MQ=255)
                                                                                                                                         aaCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGt              >  1:934283/1‑150 (MQ=255)
                                                                                                                                          aCTAAAGCGCCCTTGTGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa             >  1:553450/1‑150 (MQ=255)
                                                                                                                                          aCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa             <  2:1621211/150‑1 (MQ=255)
                                                                                                                                          aCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTa             <  1:854196/150‑1 (MQ=255)
                                                                                                                                          aCTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGAGGTTACGGGTGGTGGTGATTCTTTAATTGTAt            >  2:1289061/1‑150 (MQ=255)
                                                                                                                                            ttaaccgcgACTAGGGTGGTTTAGTTTGTTTTCCGGTATTTGTCTTTCTCTGTATCGGTTTTGTAAACAGTAGGGTTTCCCCCGACGTCTACACACCTCATTCGAGAACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTATTTGTATg           <  2:580994/145‑1 (MQ=255)
                                                                                                                                               aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa        >  1:1663267/1‑150 (MQ=255)
                                                                                                                                               aGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCa        <  1:1248553/150‑1 (MQ=255)
                                                                                                                                                  gcCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGCCGTCAACACCCCTCATTCGCGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGATTTAATTGTATGTCACCg     >  1:950292/1‑150 (MQ=255)
                                                                                                                                                    cacaAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGAc   >  2:1405618/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACG  >  NC_000913/4277791‑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: