Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,207,691 G→A 6.0% A439V (GCC→GTC)  PP_1058 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,207,6910GA 6.0% 276.6 / 4.6 94A439V (GCC→GTC) PP_1058hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base A (2/4);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGCGCGACCCAGTCGCGGTAC  >  NC_002947/1207543‑1207840
                                                                                                                                                     |                                                                                                                                                     
ggcggTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTGACTCCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGc                                                                                                                                                      <  2:63768/150‑1 (MQ=255)
      gCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCaa                                                                                                                                                <  2:303024/150‑1 (MQ=255)
      gCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGCCGACCTGCaa                                                                                                                                                <  2:1395264/150‑1 (MQ=255)
       cGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTCAGCGTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGCCGACCTGCaac                                                                                                                                               <  1:1842643/150‑2 (MQ=255)
            gAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACt                                                                                                                                          >  2:1269757/1‑150 (MQ=255)
            gAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACt                                                                                                                                          >  1:1754032/1‑150 (MQ=255)
             aaCAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtg                                                                                                                                         >  2:2431534/1‑150 (MQ=255)
              aCAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtgg                                                                                                                                        <  1:784077/150‑1 (MQ=255)
               cAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtggt                                                                                                                                       <  1:1275632/150‑1 (MQ=255)
                 aGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtggtgg                                                                                                                                     <  2:2233617/150‑1 (MQ=255)
                  ggTGAAACTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCCACTGTATGTCAGTCGTGTGCTCTTTTGCTGCAACCCCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGTCCTGCAAGGTACTGGTGGc                                                                                                                                    <  1:280508/150‑1 (MQ=255)
                   gTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGc                                                                                                                                    >  2:2364784/1‑149 (MQ=255)
                   gTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGc                                                                                                                                    <  1:2364784/149‑1 (MQ=255)
                   gTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCt                                                                                                                                   <  1:2348734/150‑1 (MQ=255)
                       aGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACa                                                                                                                               <  2:1899874/150‑1 (MQ=255)
                       aGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACa                                                                                                                               <  2:1849033/150‑1 (MQ=255)
                          tGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGt                                                                                                                            <  1:1266184/150‑1 (MQ=255)
                           gTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTc                                                                                                                           <  1:2328993/150‑1 (MQ=255)
                              gCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAg                                                                                                                        >  2:1196512/1‑150 (MQ=255)
                               cTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAgg                                                                                                                       <  2:1413102/150‑1 (MQ=255)
                               cTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAgg                                                                                                                       <  1:1289031/150‑1 (MQ=255)
                                 cACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCg                                                                                                                     <  1:708011/150‑1 (MQ=255)
                                    cTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTc                                                                                                                  >  1:2357778/1‑150 (MQ=255)
                                    cTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTATCCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTc                                                                                                                  <  1:746190/150‑1 (MQ=255)
                                      ggTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGa                                                                                                                <  2:407305/150‑1 (MQ=255)
                                          ggtggtCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtgg                                                                                                                                        <  2:2427308/122‑1 (MQ=255)
                                          ggtggtCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtgg                                                                                                                                        >  1:2427308/1‑122 (MQ=255)
                                          ggtggtCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTc                                                                                                            >  1:2026633/1‑150 (MQ=255)
                                                 ttGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGc                                                                                                     <  1:2431534/150‑1 (MQ=255)
                                                 ttGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGc                                                                                                     >  1:2447376/1‑150 (MQ=255)
                                                 ttGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGGTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGc                                                                                                     >  2:185957/1‑150 (MQ=255)
                                                    gCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAg                                                                                                                        <  2:1527489/128‑1 (MQ=255)
                                                    gCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAg                                                                                                                        >  1:1527489/1‑128 (MQ=255)
                                                    gCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCggg                                                                                                  >  1:188949/1‑150 (MQ=255)
                                                     cAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGTAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGc                                                                                                 <  1:497391/150‑1 (MQ=255)
                                                       gTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGa                                                                                               <  1:593255/150‑1 (MQ=255)
                                                        tGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGCCGACCTGCaac                                                                                                                                               >  1:786025/1‑100 (MQ=255)
                                                        tGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGCCGACCTGCaac                                                                                                                                               <  2:786025/101‑2 (MQ=255)
                                                          agcagcCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTg                                                                                            >  1:2417518/1‑150 (MQ=255)
                                                           gcagcCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGa                                                                                           <  1:1735318/150‑1 (MQ=255)
                                                            cagcCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAg                                                                                          >  1:316232/1‑150 (MQ=255)
                                                                  cTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGc                                                                                    <  1:362691/150‑1 (MQ=255)
                                                                  cTCGCCGATCAGGCCGATCTGTGCGGAAAATTCAGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGc                                                                                    <  2:188949/150‑1 (MQ=255)
                                                                     gCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGa                                                                                                                >  2:1151511/1‑119 (MQ=255)
                                                                     gCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGa                                                                                                                <  1:1151511/119‑1 (MQ=255)
                                                                        gATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAg                                                                              >  1:139230/1‑150 (MQ=255)
                                                                        gATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGACTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAg                                                                              >  1:1859657/1‑150 (MQ=255)
                                                                         aTCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAgg                                                                             <  2:540033/150‑1 (MQ=255)
                                                                          tCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGc                                                                            >  2:2270914/1‑150 (MQ=255)
                                                                                      tgtgCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGcc                                                                 >  2:796424/1‑150 (MQ=255)
                                                                                      tgtgCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGcc                                                                 >  2:1934699/1‑150 (MQ=255)
                                                                                      tgtgCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGcc                                                                 >  1:2427789/1‑150 (MQ=255)
                                                                                          cGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCt                                                                                                                                   <  2:185252/80‑1 (MQ=255)
                                                                                          cGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCt                                                                                                                                   >  1:185252/1‑80 (MQ=255)
                                                                                            gCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTg                                                           >  2:1475174/1‑150 (MQ=255)
                                                                                              aaaCTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCGGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGAt                                                         >  2:1468163/1‑150 (MQ=255)
                                                                                               aaCTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGAt                                                                                              <  2:2207150/112‑1 (MQ=255)
                                                                                               aaCTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGAt                                                                                              >  1:2207150/1‑112 (MQ=255)
                                                                                                  tCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGccc                                                     >  1:964045/1‑150 (MQ=255)
                                                                                                    cGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCtt                                                   <  1:2222265/150‑1 (MQ=255)
                                                                                                     gggTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGTGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCttg                                                  >  2:1564279/1‑150 (MQ=255)
                                                                                                            cAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTaa                                           <  2:316232/150‑1 (MQ=255)
                                                                                                            cAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTaa                                           >  2:1060400/1‑150 (MQ=255)
                                                                                                             aaGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGAtt                                                                                             <  1:2272273/99‑1 (MQ=255)
                                                                                                             aaGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGAtt                                                                                             >  2:2272273/1‑99 (MQ=255)
                                                                                                             aaGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAAt                                          >  2:862091/1‑150 (MQ=255)
                                                                                                             aaGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGATTGCGCGCCATGCTGATGCCCTTGTTGGTAAt                                          >  2:553367/1‑150 (MQ=255)
                                                                                                                 tttttCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATc                                      <  2:1859657/150‑1 (MQ=255)
                                                                                                                   tttCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCgg                                    <  2:139230/150‑1 (MQ=255)
                                                                                                                     tCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTa                                  <  1:1196512/150‑1 (MQ=255)
                                                                                                                          tACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATg                             <  2:1581224/150‑1 (MQ=255)
                                                                                                                               gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGc                                                                                                     >  2:1720916/1‑73 (MQ=255)
                                                                                                                               gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGc                                                                                                     <  1:1720916/73‑1 (MQ=255)
                                                                                                                               gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCGTGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCttgtt                                                >  2:425569/1‑126 (MQ=255)
                                                                                                                               gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCGTGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCttgtt                                                <  1:425569/126‑1 (MQ=255)
                                                                                                                               gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCCTGGGTGTCAGTTGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCAGTACGATGTTCgg                        >  1:2113235/1‑150 (MQ=255)
                                                                                                                               gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCgg                        >  1:249914/1‑150 (MQ=255)
                                                                                                                               gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCgg                        >  1:2349264/1‑150 (MQ=255)
                                                                                                                                 tAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCt                      >  1:956985/1‑150 (MQ=255)
                                                                                                                                 tAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCt                      >  2:1767048/1‑150 (MQ=255)
                                                                                                                                 tAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCt                      >  2:796308/1‑150 (MQ=255)
                                                                                                                                   aCCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTgc                    <  1:796424/150‑1 (MQ=255)
                                                                                                                                    ccTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATg                                                                                      >  2:996735/1‑83 (MQ=255)
                                                                                                                                    ccTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATg                                                                                      <  1:996735/83‑1 (MQ=255)
                                                                                                                                    ccTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTgcg                   >  2:2414485/1‑150 (MQ=255)
                                                                                                                                     cTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCttgttg                                               >  2:1437054/1‑121 (MQ=255)
                                                                                                                                     cTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCttgttg                                               <  1:1437054/121‑1 (MQ=255)
                                                                                                                                      tGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGcgcg                 >  2:1937883/1‑150 (MQ=255)
                                                                                                                                       gCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTCGGTGGGAAGTATCGGTACGATGGTCGGCTGcgcgg                >  2:391628/1‑149 (MQ=255)
                                                                                                                                       gCAGTTCCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCAAGGCAGACAACAACCAATAACGATTTATGGCGAATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGCGCGa                <  1:1475174/150‑1 (MQ=255)
                                                                                                                                          cttGCAGGCCGACCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGAt                                                         <  2:1194742/105‑1 (MQ=255)
                                                                                                                                          cttGCAGGCCGACCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGAt                                                         >  1:1194742/2‑106 (MQ=255)
                                                                                                                                             gCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCttgt                                                 <  2:1817760/111‑1 (MQ=255)
                                                                                                                                             gCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCttgt                                                 >  1:1817760/1‑111 (MQ=255)
                                                                                                                                                    ggCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGCGCGACCCAGTCGCGGTa   <  2:2417518/150‑1 (MQ=255)
                                                                                                                                                     gCCTGCAAGGTACTGGTGGCTGACAGGTCGCGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGCGCGACCCAGTCGCGGTac  <  2:725140/150‑1 (MQ=255)
                                                                                                                                                     gCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGGGCGACCTAGTCGCGGTaa  >  2:606051/1‑149 (MQ=255)
                                                                                                                                                     gCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGCGCGACCCAGTCGCGGTac  <  1:1937883/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATCTGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGCGCGACCCAGTCGCGGTAC  >  NC_002947/1207543‑1207840

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: