Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,808,237 G→C 6.7% intergenic (‑301/‑254) PP_4234 ← / → dsbD‑II bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,808,2370GC 6.7% 186.5 / 4.7 60intergenic (‑301/‑254)PP_4234/dsbD‑IIbifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (23/33);  new base C (2/2);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

GGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAG  >  NC_002947/4808089‑4808385
                                                                                                                                                    |                                                                                                                                                    
ccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcggg                                                                                                                                                     <  2:2213294‑M1/20‑1 (MQ=255)
  agttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                                   >  1:1319444‑M1/129‑150 (MQ=255)
   gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcgg                                                                                                                                                  <  2:1265859‑M1/23‑1 (MQ=255)
    ttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGggcggc                                                                                                                                                 <  2:1743147‑M1/24‑1 (MQ=255)
        tcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtaaaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                                             <  2:2392119‑M1/28‑1 (MQ=255)
           tgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                          <  1:2094483‑M1/31‑1 (MQ=255)
            gtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                         >  2:547092‑M1/119‑150 (MQ=255)
                     gcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGcggcgcttcgcgg                                                                                                                                        <  1:1445389‑M1/33‑7 (MQ=255)
                     gcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGcggcgcttcgcgg                                                                                                                                        >  2:1445389‑M1/110‑136 (MQ=255)
                            gaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGaaaa                                                                                                                         >  2:315778‑M1/103‑150 (MQ=255)
                                 ccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAAtcctc                                                                                                                    >  2:1386500‑M1/98‑150 (MQ=255)
                                       ggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCaa                                                                                                              <  2:1169410‑M1/59‑1 (MQ=255)
                                                     tcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                                              <  1:724621‑M1/27‑1 (MQ=255)
                                                     tcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                                              >  2:724621‑M1/78‑104 (MQ=255)
                                                     tcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaagagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTc                                                                                                <  2:1752334‑M1/73‑1 (MQ=255)
                                                     cgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                           >  2:39311‑M1/78‑107 (MQ=255)
                                                     cgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                           <  1:39311‑M1/30‑1 (MQ=255)
                                                         gtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                                                                        >  1:1454057‑M1/74‑138 (MQ=255)
                                                         gtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                                                                        <  2:1454057‑M1/65‑1 (MQ=255)
                                                                 cccaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCa                                                                                    <  2:706523‑M1/85‑1 (MQ=255)
                                                                  acaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCAc                                                                                   <  1:1387195‑M1/86‑1 (MQ=255)
                                                                  acaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCAc                                                                                   <  1:775869‑M1/86‑1 (MQ=255)
                                                                    aaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTc                                                                                 <  1:2287689‑M1/88‑1 (MQ=255)
                                                                      gcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                                          >  1:1431535‑M1/61‑87 (MQ=255)
                                                                      gcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                                          <  2:1431535‑M1/31‑5 (MQ=255)
                                                                        cggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTcaaag                                                                             >  2:1657255‑M1/59‑149 (MQ=255)
                                                                           tcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTc                                                                          <  1:1649117‑M1/95‑1 (MQ=255)
                                                                                     ctggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACgccgcc                                                                >  2:810856‑M1/46‑150 (MQ=255)
                                                                                        gtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGccccc                                                             <  1:229195‑M1/108‑1 (MQ=255)
                                                                                        gtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGccccc                                                             <  1:1266784‑M1/108‑1 (MQ=255)
                                                                                           caggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCa                                                          <  1:922031‑M1/111‑1 (MQ=255)
                                                                                                gggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAAtt                                                     <  1:1888063‑M1/116‑1 (MQ=255)
                                                                                                       gccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAacac                                                                                                          <  2:860747‑M1/63‑1 (MQ=255)
                                                                                                       gccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAacac                                                                                                          >  1:860747‑M1/28‑90 (MQ=255)
                                                                                                         cgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAACTTCCAAcgccg                                            <  2:1025871‑M1/125‑1 (MQ=255)
                                                                                                         cgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAcgccg                                            <  2:2131415‑M1/125‑1 (MQ=255)
                                                                                                             acaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTc                                        <  2:19031‑M1/129‑1 (MQ=255)
                                                                                                                aacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTcccc                                     <  2:737936‑M1/132‑1 (MQ=255)
                                                                                                                aacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTcccc                                     <  1:1489559‑M1/132‑1 (MQ=255)
                                                                                                                 acgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                                    <  2:2104940‑M1/133‑1 (MQ=255)
                                                                                                                 acgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCAATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                                    >  1:1468837‑M1/18‑150 (MQ=255)
                                                                                                                  cgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTc                                   >  2:780610‑M1/17‑150 (MQ=255)
                                                                                                                       gagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCAAACGCCGCACCCTCAAAGTTACAACGCCGCCTACCCTCGaaaa                              <  2:798397‑M1/139‑1 (MQ=255)
                                                                                                                       gagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACACGTAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGaaaa                              <  1:2331021‑M1/139‑1 (MQ=255)
                                                                                                                         gcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATc                            >  1:234265‑M1/10‑150 (MQ=255)
                                                                                                                            ggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa                         >  2:1489559‑M1/7‑150 (MQ=255)
                                                                                                                            ggtttaCCCGCGAAGCGCCGCGAGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAGCGCCGCCTCCCCTCGAAAATCcaa                         >  1:1782612‑M1/7‑150 (MQ=255)
                                                                                                                                   CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCACCAACaaa                  <  2:134819‑M1/150‑1 (MQ=255)
                                                                                                                                   CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa                  <  2:1925843‑M1/150‑1 (MQ=255)
                                                                                                                                   CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa                  <  1:72683‑M1/150‑1 (MQ=255)
                                                                                                                                   CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa                  >  1:737936‑M1/1‑150 (MQ=255)
                                                                                                                                   CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa                  <  1:742721‑M1/150‑1 (MQ=255)
                                                                                                                                    CGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAg                 <  1:894984‑M1/150‑1 (MQ=255)
                                                                                                                                       tAAGCTCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAcc              >  2:218211‑M1/2‑150 (MQ=255)
                                                                                                                                       GAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAATATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGacg              >  2:1649117‑M1/1‑149 (MQ=255)
                                                                                                                                         AGCGCCGCGCGGGCGTCGCTCAATCTCACAGGCGCTGAAAACCCTCCAGCAAACACCTAAAAACTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCt            >  2:894984‑M1/1‑150 (MQ=255)
                                                                                                                                          GCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                                    >  2:1266789‑M1/1‑125 (MQ=255)
                                                                                                                                          GCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                                    <  1:1266789‑M1/125‑1 (MQ=255)
                                                                                                                                              CGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCaaa       <  2:2035918‑M1/150‑1 (MQ=255)
                                                                                                                                               GCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAAc      >  1:1169410‑M1/1‑150 (MQ=255)
                                                                                                                                               GCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCCCAAAATTCCAACCCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAAc      <  2:336051‑M1/150‑1 (MQ=255)
                                                                                                                                                  ctGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCa   >  1:798397‑M1/3‑150 (MQ=255)
                                                                                                                                                  CGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCACCATCCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCa   >  1:1025871‑M1/1‑150 (MQ=255)
                                                                                                                                                   GGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGACCATGCCGGTCGCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGCCCCTGCAAACCCAg  >  1:1621832‑M1/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAG  >  NC_002947/4808089‑4808385

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