Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,754,410 A→C 8.3% Q262P (CAG→CCG) ‡ pvdR → ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,754,4100AC 8.3% 183.8 / 4.0 66Q262P (CAG→CCG) ‡pvdRABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (26/34);  new base C (2/4);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCTGATCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGATGGCGCCCTCATGGC  >  NC_002947/4754261‑4754557
                                                                                                                                                      |                                                                                                                                                   
gccgcTGATCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCa                                                                                                                                                     >  1:661314/1‑150 (MQ=255)
gccgcTGATCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCa                                                                                                                                                     >  1:578968/1‑150 (MQ=255)
       aTCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                                              <  2:942015/150‑1 (MQ=255)
                 tCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGc                                                                                                                                    <  1:17532/150‑1 (MQ=255)
                 gtgCCAAGCTGTCCCCCGTGACCGGCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGc                                                                                                                                    <  2:453718/148‑1 (MQ=255)
                     cAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                                              <  2:232524/136‑1 (MQ=255)
                     cAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                                              >  1:232524/1‑136 (MQ=255)
                      aaGCTGTCCCCTATGGCGGTCTGGGCCGAGGTCTCGGACGCCACTACCGAAAACATCAAGACCAGCATGAACGCCTACTTCACCACTCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAg                                                                                                                               <  2:518922/150‑1 (MQ=255)
                       aGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGc                                                                                                                              <  1:1263788/150‑1 (MQ=255)
                            tCCCCCATGGCCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGtgctg                                                                                                                                              <  1:676667/129‑1 (MQ=255)
                            tCCCCCATGGCCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGtgctg                                                                                                                                              >  2:676667/1‑129 (MQ=255)
                            tCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTg                                                                                                                         >  2:868918/1‑150 (MQ=255)
                             cccccATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTgg                                                                                                                        >  2:231602/1‑150 (MQ=255)
                             cccccATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTgg                                                                                                                        >  2:808561/1‑150 (MQ=255)
                               cccATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGTGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgac                                                                                                                      >  2:319227/1‑150 (MQ=255)
                                ccATGACCGTCTGGGCCCAGGTCTCGGAGGCGGATAGCTGCAAGATCCATCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgacc                                                                                                                     >  2:1019875/1‑150 (MQ=255)
                                ccATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgacc                                                                                                                     >  2:29564/1‑150 (MQ=255)
                                 cATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgacca                                                                                                                    >  1:1409817/1‑150 (MQ=255)
                                   tGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgaccaga                                                                                                                  <  2:152450/150‑1 (MQ=255)
                                      ccGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGgaccagacca                                                                                                               <  1:667014/150‑1 (MQ=255)
                                          cTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTCCTTCCCCACCCTGGCCTTCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagcca                                                                                                           <  1:882394/150‑1 (MQ=255)
                                           tGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagccag                                                                                                          <  1:231602/150‑1 (MQ=255)
                                                 cAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGggcggc                                                                                                    <  1:1061690/150‑1 (MQ=255)
                                                        cGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGcagcc                                                                                             <  1:808561/150‑1 (MQ=255)
                                                        cGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTACCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGcagcc                                                                                             <  2:578968/150‑1 (MQ=255)
                                                             gCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGcagcccagcc                                                                                        >  2:353626/1‑150 (MQ=255)
                                                                gATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCcagccagc                                                                                     >  2:1062911/1‑150 (MQ=255)
                                                                gATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCcagccagc                                                                                     <  1:29564/150‑1 (MQ=255)
                                                                    tCGGCAAGATCCAGCCTGGTATGACCGCCTACTTCACCCCCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCa                                                                                 <  2:939267/150‑1 (MQ=255)
                                                                      ggCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                                              >  1:805581/1‑87 (MQ=255)
                                                                      ggCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCTGGCGCACGGtgctg                                                                                                                                              <  2:805581/87‑1 (MQ=255)
                                                                          aGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCAccgccc                                                                           >  2:1417948/1‑149 (MQ=255)
                                                                                  ccTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCAcc                                                                   >  2:319809/1‑150 (MQ=255)
                                                                                      gCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagccag                                                                                                          <  2:8867/107‑1 (MQ=255)
                                                                                      gCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagccag                                                                                                          >  1:8867/1‑107 (MQ=255)
                                                                                        atccccGCATACTTCAACACCCCGGCCGGCGGCAACCCCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGc                                                             <  2:1084069/146‑1 (MQ=255)
                                                                                           aCCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGCGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCaa                                                          <  2:556936/150‑1 (MQ=255)
                                                                                            ccGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCAGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAg                                                         >  1:394392/1‑150 (MQ=255)
                                                                                             cGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGt                                                        >  1:780314/1‑150 (MQ=255)
                                                                                                cTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCAGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGt                                                     <  2:394392/150‑1 (MQ=255)
                                                                                                 tACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGCGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTg                                                    <  1:543389/150‑1 (MQ=255)
                                                                                                  aCTTCACCACCCTGGCCGGCGGCAAAACGCCCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGc                                                   >  1:1327847/1‑150 (MQ=255)
                                                                                                   cTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCa                                                  <  1:1417948/150‑1 (MQ=255)
                                                                                                    ttCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCaa                                                                                                                                >  2:93161/1‑71 (MQ=255)
                                                                                                    ttCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCaa                                                                                                                                <  1:93161/71‑1 (MQ=255)
                                                                                                    ttCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGTCCAACAAGGGACTCCTTCGTCTCCTCCCGCAGCCTAGTCGTGCac                                                 >  2:156027/1‑149 (MQ=255)
                                                                                                    ttCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCaa                                                 >  2:345199/1‑150 (MQ=255)
                                                                                                         caccCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATacac                                            <  2:780314/150‑1 (MQ=255)
                                                                                                            ccTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagcca                                                                                                           >  1:1060661/1‑84 (MQ=255)
                                                                                                            ccTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagcca                                                                                                           <  2:1060661/84‑1 (MQ=255)
                                                                                                            ccTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagc                                                                                                             <  1:1005518/82‑1 (MQ=255)
                                                                                                            ccTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGAccagc                                                                                                             >  2:1005518/1‑82 (MQ=255)
                                                                                                              tGGCCGGCGGCAAGCGCCGCTGGACCAGCACGGTGCGCACGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCACCGACACCGCCGGACCCACCGGTAGCCAGGTCGCTCTATACTCCGtgc                                       <  2:1156099/150‑1 (MQ=255)
                                                                                                                  cggcggCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTgctgct                                   <  1:868918/150‑1 (MQ=255)
                                                                                                                   ggcggcAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACCCCGAGctgctg                                  <  2:354577/150‑1 (MQ=255)
                                                                                                                      ggcAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGAc                               <  2:1409817/150‑1 (MQ=255)
                                                                                                                        cAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGt                             <  1:266608/150‑1 (MQ=255)
                                                                                                                        cAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGt                             <  1:353626/150‑1 (MQ=255)
                                                                                                                            gcgcCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGa                          <  1:1062911/150‑1 (MQ=255)
                                                                                                                            gcgcCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACAAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATCCACCGTGCTGCTGGACGTGGa                          >  1:350569/1‑150 (MQ=255)
                                                                                                                                     gACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGTCGGCAGCGGCCCAGCCAGCGCAACCGCCGACGCCACCGGCAGCGAAGTCCTCCAATACGCCGGGCTGCTCGACGTGGACGAcccgaa                 >  2:140784/1‑147 (MQ=255)
                                                                                                                                       ccAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGgcgc                                                                       >  1:649410/1‑94 (MQ=255)
                                                                                                                                       ccAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGgcgc                                                                       <  2:649410/94‑1 (MQ=255)
                                                                                                                                                 gcgcCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACAAGACACGCCAGGGGTTCGTCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGATGGCGCCCTCAt     <  1:319809/150‑1 (MQ=255)
                                                                                                                                                   gcCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGATGGCGCCCTCATgg   <  2:946812/150‑1 (MQ=255)
                                                                                                                                                    ccAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGAt                >  1:458936/1‑136 (MQ=255)
                                                                                                                                                    ccAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGAt                <  2:458936/136‑1 (MQ=255)
                                                                                                                                                    ccAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGATGGCGCCCTCATGGc  >  1:121949/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCCGCTGATCCTGCGCATCGCCAAGCTGTCCCCCATGACCGTCTGGGCCCAGGTCTCGGAGGCCGATATCGGCAAGATCCAGCCTGGCATGACCGCCTACTTCACCACCCTGGCCGGCGGCAAGCGCCGCTGGACCAGCACCGTGCGCCAGGTGCTGCCCATTCCGCCCAAGCCACTGGACCAGACCAGCCAGGGCGGCGGCAGCCCAGCCAGCGCCACCGCCGGCGCCACCGGCAGCCAAGTCGTGCAATACACCGTGCTGCTGGACGTGGACAACCCCGATGGCGCCCTCATGGC  >  NC_002947/4754261‑4754557

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