Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,510,295:1 +A 7.2% coding (1479/1725 nt) copA‑I ← copper resistance protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,510,2951.A 7.2% 376.2 / 21.0 97coding (1479/1725 nt)copA‑Icopper resistance protein A
Reads supporting (aligned to +/- strand):  ref base . (53/37);  new base A (2/5);  total (55/42)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGTCGGCTGCCGGGCGGGATATCGACGGTGTGCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGTAGGGGCTGCGCAGGT  >  NC_002947/2510152‑2510443
                                                                                                                                                |                                                                                                                                                    
gTGCGTCGGCTGCCGGGCGGGATATCGACGGTGTGCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCg                                                                                                                                                >  1:682477/1‑150 (MQ=255)
gTGCGTCGGCTGCCGGGCGGGATATCGACGGTGTGCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCg                                                                                                                                                >  1:1555808/1‑150 (MQ=255)
   cGTCGGCTGCCGGGCGGGATCTCGACGGTGTGCTTGCGCCCCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAg                                                                                                                                             <  2:2200019/150‑1 (MQ=255)
            ccGGGCGGGATATCGACGGTGTGCTTGCGAACCAGGAACTGGCCATGCTCGTCTTCCAGGTCACTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTTTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGc                                                                                                                                    >  2:277948/1‑150 (MQ=255)
                 cgggATATCGACGGTGTGCTTGTGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCgg                                                                                                                               >  1:464109/1‑150 (MQ=255)
                    gATATCGACGGTGTGCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTc                                                                                                                            <  2:20310/150‑1 (MQ=255)
                       atCGACGGTGTGCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGc                                                                                                                         <  2:1840522/150‑1 (MQ=255)
                       atCGACGGTGTGCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGc                                                                                                                         >  1:565402/1‑150 (MQ=255)
                               tgtgCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGa                                                                                                                 <  2:361363/150‑1 (MQ=255)
                               tgtgCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGa                                                                                                                 >  1:361363/1‑150 (MQ=255)
                                gtgCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGACGCTGaa                                                                                                                >  1:34235/1‑150 (MQ=255)
                                 tgCTTGCGCACTAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAAc                                                                                                               <  1:1309267/150‑1 (MQ=255)
                                  gCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACt                                                                                                              <  2:565402/150‑1 (MQ=255)
                                     tGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGa                                                                                                           <  2:2228118/150‑1 (MQ=255)
                                        gcACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACAACAGGACAGCTGCAGCGGTTCGGCGTCGCTGAACTTGATGc                                                                                                        >  2:1171149/1‑150 (MQ=38)
                                          aCCAGGAACTGGCCATTCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCg                                                                                                      >  2:599338/1‑150 (MQ=255)
                                          aCCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCg                                                                                                      >  1:801282/1‑150 (MQ=255)
                                            cAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAg                                                                                                                                             >  2:1849340/1‑109 (MQ=255)
                                            cAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAg                                                                                                                                             <  1:1849340/109‑1 (MQ=255)
                                            cAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTc                                                                                                    >  2:1682808/1‑150 (MQ=255)
                                            cAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTc                                                                                                    >  1:1636004/1‑150 (MQ=255)
                                             aGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCa                                                                                                                                                   <  2:2027484/102‑1 (MQ=255)
                                             aGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCa                                                                                                                                                   >  1:2027484/1‑102 (MQ=255)
                                              ggAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGa                                                                                                  >  1:1459902/1‑150 (MQ=255)
                                                      ccATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAg                                                                                          >  2:2002229/1‑150 (MQ=255)
                                                      ccATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAg                                                                                          >  2:1219704/1‑150 (MQ=255)
                                                        aTGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGc                                                                                                                                       >  1:1109933/1‑103 (MQ=255)
                                                        aTGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGc                                                                                                                                       <  2:1109933/103‑1 (MQ=255)
                                                        aTGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGc                                                                                        >  1:1114224/1‑150 (MQ=255)
                                                        aTGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTAATCCTGCTGTCGTTGAAAAAGGTGATGCCCACGCGCTAACCGTAGGTCAGCTGCAGCGGACCCGCGTCGCCGAACTTGATGCCGCCGAACGACCAgcc                                                                                        >  2:418819/1‑148 (MQ=255)
                                                           cTCGTCTTCCAGGTCGCTCCACATGCCGTGCCGGTGGATTGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGaa                                                                                     <  2:1999375/150‑1 (MQ=255)
                                                           cTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGaa                                                                                     <  2:976821/150‑1 (MQ=255)
                                                           cTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGaa                                                                                     >  2:751256/1‑150 (MQ=255)
                                                            tCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGAc                                                                                             <  1:616139/141‑1 (MQ=255)
                                                            tCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGAc                                                                                             >  2:616139/1‑141 (MQ=255)
                                                              gTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGAGCGTGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAgag                                                                                  >  1:1999375/1‑148 (MQ=255)
                                                                cTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCt                                                                                >  2:1445741/1‑150 (MQ=255)
                                                                cTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGACGAAGCGCt                                                                                <  1:1219704/150‑1 (MQ=255)
                                                                 ttCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGc                                                                                                                                 >  1:2177977/1‑100 (MQ=255)
                                                                 ttCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGc                                                                                                                                 <  2:2177977/100‑1 (MQ=255)
                                                                 ttCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGCCGAACGACCAGGCGAAGCGCTc                                                                               >  2:12871/1‑150 (MQ=255)
                                                                 ttCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGAGAAAAACACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTc                                                                               <  1:419097/150‑1 (MQ=38)
                                                                          gCTCCACATGACGCGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTGGACCAGGGTGATGCGGACGCGCTAGCCGTAGGTCAGCCGGTGCCGTTCGGCGACGCTGAACTTGATGCCGTCGAACGACCAGTCGAAGCGCTCCATATGGcc                                                                      <  2:682477/150‑1 (MQ=255)
                                                                           cTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTc                                                                                                                      >  1:2068365/1‑101 (MQ=255)
                                                                           cTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTc                                                                                                                      <  2:2068365/101‑1 (MQ=255)
                                                                            tCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCgg                                                                    >  2:2045951/1‑150 (MQ=255)
                                                                             ccACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGt                                                                   >  1:137935/1‑150 (MQ=255)
                                                                                caTGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAg                                                                <  1:1445741/150‑1 (MQ=255)
                                                                                 aTGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCa                                                                             <  2:2203281/136‑1 (MQ=255)
                                                                                 aTGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCa                                                                             >  1:2203281/1‑136 (MQ=255)
                                                                                 aTGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAgg                                                               >  1:963665/1‑150 (MQ=255)
                                                                                   gCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTg                                                             <  2:801282/150‑1 (MQ=255)
                                                                                         cAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTc                                                       <  2:300077/150‑1 (MQ=255)
                                                                                          aGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGc                                                                                                                         >  1:2348096/1‑83 (MQ=255)
                                                                                          aGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGc                                                                                                                         <  2:2348096/83‑1 (MQ=255)
                                                                                          aGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTAAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCg                                                      >  1:232509/1‑150 (MQ=255)
                                                                                            gTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGAt                                                    >  2:426182/1‑150 (MQ=255)
                                                                                             tGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCGTGGCAGGTGCAGCTCGATg                                                   >  1:1929399/1‑150 (MQ=255)
                                                                                                   gggtgggtCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGc                                                                                        <  2:483736/107‑1 (MQ=255)
                                                                                                   gggtgggtCATCATGGTGTCGTTGACAAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGc                                                                                        >  1:483736/1‑107 (MQ=255)
                                                                                                    ggtgggtCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGAc                                                                                             <  2:348952/101‑1 (MQ=255)
                                                                                                    ggtgggtCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGAc                                                                                             >  1:348952/1‑101 (MQ=255)
                                                                                                    ggtgggtCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCg                                                                     <  1:800207/125‑1 (MQ=255)
                                                                                                    ggtgggtCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCg                                                                     >  2:800207/1‑125 (MQ=255)
                                                                                                    ggtgggtCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCAcgcg                                            >  2:28181/1‑150 (MQ=255)
                                                                                                       ccaacctcatGCTGTCGTTGACCAGGGTGATGCGCCCGCGCTCACCGTAGGTCAGCCGCAGCGGGTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATg                                         <  1:112079/144‑1 (MQ=255)
                                                                                                        ggtCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATgg                                        >  1:897251/1‑150 (MQ=255)
                                                                                                          tcatcaTGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATggcg                                      >  1:1930166/1‑150 (MQ=255)
                                                                                                            atcatGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATggcggg                                    >  2:2293616/1‑150 (MQ=255)
                                                                                                              ccgctgcggctgacctacGGTGAGCGCGTGCGCATCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGt                                                                   <  2:441747‑M1/99‑1 (MQ=255)
                                                                                                              ccgctgcggctgacctacGGTGAGCGCGTGCGCATCACCGTAGGTCAGCCGCAGCGGTTCGGCGTAGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGt                                                                   >  1:441747‑M1/19‑117 (MQ=255)
                                                                                                               atGGTGTCGTTGAACAGGGTGATGCGGACGCGCTCACCGTAGGTCAGAGGCAGCGGTTGGGCGGCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGATCCATATGGCCGGTCAGGTGCAGCCCGATGTAACGCGATggagaacgg                                 >  2:112079/1‑143 (MQ=255)
                                                                                                                   tgtCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCGCGGTGGCACGCGATGGCGGGCGACCGt                             >  1:1963718/1‑150 (MQ=255)
                                                                                                                     tCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTcc                           <  2:897251/150‑1 (MQ=255)
                                                                                                                     tCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTcc                           <  1:58703/150‑1 (MQ=255)
                                                                                                                     tCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTcc                           >  1:1631589/1‑150 (MQ=255)
                                                                                                                     tCGTTGACCAGGGTGATCCGCACGGGGTCACTGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGGCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGGCCGTcc                           <  2:2196989/150‑1 (MQ=255)
                                                                                                                       gTTGACCAGGGTGATGCGCACGCGATCACCGTAGGTCAGCCGCAGAGGTTGGGCGTGGCGGAACTTGATGCCGTAGAACGACCAGGCGAAGCGCTCCATATGGCCGGTAAGGTGCAGCTCGATGTCCCGCGATGGCGGGCGACCGACCgg                         >  1:359090/1‑150 (MQ=255)
                                                                                                                            ctacGGTGAGCGCTTGCGCATCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCgg                    <  2:93284‑M1/146‑1 (MQ=255)
                                                                                                                            ctacGGTGAGCGCGTGCGCATCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCgg                    <  1:590042‑M1/146‑1 (MQ=255)
                                                                                                                            ccAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCgg                    <  1:751256/150‑1 (MQ=255)
                                                                                                                            ccAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCgg                    <  2:963665/150‑1 (MQ=255)
                                                                                                                             tacGGTGAGCGCGTGCGCATCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAg                                                                                          >  2:763385‑M1/4‑79 (MQ=255)
                                                                                                                             tacGGTGAGCGCGTGCGCATCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAg                                                                                          <  1:763385‑M1/76‑1 (MQ=255)
                                                                                                                             tacGGTGAGCGCGTGCGCATCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCggg                   <  2:462075‑M1/147‑1 (MQ=255)
                                                                                                                             cAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCTTTTAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCggg                   >  1:1869218/1‑150 (MQ=255)
                                                                                                                             cAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGa                                 >  1:26598/1‑136 (MQ=255)
                                                                                                                             cAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGa                                 <  2:26598/136‑1 (MQ=255)
                                                                                                                             cAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCggg                   <  2:1636004/150‑1 (MQ=255)
                                                                                                                             cAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCggg                   >  2:146464/1‑150 (MQ=255)
                                                                                                                               gggTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGt                 <  1:599338/150‑1 (MQ=255)
                                                                                                                                GGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGAt                                                                                                          >  1:915787‑M1/1‑60 (MQ=255)
                                                                                                                                GGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGAt                                                                                                          <  2:915787‑M1/60‑1 (MQ=255)
                                                                                                                                  TGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGTAgg              >  2:1631589‑M1/1‑150 (MQ=255)
                                                                                                                                       cgcACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGTAGGGGCTg         >  2:716792/1‑150 (MQ=255)
                                                                                                                                           cgcgCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGTAggg             >  2:203496/1‑142 (MQ=255)
                                                                                                                                           cgcgCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGTAggg             <  1:203496/142‑1 (MQ=255)
                                                                                                                                            gcgcTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTAGCTGAACCTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTGAGTGGGAGTTGGAGGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGTAGGGGCTGCGCAg    >  1:1253357/1‑150 (MQ=255)
                                                                                                                                              gcTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGTAGGGGCTGCGCAGGt  >  1:857646/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTGCGTCGGCTGCCGGGCGGGATATCGACGGTGTGCTTGCGCACCAGGAACTGGCCATGCTCGTCTTCCAGGTCGCTCCACATGCCGTGCAGGTGGATGGGGTGGGTCATCATGGTGTCGTTGACCAGGGTGATGCGCACGCGCTCACCGTAGGTCAGCCGCAGCGGTTCGGCGTCGCTGAACTTGATGCCGTCGAACGACCAGGCGAAGCGCTCCATATGGCCGGTCAGGTGCAGCTCGATGTCACGCGATGGCGGGCGACCGTCCGGGTCGGGGTAGGGGCTGCGCAGGT  >  NC_002947/2510152‑2510443

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: