Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,674 G→C 6.4% P311R (CCG→CGG)  mutY ← A/G‑specific adenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,6740GC 6.4% 195.3 / 4.1 62P311R (CCG→CGG) mutYA/G‑specific adenine glycosylase
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base C (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

ATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGG  >  NC_002947/344525‑344819
                                                                                                                                                     |                                                                                                                                                 
aTCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCg                                                                                                                                                   >  1:1347646/1‑150 (MQ=255)
 tCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCgg                                                                                                                                                  >  1:178539/1‑150 (MQ=255)
   aCGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGt                                                                                                                                                >  2:514557/1‑150 (MQ=255)
   aCGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCGTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGt                                                                                                                                                >  1:533500/1‑150 (MQ=255)
      gTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGa                                                                                                                                             <  1:405502/150‑1 (MQ=255)
       tGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGAc                                                                                                                                            <  2:1113006/150‑1 (MQ=255)
                 cGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCACGTCGACCCGGTCGGCGAGCACgtgg                                                                                                                                  <  1:546680/150‑4 (MQ=255)
                     cGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACgtggccga                                                                                                                              <  2:1347646/150‑8 (MQ=255)
                       tctcCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGc                                                                                                                            <  1:160831/150‑1 (MQ=255)
                             gCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATc                                                                                                                      <  1:815082/150‑1 (MQ=255)
                                    gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCt                                                                                                               >  1:939339/1‑150 (MQ=255)
                                    gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCt                                                                                                               >  1:1427287/1‑150 (MQ=255)
                                    gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCt                                                                                                               >  2:1211818/1‑150 (MQ=255)
                                    gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCAATCGCCAGCt                                                                                                               >  1:965535/1‑150 (MQ=255)
                                    gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCCCGATCGCCAGCt                                                                                                               >  2:59084/1‑150 (MQ=255)
                                     ttCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTg                                                                                                              >  1:519694/1‑150 (MQ=255)
                                       cgtcgGCCGCGCGTTCAAGCAGC‑TTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGaa                                                                                                           <  1:1202270/150‑1 (MQ=255)
                                         tcgGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAg                                                                                                          <  1:41046/150‑1 (MQ=255)
                                               gcgcgTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTg                                                                                                    >  2:1204774/1‑150 (MQ=255)
                                               gcgcgTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGCCTCGATCGCCAGCTGGAAGTGCCTg                                                                                                    >  2:610591/1‑150 (MQ=255)
                                                   gTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGa                                                                                                                                    >  1:436551/1‑114 (MQ=255)
                                                   gTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGa                                                                                                                                    <  2:436551/114‑1 (MQ=255)
                                                   gTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                       >  2:633212/1‑143 (MQ=255)
                                                   gTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                       <  1:633212/143‑1 (MQ=255)
                                                      cAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTAt                                                                                             <  1:514557/150‑1 (MQ=255)
                                                       aaGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATg                                                                                            >  1:1394540/1‑150 (MQ=255)
                                                       aaGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCACGTCGACCCGGTCggcga                                                                                                                                          >  1:1144244/1‑100 (MQ=255)
                                                       aaGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCACGTCGACCCGGTCGGCga                                                                                                                                          <  2:1144244/104‑2 (MQ=255)
                                                        agcagcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATgg                                                                                           <  2:1427287/150‑1 (MQ=255)
                                                          cagcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGt                                                                                         <  1:1323050/150‑1 (MQ=255)
                                                          cagcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGt                                                                                         >  2:1225407/1‑150 (MQ=255)
                                                             cTTTTTCCCCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAg                                                                                      <  2:388220/150‑1 (MQ=255)
                                                                    aCCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACAAGAAACCGCCGAAGGCAGCCGCGCGCGATCTGGGCGGTTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTcc                                                                               >  1:1242824/1‑150 (MQ=255)
                                                                        gggCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCAtgg                                                                           >  1:788427/1‑150 (MQ=255)
                                                                                     cAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                       >  2:390268/1‑109 (MQ=255)
                                                                                     cAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                       <  1:390268/109‑1 (MQ=255)
                                                                                      aGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGt                                                                                                                                                <  1:1134885/67‑1 (MQ=255)
                                                                                      aGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGt                                                                                                                                                >  2:1134885/1‑67 (MQ=255)
                                                                                      aGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGgccagcc                                                             <  2:939339/150‑1 (MQ=255)
                                                                                              cggGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAg                                                     <  1:1225407/150‑1 (MQ=255)
                                                                                                     gCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGc                                              >  2:150292/1‑150 (MQ=255)
                                                                                                        aGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTgg                                           <  2:178539/150‑1 (MQ=255)
                                                                                                           ttATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATgg                                                                                           >  2:430268/1‑99 (MQ=255)
                                                                                                           ttATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATgg                                                                                           <  1:430268/99‑1 (MQ=255)
                                                                                                                ccAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAgg                                       >  1:915557/1‑146 (MQ=255)
                                                                                                                ccAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAgg                                       <  2:915557/146‑1 (MQ=255)
                                                                                                                 cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                       <  1:1030670/81‑1 (MQ=255)
                                                                                                                 cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                       >  2:1030670/1‑81 (MQ=255)
                                                                                                                 cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAgg                                  >  2:768068/1‑150 (MQ=255)
                                                                                                                   gagCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGgtc                                >  2:346879/1‑150 (MQ=255)
                                                                                                                         gTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAg                          >  1:554871/1‑150 (MQ=255)
                                                                                                                         gTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAg                          >  2:1049050/1‑150 (MQ=255)
                                                                                                                              cctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAg                                                                                                 >  2:1264437/1‑74 (MQ=255)
                                                                                                                              cctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAg                                                                                                 <  1:1264437/74‑1 (MQ=255)
                                                                                                                               ctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGc                    >  1:1393171/1‑150 (MQ=255)
                                                                                                                                  ggccACGTGCTCGCCGACCGGGTCGACGTGCACGCGCCCGGAGTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGAt                 <  1:14273/150‑1 (MQ=255)
                                                                                                                                  ggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGAt                 <  2:1394540/150‑1 (MQ=255)
                                                                                                                                  ggccACGTGCTCGCAGACCGGGGCGACGTGCACGAGCCATGCCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGAt                 <  2:533500/150‑1 (MQ=255)
                                                                                                                                   gccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATg                <  2:327102/150‑1 (MQ=255)
                                                                                                                                      aCGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATTGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGACGacattg             <  2:25626/150‑6 (MQ=255)
                                                                                                                                       cGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGTGCGATgtcgt            >  2:1063625/1‑150 (MQ=255)
                                                                                                                                        gTGCACGTCGACCCGGTCGGCGAGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATgtcgtc           >  1:757204/1‑150 (MQ=255)
                                                                                                                                            tCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGACGACATCGCCCAGc       <  1:1204774/150‑1 (MQ=255)
                                                                                                                                                 gACCGGGTCGACGTGCACGAGACATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTACACCAGGTCGTCAATCTGGGCGATGTCGTCCAGCTCCgg  >  2:156453/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGCTCCGG  >  NC_002947/344525‑344819

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