Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,777,050 C→G 7.5% W166C (TGG→TGC PP_2429 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,777,0500CG 7.5% 295.7 / 11.1 93W166C (TGG→TGCPP_2429membrane protein
Reads supporting (aligned to +/- strand):  ref base C (43/43);  new base G (3/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTT  >  NC_002947/2776904‑2777199
                                                                                                                                                  |                                                                                                                                                     
gcgAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCAc                                                                                                                                                    >  2:2254827/1‑150 (MQ=255)
   aGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCgg                                                                                                                                                 <  1:1627321/150‑1 (MQ=255)
      cGTTGACACGTTGAGCCGAGGCCCCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGGCTGCAGGCACCggggc                                                                                                                                              <  1:721043/150‑3 (MQ=255)
        ttGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCg                                                                                                                                                  <  1:1518555/144‑1 (MQ=255)
        ttGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCg                                                                                                                                                  >  2:1518555/1‑144 (MQ=255)
         tGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTc                                                                                                                                           <  2:2251954/150‑1 (MQ=255)
               gTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCAcc                                                                                                                                     >  1:164601/1‑150 (MQ=255)
                      cGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAaca                                                                                                                              <  2:693358/150‑1 (MQ=255)
                               tCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAgg                                                                                                                     <  2:2053630/150‑1 (MQ=255)
                                cAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAggc                                                                                                                    <  2:2101823/150‑1 (MQ=255)
                                        ttGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCACCGGCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGc                                                                                                            >  2:2133058/1‑150 (MQ=255)
                                            cGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCgg                                                                                                        <  1:2254827/150‑1 (MQ=255)
                                             ggTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGc                                                                                                       >  2:51550/1‑150 (MQ=255)
                                              gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAg                                                                                                                      >  2:232414/1‑134 (MQ=255)
                                              gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAg                                                                                                                      <  1:232414/134‑1 (MQ=255)
                                               tCAAGTAGAGGCTCAGGCCAATACTGCATCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCag                                                                                                     >  1:1771847/1‑150 (MQ=255)
                                                      gagGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATc                                                                                              >  2:2224560/1‑150 (MQ=255)
                                                         gCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGc                                                                                           >  1:1782029/1‑150 (MQ=255)
                                                         gCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGc                                                                                           >  2:2013068/1‑150 (MQ=255)
                                                          cTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGc                                                                                           <  2:1306196/149‑1 (MQ=255)
                                                          cTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGc                                                                                           >  1:1306196/1‑149 (MQ=255)
                                                          cTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCa                                                                                          >  2:1337448/1‑150 (MQ=255)
                                                               gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCaa                                                                                                                                >  2:1347503/1‑107 (MQ=255)
                                                               gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCaa                                                                                                                                <  1:1347503/107‑1 (MQ=255)
                                                               gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag                                                                                                   >  1:212219/1‑136 (MQ=255)
                                                               gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag                                                                                                   <  2:212219/136‑1 (MQ=255)
                                                               gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAac                                                                                     >  2:1688534/1‑150 (MQ=255)
                                                                 cAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGGCTGCAGGCACCggg                                                                                                                                                <  1:1123761/89‑1 (MQ=255)
                                                                 cAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGGCTGCAGGCACCggg                                                                                                                                                >  2:1123761/1‑89 (MQ=255)
                                                                    tACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCAccccc                                                                                                                                  <  1:1274689/100‑1 (MQ=255)
                                                                    tACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCAccccc                                                                                                                                  >  2:1274689/1‑100 (MQ=255)
                                                                       tGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAaca                                                                                    >  2:1067311/1‑143 (MQ=255)
                                                                       tGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAaca                                                                                    <  1:1067311/143‑1 (MQ=255)
                                                                       tGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgcgc                                                                             >  1:430215/1‑150 (MQ=255)
                                                                       tGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgcgc                                                                             >  1:1843939/1‑150 (MQ=255)
                                                                         cAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCa                                                                           <  1:2400670/150‑1 (MQ=255)
                                                                               aTGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGcc                                                                     >  1:1991440/1‑150 (MQ=255)
                                                                                 gCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCcac                                                                   >  2:1096442/1‑150 (MQ=255)
                                                                                     atatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACtc                                                               >  2:57239/1‑149 (MQ=255)
                                                                                     atatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTg                                                               <  2:1782029/150‑1 (MQ=255)
                                                                                     atatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTg                                                               <  2:2085162/150‑1 (MQ=255)
                                                                                          aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                                          >  2:1926251/1‑150 (MQ=255)
                                                                                          aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                                          >  2:2052403/1‑150 (MQ=255)
                                                                                           acaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacg                                                         >  2:1676279/1‑150 (MQ=255)
                                                                                            caaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCCCCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacga                                                        <  2:1771847/150‑1 (MQ=255)
                                                                                            caaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacga                                                        <  1:1096442/150‑1 (MQ=255)
                                                                                             aaaaCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag                                                                                                   >  1:1310924/1‑106 (MQ=255)
                                                                                             aaaaCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag                                                                                                   <  2:1310924/106‑1 (MQ=255)
                                                                                                  ccccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcgc                                                  <  2:430215/150‑1 (MQ=255)
                                                                                                  ccccGCCAACACTTGGTTGAACACCCGCAACCACCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCCCGGCAGAGCGATCAGCAAGAACACCAGCCCCAGCTGCCACACTGGCAACGACGgcgc                                                  >  1:5880/1‑150 (MQ=255)
                                                                                                    ccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCgg                                                >  2:1606448/1‑150 (MQ=255)
                                                                                                     cGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGcgcgtc                                               <  1:1676279/150‑3 (MQ=255)
                                                                                                     cGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGc                                               <  1:2133058/150‑1 (MQ=255)
                                                                                                        cAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa                                                                                         >  1:694643/1‑105 (MQ=255)
                                                                                                        cAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa                                                                                         <  2:694643/105‑1 (MQ=255)
                                                                                                            acTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCacac                                                                 >  2:2204750/1‑125 (MQ=255)
                                                                                                            acTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCacac                                                                 <  1:2204750/125‑1 (MQ=255)
                                                                                                             cTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAccc                                       <  2:1991440/150‑1 (MQ=255)
                                                                                                              ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACCCCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAcccc                                      <  2:161139/150‑1 (MQ=255)
                                                                                                                  ttGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGacc                                  >  2:2207716/1‑150 (MQ=255)
                                                                                                                    gAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGa                                                                                                <  1:2349783/86‑1 (MQ=255)
                                                                                                                    gAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGa                                                                                                >  2:2349783/1‑86 (MQ=255)
                                                                                                                     aaCACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAcccc                                      <  2:835046/143‑1 (MQ=255)
                                                                                                                     aaCACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAcccc                                      >  1:835046/1‑143 (MQ=255)
                                                                                                                     aaCACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGaccacc                               >  1:2282007/1‑150 (MQ=255)
                                                                                                                       cacCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGc                             >  1:971884/1‑150 (MQ=255)
                                                                                                                          ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCa                                                                                             <  1:693330/83‑1 (MQ=255)
                                                                                                                          ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCa                                                                                             >  2:693330/1‑83 (MQ=255)
                                                                                                                          ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCaaa                                            <  2:159382/132‑1 (MQ=255)
                                                                                                                          ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCaaa                                            >  1:159382/1‑132 (MQ=255)
                                                                                                                          ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAccc                                       <  1:139252/137‑1 (MQ=255)
                                                                                                                          ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAccc                                       >  2:139252/1‑137 (MQ=255)
                                                                                                                            gCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACgg                        <  1:2207716/150‑1 (MQ=255)
                                                                                                                               aCCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCa                     <  2:2282007/150‑1 (MQ=255)
                                                                                                                                ccGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCAt                    <  1:1606448/150‑1 (MQ=255)
                                                                                                                                  gCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag                                                                                                   <  1:1721423/69‑1 (MQ=255)
                                                                                                                                  gCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag                                                                                                   >  2:1721423/1‑69 (MQ=255)
                                                                                                                                  gCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCa                  >  2:1661314/1‑150 (MQ=255)
                                                                                                                                    cccGGTGGCTGCAGGCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCaa                                                           <  2:1219988/107‑1 (MQ=255)
                                                                                                                                    cccGGTGGCTGCAGGCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCaa                                                           >  1:1219988/1‑107 (MQ=255)
                                                                                                                                    cccGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCg                                                                                                 <  2:1788075/69‑1 (MQ=255)
                                                                                                                                    cccGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCg                                                                                                 >  1:1788075/1‑69 (MQ=255)
                                                                                                                                     ccGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAcc               <  1:1926251/150‑1 (MQ=255)
                                                                                                                                      cGGTGGCTGCAGGCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGc                                                                                                       >  2:1368046/1‑61 (MQ=38)
                                                                                                                                      cGGTGGCTGCAGGCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGc                                                                                                       <  1:1368046/61‑1 (MQ=38)
                                                                                                                                             tGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACt                                                                >  2:1387358/1‑93 (MQ=255)
                                                                                                                                             tGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACt                                                                <  1:1387358/93‑1 (MQ=255)
                                                                                                                                             tGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTg       <  2:2038146/150‑1 (MQ=255)
                                                                                                                                             tGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTg       >  1:1820774/1‑150 (MQ=255)
                                                                                                                                               cAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCCCGGCCATAAACCCGACCTTggg     <  2:1843939/150‑1 (MQ=255)
                                                                                                                                               cAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTggg     <  1:1337448/150‑1 (MQ=255)
                                                                                                                                               cAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTggg     <  2:971884/150‑1 (MQ=255)
                                                                                                                                                aGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTgggg    <  1:1661314/150‑1 (MQ=255)
                                                                                                                                                  ccACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCACACCTTGGGGtt  >  2:1126627/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGTT  >  NC_002947/2776904‑2777199

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