Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,668 C→T 6.3% G313D (GGC→GAC)  mutY ← A/G‑specific adenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,6680CT 6.3% 179.2 / 2.4 63G313D (GGC→GAC) mutYA/G‑specific adenine glycosylase
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base T (2/2);  total (33/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 = 6.95e-01

GGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGC  >  NC_002947/344520‑344814
                                                                                                                                                    |                                                                                                                                                  
ggCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGcc                                                                                                                                                   >  2:611946/1‑150 (MQ=255)
     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:1427287/1‑150 (MQ=255)
                                         gTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCt                                                                                                          >  1:939339/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                                                                                                  <  1:633212/143‑1 (MQ=255)
                                                        gTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                  >  2:633212/1‑143 (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                                                                                    >  2:1225407/1‑150 (MQ=255)
                                                               cagcTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGt                                                                                    <  1:1323050/150‑1 (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                                  <  2:915557/146‑1 (MQ=255)
                                                                                                                     ccAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAgg                                  >  1:915557/1‑146 (MQ=255)
                                                                                                                      cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                  >  2:1030670/1‑81 (MQ=255)
                                                                                                                      cAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTgg                                                                                                  <  1:1030670/81‑1 (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                                                                                            <  1:1264437/74‑1 (MQ=255)
                                                                                                                                   cctcggccACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAg                                                                                            >  2:1264437/1‑74 (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)
                                                                                                                                                    |                                                                                                                                                  
GGCACATCACGGTGCGGGTCATCGGGCGTCTCCTGCGATCAGTTCGTCGGCCGCGCGTTCAAGCAGCTTTTTCACCGGGGCGGCGAGGCCCAGGCGCGGCGGGGTGGCGAGGTTATACCAGAGCCAGTCGGCCTCGGCCACGTGCTCGCCGACCGGGTCGACGTGCACGAGCCATGGCTCGATCGCCAGCTGGAAGTGGCTGAAGGTATGGGTCAGGCCGTCCATGGCGTGGCGGCCAGCCAGGCGCAGGCCGTGCTGGTAGGCCAGGTCGTCGAGCTGGGCGATGTCGTCCAGC  >  NC_002947/344520‑344814

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: ‑