Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7740GA31.1% 79.5 / 12.1 45noncoding (7/68 nt)rydBsmall regulatory RNA RydB
Reads supporting (aligned to +/- strand):  ref base G (28/3);  new base A (0/14);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 4.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAATTCAAGTTC  >  NC_000913/1764531‑1765020
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
aaCAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACACCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAAtc                                                                                                                                                                                                                                                  >  2:5399/1‑251 (MQ=255)
 aCAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAAtct                                                                                                                                                                                                                                                 <  1:702928/251‑1 (MQ=255)
       aCGCTAAGAATTAGGTATCTCATTGTTAGATAATGGATATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAACAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCCGTGGCTT‑AGCTCTTGAAGGGGCGATAAGACTAATCTCATaa                                                                                                                                                                                                                                            >  2:101353/1‑250 (MQ=255)
         gCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGGGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc                                                                                                                                                                                                                                         >  2:513233/1‑251 (MQ=255)
            aaGAATTAGGTATCTCATTGTTAGAGAATGGTTATTATCTAAGGAACATTAAGGCTGGGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGGCTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTGAGGGTTT‑CTACCCCATCCGGCGCTTATAGCCGGCACTCTCAGTGGTGT‑AGCTCTTGAAGTGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                                                                                                        <  2:493547/249‑1 (MQ=255)
              gAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGACGGGGCGATAAGAATAATCTCATAAAGCTAAc‑c                                                                                                                                                                                                                                   >  2:479001/1‑251 (MQ=255)
                              tGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAAcaca                                                                                                                                                                                                                    >  1:595634/1‑250 (MQ=255)
                               gTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACaa                                                                                                                                                                                                                   >  1:124298/1‑250 (MQ=255)
                                           ttataaTCTAGGGAGAATAAAGGCTGTGATGAAAAAGAAAAAGTTTT‑ATGAAAAATGACTATTTAAGA‑TATAGATAAAAAAGTGCAGGTTTAAAAAAAAAGAAACCCCTGAA‑GTTGGTG‑CAACCCGGGTGAAAGGTTTT‑CACCCCCACCCGGCGTTTATCCCCGGCACTCTCGGTGGTTT‑CGTTCTGGAGGGGGAGATAAAAATATTTTAATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGt                                                                                                                                                                                                       <  1:261822/250‑1 (MQ=255)
                                             attatCAAAGGGACATTAAGGTTGTGAAAAAAAAAAAAAATTGGT‑ATGAAAAATGATTATTAAAGA‑GAAAGGAAAAAAAGGGAAGCTTTAAAAAAAAAAAAAACCCTAAA‑G‑CTGGGGCGACACAGGTTTAGGGTTTT‑CTCCCCCACCCGGCGTTTATCCCCGGCCCTCTCAGGGGTTT‑AGCTTTTGAAGGGGCAAAAAAAAAAATCTCAAAAAGCAAAC‑CCGCCGTTTTAACACAAACTGCGATTAGtat                                                                                                                                                                                                     <  1:72360/250‑1 (MQ=255)
                                                      ggAGCATTAAGGCTGTGAAAAAAAAGAAAAATTTTT‑ATGGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGGGCAGCGTTAAAAAAAAAAAAACCCCTTAA‑GTCTGTG‑CGACACAGGCTAAGGGGTTT‑CTCCCCCACCCGGCGTTTATCCCCGGCCCTCTCAGTGGTTT‑AGCTCTTGAGGGGGCGAAAAAAATAATCTAATAAAGCAACC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                                                                                                           <  1:166635/251‑1 (MQ=255)
                                                        agtaTTAAGGTTGTGAT‑AAAAAAAAAAAATTGTTATTGAAAGGGATTATT‑AAGAAGAAATGAAAAAAAGTGCGGCGTTCAAAAAAAAAAAACCCCTTAT‑GTCTGGG‑CTACACAGGTTTATGGTTTC‑CCACCCCATCCGGCTTTTACCTCGGGCACTCTCA‑TGGTTTTACCTCTTGAAGGGGGAATAAAAAATATCCCATAAAGCTACC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa                                                                                                                                                                                         <  1:270300/248‑1 (MQ=255)
                                                          cATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAATTGAAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAat                                                                                                                                                                                       >  2:625572/1‑251 (MQ=255)
                                                             aaaggctggaataaaaaaaaaaaagTTGA‑ATTAAAAATAATTTTTAAAGA‑GATAGGAAAAAAAGGGCAGCGTTAAAAAATAAGAAACCCCTTAA‑GTTTGTG‑CGACCCAGGTTTAAGGTTTT‑TCACCCCACCGGGCGCTTATCCCCGGCACTCTGGGTGTCTT‑AGCTTTGGAGGGGGCGATAAAAATAATCTCAAAAATATACC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATc                                                                                                                                                                                    <  1:684985/238‑1 (MQ=255)
                                                                ggcggtgagaaaaaaaaaaaagTTGT‑TATGAAAATGATTATTTAAAA‑GAAAGGTAAAAAGGTGCAGCGTT‑AAAAAAAAAAAAACCCTTAAAGTTGGTG‑CGACAC‑GGGTTAAGGGTTTTCTCCCCCACCCGGCGTTTATCTCCGCCACTCTCAGTGGTTT‑AGCTCTTGAGGGGGCAAAAAAAATAATCTAAAAAAGCAAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAgg                                                                                                                                                                                 <  1:150639/241‑1 (MQ=255)
                                                                  cTGTGATGAAAAAGAAACAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATAAGGCg                                                                                                                                                                               >  2:214119/1‑251 (MQ=255)
                                                                  cTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCg                                                                                                                                                                               >  2:111107/1‑251 (MQ=255)
                                                                                        gT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaaga                                                                                                                                                           >  2:138998/1‑249 (MQ=255)
                                                                                               aaaaTGACTATTAAAGA‑AATAGTAAAAAAAGTGCAGCTTTCAAAAAAAAAAAACCCCTTAA‑TGCTGTG‑CGACACGGGCTGAGGGTTTT‑CTCCCCCTTCCGGCGCTTTCCCCCGGCCCTCTCGTTGGTTT‑ACTTCTTGAAGGTGCAAAAAAAATAATCTAAAAAACCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATc                                                                                                                                                   <  1:261457/251‑1 (MQ=255)
                                                                                                                            aaGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCACATCCGGCGCTTATCTCCGGCACTCTCAATGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGCATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGACGGTCATCGAACCTGGAAGGAGGTTAATCCAGGTCg                                                                                                                       >  2:419982/1‑251 (MQ=255)
                                                                                                                              gTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGAt                                                                                                                     >  2:674832/1‑251 (MQ=255)
                                                                                                                                        aaaaaaaaaaaaaccccTTAA‑GTCGGGG‑CGCCACAGGATTAAGGTTTT‑CTCCCCCCCCCGGCGCTTATCCCCGGCCCTCTCAGTGGCTT‑GGCTTTTGAGGGGGAGATAAAAATAATCTCAAAAAGCTACC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACt                                                                                                           <  1:138998/237‑1 (MQ=255)
                                                                                                                                         aGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt                                                                                                          >  1:115531/1‑251 (MQ=255)
                                                                                                                                              tAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACACCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGACTAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc                                                                                                     >  1:642156/1‑251 (MQ=255)
                                                                                                                                                 gAAAACCCTTAA‑GTATGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCAtt                                                                                                  >  2:487638/1‑251 (MQ=255)
                                                                                                                                                   aaaCCCTTAA‑GTCTGTG‑CTACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTa                                                                                                 >  2:642795/1‑250 (MQ=255)
                                                                                                                                                          aaa‑GTCTGTG‑CGACACAGGTTTAAGGTTTC‑CTCCCCCCCCCGTCCCTTCTCCCCGCCACTCTCTGTGGTTA‑ATCTCTTGAGGGGGAGAAAAGAAAAACCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAg                                                                                           <  1:172996/248‑1 (MQ=255)
                                                                                                                                                                   tg‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGc                                                                                  >  1:476257/1‑250 (MQ=255)
                                                                                                                                                                          acgggctaaggggTTT‑CTCCCCCCCCCGGCGTTTTTCCCCGGCACTCCCAGGGGTTA‑AGCTTTTGAGGGGGGAAAAAAAATAATCCCAAAAAGC‑AACACCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCa                                                                           <  1:513233/241‑1 (MQ=255)
                                                                                                                                                                           cgggattatgggTTT‑CTCCCCCCCCCGGCGCTTATCCCCGGCACTCTGAGTGGTTT‑GGCTCTGGAAGGGGAAATAAAAAAAATCTAAAAAAGC‑AACACCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAg                                                                          <  2:124298/242‑1 (MQ=255)
                                                                                                                                                                            aGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAgc                                                                         >  2:232943/1‑250 (MQ=255)
                                                                                                                                                                             ggCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGAc                                                                        >  2:41660/1‑251 (MQ=255)
                                                                                                                                                                                               cccATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGAtt                                                       >  1:462974/1‑251 (MQ=255)
                                                                                                                                                                                                     cGGCGCTTATCCCCGGCCCTCCCAGTGGTTT‑AGCTCTGGAGGGGGCAAAAAAAATAATCTAAAAAAGCTACC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTtcttct                                                 <  1:173981/251‑1 (MQ=255)
                                                                                                                                                                                                         gcttacccccGGCCCTCTGAGGGGCTT‑AGCTCTGGAAGGGGAGATAAAAATAACCTAAAAAAGCAACC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTa                                             <  1:487638/243‑1 (MQ=255)
                                                                                                                                                                                                          cTTATCTACGGCACTCTCATGGGCTT‑AGCTCTTGAAGGGGCGATACGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACATGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGAAAGGTAGATTTCTTATTTTa                                             >  2:521215/1‑250 (MQ=255)
                                                                                                                                                                                                                                     aGCTCTTGAAGGGGCGAAAAAAATATTCTAAAAATACTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGt                  <  1:111107/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCaa            >  1:645459/1‑250 (MQ=255)
                                                                                                                                                                                                                                             agggggCAATAAAAATAATCTCATAATGCAAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAAt          <  1:642795/249‑1 (MQ=255)
                                                                                                                                                                                                                                             aaGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAAt          >  1:104339/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggggCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAATTc        >  1:678304/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggggCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAATTc        >  1:325608/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gggCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAATTc        >  1:11871/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gggCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCAAGGTCGATTTGGCGAACTTGCGGCATTAAGTAAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAATTc        >  2:404967/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     aTAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAATTCAAGTTc  >  2:618990/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
AACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGT‑ATTGAAAATGACTATTTAAGA‑GATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTG‑CGACACAGGCTTAAGGGTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTT‑AGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAC‑CCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAGATTTCTTCTTTTAATTCAGCAGAAGAGAGAGTAAGCAGGTCAGCAAATTCAAGTTC  >  NC_000913/1764531‑1765020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑