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.98.2% 1568.2 / 22.4 331intergenic (‑192/+27)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base . (197/127);  total (199/132)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 375 total aligned reads displayed.

CAAGGAGTGAATGTTTTCTTGGCTTGCCTGGTGGAGAGGCGAACGAACGCCTGCCCGCGCGGCAAGGGGCGACATTCTAGAGAAGGGCAACCGGCGAAGCAAGGAAGATTGGCGAATATTTGTACAGCGGATGGAATGGGCGCGCGGGGCCCCCCCCGGGGGGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTTGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGCAGTGC  >  NC_008463/5809157‑5809450
                                                                                                                                                     |                                                                                                                                                
aggcaaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGc                                                                                                                                                  >  2:9734632‑M2/143‑150 (MQ=255)
 ggcaaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggactggCGCGGGGcc                                                                                                                                                 >  2:9465320‑M2/142‑150 (MQ=255)
  gcaaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGgccg                                                                                                                                                >  2:3775617‑M2/141‑149 (MQ=255)
  gcaaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGgccg                                                                                                                                                >  1:7273052‑M2/141‑149 (MQ=255)
    aaggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGgggccgcg                                                                                                                                              >  1:9092379‑M2/139‑147 (MQ=255)
     aggagtgtatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCggggccgcgc                                                                                                                                             <  1:10032847‑M2/13‑5 (MQ=255)
     aggagtgaatgttttcttggcttgcctggtggagaggcgaacgaccgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgactatttgtacagcggatggaatggCGCggggccgcgc                                                                                                                                             >  2:5848924‑M2/138‑146 (MQ=255)
     aggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCggggccgcgc                                                                                                                                             >  1:4584027‑M2/138‑146 (MQ=255)
      ngagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGcggggccgcgcc                                                                                                                                            >  1:10236101‑M2/137‑145 (MQ=255)
      ggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGcggggccgcgcc                                                                                                                                            >  2:9045819‑M2/137‑145 (MQ=255)
      ggagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGcggggccgcgcc                                                                                                                                            >  1:1138425‑M2/137‑145 (MQ=255)
       gagtgaatgttttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCgcggggccgcgccg                                                                                                                                           <  2:824810‑M2/15‑7 (MQ=255)
       gagtgaatgttttcttggcttgcctgggggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcacggacgattggcgaatatttgtacagcggctggagtggCgcggggccgcgccg                                                                                                                                           >  1:3772016‑M2/136‑144 (MQ=255)
                  ttcttggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCg                                                                                                                               >  1:1771079‑M2/125‑150 (MQ=255)
                      tggcttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcacggaagattggcgaatctgtgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAcg                                                                                                                           >  2:8608142‑M2/121‑150 (MQ=255)
                         cttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACgcgc                                                                                                                        >  2:4613258‑M2/118‑150 (MQ=255)
                         cttgcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACgcgc                                                                                                                        >  2:5299197‑M2/118‑150 (MQ=255)
                          tagcctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctatagaagggcaaccggcgaagcaaggaagattggcgaatagttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCAcgagca                                                                                                                       >  1:1705291‑M2/117‑149 (MQ=255)
                             cctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCAc                                                                                                                    >  1:7521669‑M2/114‑150 (MQ=255)
                             cctggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCAc                                                                                                                    <  2:10674278‑M2/37‑1 (MQ=255)
                                ggtggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTcc                                                                                                                 >  2:3953090‑M2/111‑150 (MQ=255)
                                  tggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAg                                                                                                               >  1:3240519‑M2/109‑150 (MQ=255)
                                   ggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGc                                                                                                              <  2:3323409‑M2/43‑1 (MQ=255)
                                   ggagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaagatttgtagagcggatggaatggCGCGGGGCCGCGCCGGGGCGGCCCCGGACGCGCTCACTCCAGc                                                                                                              >  1:989219‑M2/108‑150 (MQ=255)
                                    gagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCg                                                                                                             >  2:6086976‑M2/107‑150 (MQ=255)
                                    gagaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCg                                                                                                             >  2:6240812‑M2/107‑150 (MQ=255)
                                     agaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCgg                                                                                                            >  1:10429988‑M2/106‑150 (MQ=255)
                                     agaggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCgg                                                                                                            >  2:6448762‑M2/106‑150 (MQ=255)
                                       aggcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggactggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCgggg                                                                                                          >  2:9241012‑M2/104‑150 (MQ=255)
                                         gcgaacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTg                                                                                                        >  1:7606756‑M2/102‑150 (MQ=255)
                                            aacgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTgcgc                                                                                                     >  1:3439521‑M2/99‑150 (MQ=255)
                                            aacgaacgccggcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatattggtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTgcgc                                                                                                     <  1:463194‑M2/52‑1 (MQ=255)
                                             acgaacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCt                                                                                                    <  2:5098482‑M2/53‑1 (MQ=255)
                                                aacgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGc                                                                                                 <  1:1454504‑M2/56‑1 (MQ=255)
                                                 acgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgactatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCg                                                                                                <  2:4698734‑M2/57‑1 (MQ=255)
                                                  cgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCgg                                                                                               <  1:1625710‑M2/58‑1 (MQ=255)
                                                  cgcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCgg                                                                                               <  2:3114683‑M2/58‑1 (MQ=255)
                                                   gcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGGgga                                                                                              >  1:6463340‑M2/92‑150 (MQ=255)
                                                   gcctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCgga                                                                                              <  1:6406662‑M2/59‑1 (MQ=255)
                                                    cctgcccgcgcgtcaaggggcgacattctagagacgggcaaccggcgaagcaaggaagattggcgactatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCggag                                                                                             <  2:364975‑M2/60‑1 (MQ=255)
                                                     ctgcccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCggtgg                                                                                            >  2:7877913‑M2/90‑150 (MQ=255)
                                                        cccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATc                                                                                         >  2:2019531‑M2/87‑150 (MQ=255)
                                                         ccgcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCt                                                                                        <  1:2579892‑M2/65‑1 (MQ=255)
                                                           gcgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGGGCGCTGGCGGAGGAtctcg                                                                                      >  2:9877913‑M2/84‑148 (MQ=255)
                                                            cgcggcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGa                                                                                     <  2:6835267‑M2/68‑1 (MQ=255)
                                                                gcaaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAg                                                                                 <  2:4312172‑M2/72‑1 (MQ=255)
                                                                  aaggggcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGc                                                                               <  1:7021169‑M2/74‑1 (MQ=255)
                                                                   aggggcgaccttctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGTGCCGCGCCGCGGCGGCCCCGCACGCGCCCACTCCAGCCGGGGGCGCTGGCGGAGGATCTCGAACAGGCa                                                                              >  1:6591433‑M2/76‑150 (MQ=255)
                                                                    ggggcgacattctagagacgggcaaccggcgaagcaaggaagattggcgaatatgtgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACCGGCAg                                                                             >  1:6099766‑M2/75‑150 (MQ=255)
                                                                       gcgacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACg                                                                          >  2:10930778‑M2/72‑150 (MQ=255)
                                                                        cgaccttctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtccagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGc                                                                         <  2:9595265‑M2/80‑1 (MQ=255)
                                                                         gacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGcc                                                                        >  2:7997169‑M2/70‑150 (MQ=255)
                                                                         gacattctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGcc                                                                        >  2:2801858‑M2/70‑150 (MQ=255)
                                                                            attctagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGt                                                                     >  1:7959115‑M2/67‑150 (MQ=255)
                                                                                tagagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCg                                                                 <  2:7250171‑M2/88‑1 (MQ=255)
                                                                                  gagaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGAc                                                               <  1:5186488‑M2/90‑1 (MQ=255)
                                                                                    gaagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCg                                                             >  1:5936337‑M2/59‑150 (MQ=255)
                                                                                     aagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCCGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCga                                                            >  1:9396347‑M2/58‑150 (MQ=255)
                                                                                     aagggcaaccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCga                                                            <  1:2440723‑M2/93‑1 (MQ=255)
                                                                                           aaccggcgacgcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCGGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGacggt                                                      >  2:9076000‑M2/52‑148 (MQ=255)
                                                                                           aaccggcgaagcaaggaagattggcgaatattcgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACgtg                                                      >  1:8900895‑M2/52‑149 (MQ=255)
                                                                                             ccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACCGGCAGACGCCGGTGGCGACCGAGACGTTCa                                                    <  2:2489938‑M2/101‑1 (MQ=255)
                                                                                             ccggcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCa                                                    >  2:9372044‑M2/50‑150 (MQ=255)
                                                                                                gcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgc                                                 >  2:6572601‑M2/47‑150 (MQ=255)
                                                                                                gcgaagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgc                                                 >  2:4895206‑M2/47‑150 (MQ=255)
                                                                                                   aagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctgc                                              >  1:8179165‑M2/44‑150 (MQ=255)
                                                                                                   aagcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctgc                                              <  1:6814727‑M2/107‑1 (MQ=255)
                                                                                                    agcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGgctgct                                             <  2:1613451‑M2/108‑1 (MQ=255)
                                                                                                    agcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTAACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACCCCGGTGGCGACCGAGACGTTCAGgctgct                                             >  1:8757202‑M2/43‑150 (MQ=255)
                                                                                                    aaGGAAGATTGGCGAATATTTGGACAGCGGATGGAATG‑‑‑GCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACACCAGCGGGGTGCGCTGGGGGGGGATCTCGACCAGGCAGACGCCGGTGGCGACCGAGACGTGCAGGCTGCTGAc                                          >  1:9716149/1‑150 (MQ=255)
                                                                                                     gcaaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGctgctg                                            >  1:848834‑M2/42‑150 (MQ=255)
                                                                                                      caaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGa                                           <  2:4251754‑M2/110‑1 (MQ=255)
                                                                                                      caaggaagatgggcgaatatttgtccagcggatggaaaggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGa                                           <  1:4375183‑M2/110‑1 (MQ=255)
                                                                                                       aaggaagcttggcgaatatttgtccagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGAc                                          <  2:6915442‑M2/111‑1 (MQ=255)
                                                                                                       aaggaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCCCACGCGCTCACACCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCCAGACGTTCAGGCTGCTGAc                                          >  1:9573457‑M2/40‑150 (MQ=255)
                                                                                                         ggacgattggcgaatattggtacagcggagggactggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGc                                        <  1:4034079‑M2/113‑1 (MQ=255)
                                                                                                          gaagattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCt                                       >  2:10241456‑M2/37‑150 (MQ=255)
                                                                                                           aagattggcgaatatttgaacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTg                                      <  1:11071165‑M2/115‑1 (MQ=255)
                                                                                                            agattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGc                                     >  2:1691876‑M2/35‑150 (MQ=255)
                                                                                                             gattggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGcc                                    >  2:5664125‑M2/34‑150 (MQ=255)
                                                                                                              attggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCCCACGCGCTCACACCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCAGCTGACGCTGCCg                                   >  1:1844830‑M2/33‑150 (MQ=255)
                                                                                                                tggcgaatatttgtaccgcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGc                                 <  2:2866973‑M2/120‑1 (MQ=255)
                                                                                                                 ggcgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCCGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGcc                                >  2:1375537‑M2/30‑150 (MQ=255)
                                                                                                                   cgaatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCAt                              <  1:663814‑M2/123‑1 (MQ=255)
                                                                                                                     atTTGTACAGCGGATGGAATG‑‑‑GCGCGGGGCCGCGCCGCGGCGGCCCCCCCCGCGCTCACACCAGCGGGGTGCGCAGGCGGAGGAGATCGAACAGGCAGACGCCGGGGGCGACCGAGACGTTCAGGCTTCTGACGCTGCCCGCCAACGGCa                         >  1:2096820/1‑150 (MQ=255)
                                                                                                                     aatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCg                            >  1:10629270‑M2/26‑150 (MQ=255)
                                                                                                                     aatatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTACAGGCTGCTGACGCTGCCGGCCATCg                            >  2:4161145‑M2/26‑150 (MQ=255)
                                                                                                                      gttGTACAGCGGATGGAATG‑‑‑GCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCCGCGGGGTGCGCGGGCGGAGGCTCTCGAACATGCAGGCGCGGGTGGCGACCGAGACGGTCAGGCTGGTGTCGCTGCCGGCCCTCGGCaa                        >  2:10150754/2‑150 (MQ=255)
                                                                                                                      atatttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCgg                           >  2:10110287‑M2/25‑150 (MQ=255)
                                                                                                                         gTACAGCGGATGGAATG‑‑‑GCGCGGGGCCGCGCCCCGGCGGCCCCCCACGCGCTCACACCAGCGGGGTGCGCTGGCGGTGGATCTCGAACAGGCAGAAGCCGGTGGCGACCCAGAAGTTCAGGCTGCTGAAGCTGCCGGCCATAGGCAACCt                     >  1:3678925/1‑150 (MQ=255)
                                                                                                                          ttgtacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAAc                       <  1:2930627‑M2/130‑1 (MQ=255)
                                                                                                                             tacagcggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTg                    >  2:8070356‑M2/18‑150 (MQ=255)
                                                                                                                              acagcggatggactggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTgg                   >  2:11082624‑M2/17‑150 (MQ=255)
                                                                                                                              acagcggatggactggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTgg                   <  1:4188095‑M2/134‑1 (MQ=255)
                                                                                                                                  cggatggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAg               <  2:10017782‑M2/138‑1 (MQ=255)
                                                                                                                                      tggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCATCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGTAg           >  2:1830874‑M2/9‑150 (MQ=255)
                                                                                                                                       ggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGACGt          >  2:5013191‑M2/8‑150 (MQ=255)
                                                                                                                                       ggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGt          <  2:5899709‑M2/143‑1 (MQ=255)
                                                                                                                                       ggaatggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCAGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGt          >  2:7599570‑M2/8‑150 (MQ=255)
                                                                                                                                          atggCGCGGGGCCGCGCCGCGGCGGCCCCGCACGCGCTCACTCCAGCGGGGTGCGCTGGCGGAGGATCTCGAACAGGCAGACGCCGGTGGCGACCGAGACGTTCAGGCTGCTGACGCTGCCGGCCATCGGCAACCTGGCGAGGAAGTCGc       >  1:4506376‑M2/5‑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: