Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,809,311 C→G intergenic (‑197/+22) 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,3110CG97.7% 1414.6 / 21.1 394intergenic (‑197/+22)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base G (224/159);  minor base . (7/2);  total (233/161)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

Only 100 of 394 total aligned reads displayed.

CTTGGCTTGCCTGGTGGAGAGGCGAACGAACGCCTGCCCGCGCGGCAAGGGGCGACATTCTAGAGAAGGGCAACCGGCGAAGCAAGGAAGATTGGCGAATATTTGTACAGCGGATGGAATGGGCGCGCGGGGCCCCCCCCGGGGGGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTTGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTGCT  >  NC_008463/5809174‑5809451
                                                                                                                                         |                                                                                                                                            
tttcttggcttgcctggtggagaggcgaacgaccgcctgcccgcgcggcaaggggcgacattctagagacgggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCg                                                                                                                                <  2:2909572‑M2/26‑1 (MQ=255)
     tgtcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgagattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAcgc                                                                                                                           >  1:11356109‑M2/121‑151 (MQ=255)
      ggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAcgcg                                                                                                                          <  2:7251160‑M2/32‑1 (MQ=255)
        cttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCt                                                                                                                        >  1:6881743‑M2/118‑151 (MQ=255)
         ttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTc                                                                                                                       <  2:11081478‑M2/35‑1 (MQ=255)
             ctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCAATc                                                                                                                   >  1:89111‑M2/113‑151 (MQ=255)
              tggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtaccgcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTcc                                                                                                                  <  1:2326073‑M2/40‑1 (MQ=255)
                 tggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagacgggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGc                                                                                                               <  2:7876791‑M2/43‑1 (MQ=255)
                     gaggcgaacgaccgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcacggaagattggcgaatatttgtacagcggatggactggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCgggg                                                                                                           <  2:8618841‑M2/47‑1 (MQ=255)
                     gaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCgggg                                                                                                           >  1:5254362‑M2/105‑151 (MQ=255)
                       ggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTg                                                                                                         <  2:2390970‑M2/49‑1 (MQ=255)
                         cgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTgcg                                                                                                       >  1:9849056‑M2/101‑151 (MQ=255)
                           aacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCt                                                                                                     <  2:3563348‑M2/53‑1 (MQ=255)
                           aacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCt                                                                                                     >  1:11222771‑M2/99‑151 (MQ=255)
                              gaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGc                                                                                                  >  1:8049472‑M2/96‑151 (MQ=255)
                              gaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGc                                                                                                  >  1:6520619‑M2/96‑151 (MQ=255)
                              gaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGc                                                                                                  >  1:5675431‑M2/96‑151 (MQ=255)
                               aacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCg                                                                                                 >  2:7622226‑M2/95‑151 (MQ=255)
                                 cgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCgga                                                                                               >  1:10961888‑M2/93‑151 (MQ=255)
                                  gcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACACCAGCGGGGTGCGCTGGCggag                                                                                              >  2:11527413‑M2/92‑151 (MQ=255)
                                  gccggcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCggag                                                                                              <  2:11725176‑M2/60‑1 (MQ=255)
                                    ctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCCGCGGGGTGCGCTGGCggagga                                                                                            >  2:5207107‑M2/90‑151 (MQ=255)
                                    ctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGATGGCggagga                                                                                            >  1:6731244‑M2/90‑151 (MQ=255)
                                     tgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGAt                                                                                           >  2:1256801‑M2/89‑151 (MQ=255)
                                     ggcccgcgcggcaaggggcgacattctagataagggcacccggcgaatcaaggaagtttggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGAt                                                                                           <  1:8886124‑M2/63‑1 (MQ=255)
                                        ccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATctc                                                                                        <  2:4428076‑M2/66‑2 (MQ=255)
                                        ccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATctc                                                                                        >  1:955886‑M2/86‑150 (MQ=255)
                                         cgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGAtctcg                                                                                       <  1:4613051‑M2/67‑3 (MQ=255)
                                           cgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGaa                                                                                     <  1:11632226‑M2/69‑1 (MQ=255)
                                            gcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAAc                                                                                    <  1:3732174‑M2/70‑1 (MQ=255)
                                                 aaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCa                                                                               >  1:10848484‑M2/77‑151 (MQ=255)
                                                  aggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAggccg                                                                              >  1:11050949‑M2/76‑149 (MQ=255)
                                                  aggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAg                                                                              >  2:8215122‑M2/76‑151 (MQ=255)
                                                    gggcgacattctagagaagggccaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCACGAACAGGCAGAc                                                                            <  2:4309700‑M2/78‑1 (MQ=255)
                                                    gggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGAc                                                                            <  2:6749564‑M2/78‑1 (MQ=255)
                                                    gggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGAc                                                                            >  2:9274804‑M2/74‑151 (MQ=255)
                                                        gacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggttggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCg                                                                        <  2:1387827‑M2/82‑1 (MQ=255)
                                                           attctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGTTCTCGAACAGGCAGACGCCGGTg                                                                     >  1:10621132‑M2/67‑151 (MQ=255)
                                                           attctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTg                                                                     <  1:7351293‑M2/85‑1 (MQ=255)
                                                             tctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGc                                                                   <  2:7122422‑M2/87‑1 (MQ=255)
                                                               tagagacgggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggagggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGa                                                                 <  1:6113221‑M2/89‑1 (MQ=255)
                                                               tagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGa                                                                 >  2:9041031‑M2/63‑151 (MQ=255)
                                                                  agaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCg                                                              >  2:6363677‑M2/60‑151 (MQ=255)
                                                                     cgggcaaccggcgacgcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCgaga                                                           <  1:389532‑M2/95‑1 (MQ=255)
                                                                     agggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCgaga                                                           <  2:9765094‑M2/95‑1 (MQ=255)
                                                                      gggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGAc                                                          >  2:5485951‑M2/56‑151 (MQ=255)
                                                                       taCCGGCGAAGCAAGGTAGATTGGCGAATGTTTGTCCAGCGGATGGAATG‑‑‑GCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGCTGTCGAACAGGCAGTCGGCGGTGGCGGCCGCGACGtgc                                                      >  2:9138881/2‑149 (MQ=255)
                                                                         taaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCGCACTCCAGCGGGGTGCGCTGGCGGAGGCTCTCGAACAGGCAGACGCCGGTGGCGACCGAGacggt                                                       >  1:4969863‑M2/53‑149 (MQ=255)
                                                                         caaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCGCGAACAGGCAGACGCCGGTGGCGACCGAGACGtt                                                       >  2:10010073‑M2/53‑151 (MQ=255)
                                                                          aaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTc                                                      >  2:5436871‑M2/52‑151 (MQ=255)
                                                                           cccggcgaagcaaggaagattggcgaatatttgtaccgcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCa                                                     <  1:3339053‑M2/101‑1 (MQ=255)
                                                                            ccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAg                                                    >  2:10358277‑M2/50‑151 (MQ=255)
                                                                              ggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgc                                                  >  1:3281512‑M2/48‑151 (MQ=255)
                                                                              ggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgc                                                  <  2:7961100‑M2/104‑1 (MQ=255)
                                                                               gcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgct                                                 <  1:6645808‑M2/105‑1 (MQ=255)
                                                                               gcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGGTCAGgct                                                 >  2:11000112‑M2/47‑151 (MQ=255)
                                                                                cgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACACCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctg                                                >  2:2018428‑M2/46‑151 (MQ=255)
                                                                                 gaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctgc                                               >  2:1233401‑M2/45‑151 (MQ=255)
                                                                                 gaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctgc                                               >  2:740427‑M2/45‑151 (MQ=255)
                                                                                 gaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctgc                                               >  1:4381502‑M2/45‑151 (MQ=255)
                                                                                   agcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGctgctg                                             >  2:6589730‑M2/43‑151 (MQ=255)
                                                                                    gcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGCCGCCGGTGGCGACCGAGACGTTCAGGCTGCTGa                                            >  2:10795058‑M2/42‑151 (MQ=255)
                                                                                    gcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGa                                            <  1:1610868‑M2/110‑1 (MQ=255)
                                                                                      aaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACg                                          >  2:5655794‑M2/40‑151 (MQ=255)
                                                                                      aaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACg                                          <  2:4929401‑M2/112‑1 (MQ=255)
                                                                                       aggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTACAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGc                                         >  1:9126747‑M2/39‑151 (MQ=255)
                                                                                        ggaagattggcgaatatttggacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCt                                        <  2:2642865‑M2/114‑1 (MQ=255)
                                                                                          aagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCGCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGc                                      >  1:2234381‑M2/36‑151 (MQ=255)
                                                                                           agattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGGGGATCTCGAACAGGCACACGCCGGTGGCGACCGAGACGGTAAGGCTGCTGACGCTGcc                                     >  1:7864540‑M2/35‑151 (MQ=255)
                                                                                           agattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCCGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGcc                                     >  1:6399259‑M2/35‑151 (MQ=255)
                                                                                           agattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGcc                                     >  2:10469848‑M2/35‑151 (MQ=255)
                                                                                           agattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGcc                                     >  2:6123045‑M2/35‑151 (MQ=255)
                                                                                            gattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCACGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCg                                    <  1:2970504‑M2/118‑1 (MQ=255)
                                                                                               tggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCCGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGcc                                 >  2:6885409‑M2/31‑151 (MQ=255)
                                                                                                ggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCa                                >  2:9528050‑M2/30‑151 (MQ=255)
                                                                                                 gcgaatctatgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGCCGTTCAGGCTGCTGACGCTGCCGGCCAt                               <  1:5616288‑M2/123‑1 (MQ=255)
                                                                                                 gcgaatatttgtacagcggatggtatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCAt                               <  1:2036288‑M2/123‑1 (MQ=255)
                                                                                                 gcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCAt                               <  2:4706324‑M2/123‑1 (MQ=255)
                                                                                                  cgaatatttgtacagcggatggactggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGccagc                              >  1:1241932‑M2/28‑149 (MQ=255)
                                                                                                    aatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCgg                            <  1:10120632‑M2/126‑1 (MQ=255)
                                                                                                     tttGTTCAGCGGATGGTCTG‑‑‑GCGCGGGGCCGCGCCGGGGCGGCGCCGCACGCGGTGACTCCAGCGGGGTGCGCTGGCGGTGGATCGGGAAGTGGGCGTCGCCGGTGGCGGCCGTGTCGTTCAGGGTGCGGTCGGTGCCGGGCAGCGGCacc                        >  2:10116828/1‑149 (MQ=255)
                                                                                                     atatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGc                           <  2:655036‑M2/127‑1 (MQ=255)
                                                                                                        tttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAAc                        <  2:5547760‑M2/130‑1 (MQ=255)
                                                                                                           gtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTg                     >  1:9648748‑M2/19‑151 (MQ=255)
                                                                                                            tacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGGGGAGGATCTCGAACAGGCAGAAGCCGGTGGCGACCGAGAAGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCAgg                    >  1:8324596‑M2/18‑151 (MQ=255)
                                                                                                            tacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTgg                    >  2:3172367‑M2/18‑151 (MQ=255)
                                                                                                            aGCGGATGGAATG‑‑‑GCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACGCCAGCGGGGTGCGCTGGCGGAGGTTCGCGTACCGGCAGACGCCGGTGGCGACCGTGACGTTCAGGCTGCTGTCGCTGCCGGCCAGCGGCAACCTGGCGa                 >  1:10347976/1‑151 (MQ=255)
                                                                                                              cagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACCGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCg                  <  2:3853140‑M2/136‑1 (MQ=255)
                                                                                                              cagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGAATACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCTACCTGGCg                  <  2:8405800‑M2/136‑1 (MQ=255)
                                                                                                               agcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGa                 >  2:3663780‑M2/15‑151 (MQ=255)
                                                                                                                     tggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTGCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGCAGt           >  2:1874867‑M2/9‑151 (MQ=255)
                                                                                                                        aatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGc        >  2:9211447‑M2/6‑151 (MQ=255)
                                                                                                                         atggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCa       >  2:2220635‑M2/5‑151 (MQ=255)
                                                                                                                          tggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAg      >  1:6935294‑M2/4‑151 (MQ=255)
                                                                                                                          tggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAg      >  1:7161048‑M2/4‑151 (MQ=255)
                                                                                                                          tggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAg      >  2:1653903‑M2/4‑151 (MQ=255)
                                                                                                                           ggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGt     <  1:7279559‑M2/149‑1 (MQ=255)
                                                                                                                           ggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGt     >  2:8314054‑M2/3‑151 (MQ=255)
                                                                                                                              gcgGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCTGTGCt  >  1:3138182/1‑151 (MQ=255)
                                                                                                                              gcgGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTGCt  <  2:3320158/151‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
CTTGGCTTGCCTGGTGGAGAGGCGAACGAACGCCTGCCCGCGCGGCAAGGGGCGACATTCTAGAGAAGGGCAACCGGCGAAGCAAGGAAGATTGGCGAATATTTGTACAGCGGATGGAATGGGCGCGCGGGGCCCCCCCCGGGGGGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTTGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTGCT  >  NC_008463/5809174‑5809451

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: