Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,610,313 Δ1 bp 6.2% intergenic (‑38/+18) PP_3183 ← / ← mcoA SCO1/SenC family protein/Mn(II) copper oxidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,610,3130G. 6.2% 288.3 / 14.3 97intergenic (‑38/+18)PP_3183/mcoASCO1/SenC family protein/Mn(II) copper oxidase A
Reads supporting (aligned to +/- strand):  ref base G (48/39);  new base . (3/3);  total (54/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

CCAGGGAGAGGCCGCTGCTTCGGCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGCCCCCGTAGTCGTGGCA  >  NC_002947/3610165‑3610460
                                                                                                                                                         |                                                                                                                                                   
ccAGGGAGAGGCCGCTGCTTCGGCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATgg                                                                                                                                                    >  1:1131809/1‑150 (MQ=255)
   gggAGAGGCCGCTGCTTCGGCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATggcgg                                                                                                                                                 >  2:1888713/1‑150 (MQ=255)
    ggAGAGGCCGCTGCTTCGGCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATggcggc                                                                                                                                                >  1:1032309/1‑150 (MQ=255)
       gagGCCGCTGCTTCGGCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCCATGCGGcccc                                                                                                                                             <  1:2208741/150‑1 (MQ=255)
          gCCGCTGCTTCGGCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACACC‑‑GAATCGTTCTGAGCAGCCGGTCATGATTCGGTACTCGTAGGTTGGCTGATGCGTCGCGCATGGCGGCCCCGCt                                                                                                                                          >  1:44278/1‑150 (MQ=255)
             gctgctTCGGCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTg                                                                                                                                       >  1:1682198/1‑150 (MQ=255)
                      gCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACtt                                                                                                                              >  1:2278523/1‑150 (MQ=255)
                        cACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGt                                                                                                                            >  1:1817341/1‑150 (MQ=255)
                          ccGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAt                                                                                                                          <  2:1998132/150‑1 (MQ=255)
                           cgttcgttGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAtg                                                                                                                         <  2:1342136/150‑1 (MQ=255)
                            gttcgttGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAtgt                                                                                                                        <  1:1190332/150‑1 (MQ=255)
                            gttcgttGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAtgt                                                                                                                        >  1:375245/1‑150 (MQ=255)
                            gttcgttGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAtgt                                                                                                                        <  1:839182/150‑1 (MQ=255)
                                 ttGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCttt                                                                                                                                        <  1:1895469/129‑1 (MQ=255)
                                 ttGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCttt                                                                                                                                        >  2:1895469/1‑129 (MQ=255)
                                   ggggTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGgc                                                                                                                 <  2:959947/150‑1 (MQ=255)
                                    gggTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCCGAGATTGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGCCTGgcg                                                                                                                >  1:2331648/1‑150 (MQ=255)
                                      gTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCt                                                                                                              <  2:1682198/150‑1 (MQ=255)
                                      gTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCt                                                                                                              >  1:2204883/1‑150 (MQ=255)
                                       ttGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGTAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAt                                                                                                                          >  2:666154/1‑137 (MQ=255)
                                       ttGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGTAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAt                                                                                                                          <  1:666154/137‑1 (MQ=255)
                                           aaCGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCAc                                                                                                         >  2:2251089/1‑150 (MQ=255)
                                                ccAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCCGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCt                                                                                                    >  1:1953551/1‑150 (MQ=255)
                                                 cAGAACGGAATGCTGGCGGCGAACAACGCCAGTAGCAGCAGGAATCG‑‑NCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCACGCGACACACTTGTATGTGACTGGCGCAGACAccagcag                                                                                                   >  1:738958/1‑145 (MQ=255)
                                                 cAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCCATGCGGccccc                                                                                                                                            <  1:459486/109‑2 (MQ=255)
                                                 cAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCCATGCGGccccc                                                                                                                                            >  2:459486/1‑108 (MQ=255)
                                                    aaCGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGcc                                                                                                <  1:1888713/150‑1 (MQ=255)
                                                      cGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACtt                                                                                                                              <  1:2144731/118‑1 (MQ=255)
                                                      cGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACtt                                                                                                                              >  2:2144731/1‑118 (MQ=255)
                                                          aTGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGt                                                                                          <  2:302944/150‑1 (MQ=255)
                                                           tGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACtt                                                                                                                              <  2:669437/113‑1 (MQ=255)
                                                           tGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACtt                                                                                                                              >  1:669437/1‑113 (MQ=255)
                                                           tGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAt                                                                                                                          <  1:318601/117‑1 (MQ=255)
                                                           tGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAt                                                                                                                          >  2:318601/1‑117 (MQ=255)
                                                             cTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGa                                                                                       <  2:1424947/150‑1 (MQ=255)
                                                              tGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTa                                                                                                                           >  2:1311612/1‑113 (MQ=255)
                                                              tGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTa                                                                                                                           <  1:1311612/113‑1 (MQ=255)
                                                              tGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGAt                                                                                      >  1:2125606/1‑150 (MQ=255)
                                                              tGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGAt                                                                                      >  2:698077/1‑150 (MQ=255)
                                                                  ggcgAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCa                                                                                  >  2:1160168/1‑150 (MQ=255)
                                                                     gAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCG                                                                               >  1:1319914‑M1/1‑150 (MQ=255)
                                                                      aacaacGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGT                                                                              >  1:364862‑M1/1‑150 (MQ=255)
                                                                         aacGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGAC                                                                           >  1:29446‑M1/1‑150 (MQ=255)
                                                                         aacGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGAC                                                                           >  1:52870‑M1/1‑150 (MQ=255)
                                                                         aacGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCAGCCTCGCCTTGAACTCAGAGGTAGGTGGCGGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGAC                                                                           >  1:1250820‑M1/1‑150 (MQ=255)
                                                                           cGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCg                                                                         >  1:1425366/1‑150 (MQ=255)
                                                                            gCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGt                                                                        >  1:881411/1‑150 (MQ=255)
                                                                              cAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCg                                                                      <  2:375245/150‑1 (MQ=255)
                                                                                ttACCGGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCt                                                                    >  2:1776719/2‑150 (MQ=255)
                                                                                gTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCCGCGTCGATTTCACCGTGACCGTCGCt                                                                    <  2:44278/150‑1 (MQ=255)
                                                                                gTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCt                                                                    >  2:2361594/1‑150 (MQ=255)
                                                                                 tACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTg                                                                   >  2:1068635/1‑150 (MQ=255)
                                                                                      gcagGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTg                                                              >  2:2183054/1‑150 (MQ=255)
                                                                                      gcagGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTg                                                              >  1:1381917/1‑150 (MQ=255)
                                                                                       cagGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTgg                                                             <  1:1705487/150‑1 (MQ=255)
                                                                                           aaTCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGc                                                                                                                                           >  1:1130822/1‑68 (MQ=255)
                                                                                           aaTCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGc                                                                                                                                           <  2:1130822/68‑1 (MQ=255)
                                                                                           aaTCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCgg                                                         <  1:1776719/150‑1 (MQ=255)
                                                                                           aaTCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCgg                                                         <  2:2125606/150‑1 (MQ=255)
                                                                                            aTCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTAACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCggg                                                        <  2:408973/150‑1 (MQ=255)
                                                                                             tCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTg                                                                                                                             <  1:1180343/80‑1 (MQ=255)
                                                                                             tCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTg                                                                                                                             >  2:1180343/1‑80 (MQ=255)
                                                                                                    cATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGTGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCaaa                                                   >  1:554765/1‑150 (MQ=255)
                                                                                                     aTGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAg                                                  >  1:1343061/1‑150 (MQ=255)
                                                                                                      tGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTg                                                                                                                                       >  1:1758606/1‑64 (MQ=255)
                                                                                                      tGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTg                                                                                                                                       <  2:1758606/64‑1 (MQ=255)
                                                                                                        aGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCg                                               <  1:2361594/150‑1 (MQ=255)
                                                                                                           cGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCggtg                                            <  2:1032309/150‑1 (MQ=255)
                                                                                                            gggCATGATGCGGTCCTCGTGGGTTGGCTGGTGGGGGCACCGTATGGGCCCTGACCTGCACTCACTGGTAGGTGACGGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGTGCCAAAGCCggtgt                                           <  2:1131809/150‑2 (MQ=255)
                                                                                                              gtattatgCGGTACTAGTGGGTTGGCTGATGGGGGCCTCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCggtggtg                                         <  2:1343061/145‑1 (MQ=255)
                                                                                                              gCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTg                                         <  2:2405175‑M2/150‑2 (MQ=255)
                                                                                                               cATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTgg                                        <  1:1248521‑M2/150‑3 (MQ=255)
                                                                                                                   atgAGGTCTTCGTGGGTTGTCGGTTGGGGGCCGTATGGCGGCCCCGCTTTGCACTCGCTTGTATGTGACTGGCGCTGTCACTGTCTGGCCCAGCGTCGATTTCAGCGTGACCGTCGCTGGCGTTTCCGGGCCAAAGCCGGTGGTGGTCAc                                    <  2:989120/150‑1 (MQ=255)
                                                                                                                    tgCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCAcc                                   >  2:514272/1‑150 (MQ=255)
                                                                                                                       ggTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCt                                >  2:278221/1‑150 (MQ=255)
                                                                                                                          ccTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTAcc                             <  2:1319914/150‑1 (MQ=255)
                                                                                                                           cTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTAGGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTAccc                            >  2:1782826/1‑150 (MQ=255)
                                                                                                                            tCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGAt                                                                                      <  1:1658217/91‑1 (MQ=255)
                                                                                                                            tCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGAt                                                                                      >  2:1658217/1‑91 (MQ=255)
                                                                                                                                  gTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACg                     >  2:670841/1‑150 (MQ=255)
                                                                                                                                  gTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACg                     <  2:1521009/150‑1 (MQ=255)
                                                                                                                                        tGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGcccc               >  1:707574/1‑150 (MQ=255)
                                                                                                                                        tGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGcccc               <  1:278221/150‑1 (MQ=255)
                                                                                                                                        aaagcggggCCGCCATGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGcccc               >  2:2404190/6‑150 (MQ=255)
                                                                                                                                        aaagcggggCCGCCATGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGcccc               >  1:383821/6‑150 (MQ=255)
                                                                                                                                          agcggggccgCCATGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTg                                         <  1:1809191‑M2/119‑2 (MQ=255)
                                                                                                                                          agcggggccgCCATGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTg                                         >  2:1809191‑M2/4‑121 (MQ=255)
                                                                                                                                            gggggCCGCATGGGGACCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCaccaccggctttggcccggccacgccagcgacggtcacg                                   <  1:127205‑M1/126‑40 (MQ=255)
                                                                                                                                            gggggCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCaccaccggctttggcccggccacgccagcgacggtcacg                                   >  2:127205‑M1/1‑87 (MQ=255)
                                                                                                                                             ggggCCGCCATGCGGCCCCGCTTTGCACTCAATTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGcc                                                           <  2:1262488/101‑1 (MQ=255)
                                                                                                                                             ggggCCGCCATGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTggtcagc                                   >  2:405931‑M2/1‑118 (MQ=255)
                                                                                                                                             ggggCCGCCATGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCAcc                                   >  1:405931/1‑125 (MQ=255)
                                                                                                                                             ggggCCGCCATGCGGCCCCGCTTTGCACTCAATTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGcc                                                           >  1:1262488/1‑101 (MQ=255)
                                                                                                                                                gCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCaccaccggctttggcccggccacgccagcgacggtcacg                                   >  2:1675612‑M1/1‑83 (MQ=255)
                                                                                                                                                gCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCaccaccggctttggcccggccacgccagcgacggtcacg                                   <  1:1675612‑M1/122‑40 (MQ=255)
                                                                                                                                                gCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGCCCCCGTAGTCg       >  1:667238/1‑150 (MQ=255)
                                                                                                                                                gCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGCCCCCGTAGTCg       >  1:860248/1‑150 (MQ=255)
                                                                                                                                                   gCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGCCCCCGTAGTCGTgg    <  2:1817341/150‑1 (MQ=255)
                                                                                                                                                       gtGGCTGCCCTGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCTCCGTCTGGCCCAGCGTCGATGTCACCGTGACCGCCTCTGGCGTGGACGGGCCAAAGCCGGTTGTGGGCACCGCTACCCTCCAACGTGCCCCCGTAGGCGTGGCa  <  2:738958/149‑1 (MQ=255)
                                                                                                                                                       aTGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGCCCCCGTAGTCGTGGCa  >  1:1193455/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
CCAGGGAGAGGCCGCTGCTTCGGCCACCGTTCGTTGGGGTTGCAACGGCCAGAACGGAATGCTGGCGGCGAACAACGCCAGTACCAGCAGGAATCG‑‑CCATGAGCCGGGCATGATGCGGTCCTCGTGGGTTGGCTGATGGGGGCCGCATGGCGGCCCCGCTTTGCACTCACTTGTATGTGACTGGCGCTGTCACCGTCTGGCCCAGCGTCGATTTCACCGTGACCGTCGCTGGCGTGGCCGGGCCAAAGCCGGTGGTGGTCACCGCTACCCGCCAACGGGCCCCCGTAGTCGTGGCA  >  NC_002947/3610165‑3610460

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