Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,809,308 Δ1 bp intergenic (‑194/+25) rpsF ← / ← rlmB 30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,809,3060C.97.4% 1640.9 / 35.9 352intergenic (‑192/+27)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base C (7/2);  new base . (187/156);  total (194/158)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 383 total aligned reads displayed.

CAAGGAGTGAATGTTTTCTTGGCTTGCCTGGTGGAGAGGCGAACGAACGCCTGCCCGCGCGGCAAGGGGCGACATTCTAGAGAAGGGCAACCGGCGAAGCAAGGAAGATTGGCGAATATTTGTACAGCGGATGGAATGGGCGCGCGGGGCCCCCCCCGGGGGGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTTGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTGC  >  NC_008463/5809157‑5809450
                                                                                                                                                     |                                                                                                                                                
aggcaaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGc                                                                                                                                                  <  2:10131432‑M2/8‑1 (MQ=255)
aggcaaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGc                                                                                                                                                  <  2:6877910‑M2/8‑1 (MQ=255)
  gcaaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGgccg                                                                                                                                                <  2:6450948‑M2/10‑2 (MQ=255)
   caaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGggccgc                                                                                                                                               <  2:4789088‑M2/11‑3 (MQ=255)
    aaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagtcgggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGgggccgcg                                                                                                                                              >  2:10203814‑M2/139‑147 (MQ=255)
     aggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCggggccgcgc                                                                                                                                             >  1:336963‑M2/138‑146 (MQ=255)
     aggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCggggccgcgc                                                                                                                                             >  2:3716842‑M2/138‑146 (MQ=255)
     aggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatgtgtacagcggatggaatggCGCggggccgcgc                                                                                                                                             >  1:4876378‑M2/138‑146 (MQ=255)
     aggagtgaatgatttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggctcggggcgacattctagagaagggcacccggcgaagcacggacgattggcgaatatttgtcaagcggatggaatggCGCggggccgcgc                                                                                                                                             >  1:6266748‑M2/138‑146 (MQ=255)
      ggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGcggggccgcgcc                                                                                                                                            >  2:6593415‑M2/137‑145 (MQ=255)
      ggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatagttgtacagcggatggaatggCGcggggccgcgcc                                                                                                                                            >  2:6084758‑M2/137‑145 (MQ=255)
      ggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaagatttgtacagcggatggaatggCGcggggccgcgcc                                                                                                                                            >  1:5886687‑M2/137‑145 (MQ=255)
       gagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCgcggggccgcgccg                                                                                                                                           >  1:6134141‑M2/136‑144 (MQ=255)
        agtgaatgttttcttggcttgcctggtggagaggcgaccgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggcgcggggccgcgccgc                                                                                                                                          <  2:10871581‑M2/16‑8 (MQ=255)
                   tcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggacgtttggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGc                                                                                                                              <  1:3440260‑M2/27‑1 (MQ=255)
                      tggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAcg                                                                                                                           >  2:7714335‑M2/121‑150 (MQ=255)
                       ggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAcgc                                                                                                                          >  1:1648281‑M2/120‑150 (MQ=255)
                       ggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAcgc                                                                                                                          <  2:8615553‑M2/31‑1 (MQ=255)
                       ggctggcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagccaaggccaaccggcgaagcacggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAcgc                                                                                                                          <  2:1907578‑M2/31‑1 (MQ=255)
                         cttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACgcgc                                                                                                                        >  1:7375094‑M2/118‑150 (MQ=255)
                           ggcctggtggagaggcgaacgaacgcctgcccgcgcggcatggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTc                                                                                                                      <  2:5457241‑M2/35‑1 (MQ=255)
                            gcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCa                                                                                                                     >  1:3029299‑M2/115‑150 (MQ=255)
                             cctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagacgggcaaccggcgaagcacggaagattggcgaatatttgtacagcggagggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCAc                                                                                                                    <  2:9259494‑M2/37‑1 (MQ=255)
                              ctggtggagaggcgaacgaccgcctgcccgcgcggcacggggcgacattctagcgaagggcaaccggcgaagcacggaagtttggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACt                                                                                                                   <  1:2014078‑M2/38‑1 (MQ=255)
                                ggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTcc                                                                                                                 <  1:10144046‑M2/40‑1 (MQ=255)
                                 gtggagaggcgaccgaacgcctgcccgcgcggcacggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatgtgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCGCACTCCa                                                                                                                >  1:10915571‑M2/110‑150 (MQ=255)
                                 gtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCa                                                                                                                >  1:2163058‑M2/110‑150 (MQ=255)
                                    gagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCg                                                                                                             <  2:5084018‑M2/44‑1 (MQ=255)
                                        ggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgactatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGt                                                                                                         >  1:5116647‑M2/103‑150 (MQ=255)
                                            aacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTgcgc                                                                                                     >  2:9030164‑M2/99‑150 (MQ=255)
                                             acgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggacgattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCt                                                                                                    >  2:7910811‑M2/98‑150 (MQ=255)
                                              cgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTg                                                                                                   <  2:4978278‑M2/54‑1 (MQ=255)
                                               gaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTgg                                                                                                  >  1:55388‑M2/96‑150 (MQ=255)
                                                        cccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATc                                                                                         >  2:3526916‑M2/87‑150 (MQ=255)
                                                           gcgcggcaaggggcgacattctagagaagggcaaccggcgacgcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGAtctcg                                                                                      <  1:6400201‑M2/67‑3 (MQ=255)
                                                             gcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGaa                                                                                    <  1:4491575‑M2/69‑1 (MQ=255)
                                                              cggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAAc                                                                                   <  1:6514845‑M2/70‑1 (MQ=255)
                                                               ggcaaggggcgacattctagagaagggccaccggcgaagcaaggaagattggcgaatctttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACa                                                                                  <  1:689806‑M2/71‑1 (MQ=255)
                                                                   aggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCa                                                                              >  2:5424837‑M2/76‑150 (MQ=255)
                                                                   aggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCa                                                                              >  2:6228354‑M2/76‑150 (MQ=255)
                                                                    ggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAggccg                                                                             >  1:10320420‑M2/75‑148 (MQ=255)
                                                                     cGAGATTCTAGAGCAGGTCAACCGGCGACGCAAGGAAGATTGGCGAATATTTCTACAGCGGATGGAATG‑‑‑GCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGGTGGCGGACGAACTCGAACAGGCAGACGc                                                                         <  2:8460933/150‑1 (MQ=255)
                                                                           cattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTAACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCgg                                                                      >  1:9145789‑M2/68‑150 (MQ=255)
                                                                           cattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggagggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCgg                                                                      >  2:5957428‑M2/68‑150 (MQ=255)
                                                                            attctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtaccgcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGt                                                                     <  2:6767284‑M2/84‑1 (MQ=255)
                                                                            attctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGGAGACGCCGGt                                                                     >  1:4541071‑M2/67‑150 (MQ=255)
                                                                            attctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGt                                                                     >  1:4162775‑M2/67‑150 (MQ=255)
                                                                            attctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGt                                                                     <  2:3010271‑M2/84‑1 (MQ=255)
                                                                              tctagagcagggcaaccggcgacgcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTgg                                                                   <  2:4560587‑M2/86‑1 (MQ=255)
                                                                               ctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGc                                                                  >  2:3244506‑M2/64‑150 (MQ=255)
                                                                                tagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCg                                                                 >  1:10012678‑M2/63‑150 (MQ=255)
                                                                                  gagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGAc                                                               >  1:324107‑M2/61‑150 (MQ=255)
                                                                                   agaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGAcc                                                              >  2:6929176‑M2/60‑150 (MQ=255)
                                                                                    gaagggcaaccggcgaagcaaggaagattggcgcatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACCGGCCGACGCCGGTGGCGACCg                                                             >  1:5291654‑M2/59‑150 (MQ=255)
                                                                                         gcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACg                                                        <  2:11018024‑M2/97‑1 (MQ=255)
                                                                                            accggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGTGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTc                                                     >  1:8717873‑M2/51‑150 (MQ=255)
                                                                                            accggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTc                                                     >  1:1329046‑M2/51‑150 (MQ=255)
                                                                                               ggcgaagcaaggacgattggcgaatatttgtacagcggatggaagggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAgg                                                  <  1:806155‑M2/103‑1 (MQ=255)
                                                                                                 cgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAggcg                                                >  2:9367555‑M2/46‑149 (MQ=255)
                                                                                                   aagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctgc                                              <  1:1909886‑M2/107‑1 (MQ=255)
                                                                                                   aagcaaggaagattggcgaatatttgtacagcggatggaagggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctgc                                              <  1:7501186‑M2/107‑1 (MQ=255)
                                                                                                    agcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCCCGCGCTCACTCCAGCGGGGTGCGCGGGCGGGGGATCTCGAACAGGCAGACGCCGGGGGCGACCGAGACGTTCAggctgcg                                             >  1:3859782‑M2/43‑149 (MQ=255)
                                                                                                     gcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCCCTCACACCAGCGGGGTGCGCGGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCacgctgcag                                            >  1:2796500‑M2/42‑148 (MQ=255)
                                                                                                        aggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACg                                         >  1:5982263‑M2/39‑150 (MQ=255)
                                                                                                        aggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACg                                         <  1:2357589‑M2/112‑1 (MQ=255)
                                                                                                          gaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCt                                       <  2:1652439‑M2/114‑1 (MQ=255)
                                                                                                             gattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGctgac                                    <  1:3117764‑M2/117‑3 (MQ=255)
                                                                                                              attggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCg                                   <  2:4700872‑M2/118‑1 (MQ=255)
                                                                                                               ttggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGACgg                                  >  2:2326915‑M2/32‑150 (MQ=255)
                                                                                                                tggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGc                                 >  2:7517599‑M2/31‑150 (MQ=255)
                                                                                                                 ggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGcc                                >  2:7114853‑M2/30‑150 (MQ=255)
                                                                                                                   cgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCCCGGCGGCCCCGCACGCGCTCACACCAGCGGGGGGCGCGGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCggccct                              >  1:3708690‑M2/28‑148 (MQ=255)
                                                                                                                    gaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATc                             <  2:3803857‑M2/124‑1 (MQ=255)
                                                                                                                     aatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAAGCGCTCACACCAGCGGGGTGCGCTGGGGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCCTCg                            >  1:10473930‑M2/26‑150 (MQ=255)
                                                                                                                      ctatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCgg                           <  2:4025687‑M2/126‑1 (MQ=255)
                                                                                                                       tatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCCCACGCGCGCACACCAGCGGGGTGCGCGGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGc                          >  1:9229612‑M2/24‑150 (MQ=255)
                                                                                                                         tttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCaa                        <  1:1228305‑M2/129‑1 (MQ=255)
                                                                                                                          ttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAAc                       <  2:5670258‑M2/130‑1 (MQ=255)
                                                                                                                          ttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAAc                       >  2:10616315‑M2/21‑150 (MQ=255)
                                                                                                                           tgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAAcc                      >  2:8074543‑M2/20‑150 (MQ=255)
                                                                                                                            gtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCt                     <  2:2893338‑M2/132‑1 (MQ=255)
                                                                                                                            gtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCt                     <  2:834654‑M2/132‑1 (MQ=255)
                                                                                                                             tacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTTCTGACGCTGCCGGCCATCGGCaaccag                    >  2:2011090‑M2/18‑148 (MQ=255)
                                                                                                                              acagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTgg                   >  2:8126648‑M2/17‑150 (MQ=255)
                                                                                                                               cagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGc                  <  1:244579‑M2/135‑1 (MQ=255)
                                                                                                                               cagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACACCCGCGGGGTGCGCTGGCGGAGGATCTCGCAAAGGCAGACGCCGGTGGCGACCGAGCCGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGc                  >  1:10754802‑M2/16‑150 (MQ=255)
                                                                                                                                agcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCg                 <  2:10326710‑M2/136‑1 (MQ=255)
                                                                                                                                    gatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGa             <  1:1125548‑M2/140‑1 (MQ=255)
                                                                                                                                     atggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCAGGCGAGGaa            >  2:4336981‑M2/10‑150 (MQ=255)
                                                                                                                                       ggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGt          <  1:11077188‑M2/143‑1 (MQ=255)
                                                                                                                                        gaatggcgcgGGGCCGCGCCGCGGCGGCCCCCCACGCGCACACACCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCCGTGGCGACCGAGACCTTCAGGCTGCTGACGCAGCCGGCCAACGGCAACCTGGCGAGGAAGTc         >  1:7725020/6‑150 (MQ=255)
                                                                                                                                        gaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGACGGTGGCGACCGAGACGTTCAGGCTGCGGACGCTGCCGGCCAAAGGCAACCTGGCGAGGAAGTc         >  2:1387090‑M2/7‑150 (MQ=255)
                                                                                                                                         aatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCg        <  1:8531809‑M2/145‑1 (MQ=255)
                                                                                                                                          atggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGc       <  2:3302231‑M2/146‑1 (MQ=255)
                                                                                                                                           tggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCa      >  1:7005652‑M2/4‑150 (MQ=255)
                                                                                                                                            ggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAg     >  2:8730346‑M2/3‑150 (MQ=255)
                                                                                                                                             gCGCGCGGCCGCCCCGCGGCGGCCCCGCACGCGCTCACACCAGCCGGCTGCGCTCGCGGAGGATCACGCACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACCCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGt    <  1:825171‑M2/149‑1 (MQ=255)
                                                                                                                                              cgcgGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTg   >  1:9512304/1‑150 (MQ=255)
                                                                                                                                               gcgGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTGc  >  2:2745897/1‑150 (MQ=255)
                                                                                                                                               gcgGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTGc  >  1:2684608/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CAAGGAGTGAATGTTTTCTTGGCTTGCCTGGTGGAGAGGCGAACGAACGCCTGCCCGCGCGGCAAGGGGCGACATTCTAGAGAAGGGCAACCGGCGAAGCAAGGAAGATTGGCGAATATTTGTACAGCGGATGGAATGGGCGCGCGGGGCCCCCCCCGGGGGGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTTGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTGC  >  NC_008463/5809157‑5809450

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: