Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,207,688 T→C 6.1% D440G (GAC→GGC)  PP_1058 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,207,6880TC 6.1% 287.8 / 6.3 98D440G (GAC→GGC) PP_1058hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (47/45);  new base C (2/4);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

CATCGGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGCGCGA  >  NC_002947/1207539‑1207826
                                                                                                                                                      |                                                                                                                                          
cATCGGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGt                                                                                                                                            <  2:2426709/150‑1 (MQ=255)
 aTCGGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTc                                                                                                                                           >  1:2233617/1‑150 (MQ=255)
 aTCGGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTc                                                                                                                                           >  2:1735318/1‑150 (MQ=255)
  tCGGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCg                                                                                                                                          <  2:538960/150‑1 (MQ=255)
  tCGGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCg                                                                                                                                          <  2:538959/150‑1 (MQ=255)
    ggcggTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTGACTCCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGc                                                                                                                                        <  2:63768/150‑1 (MQ=255)
          gCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCaa                                                                                                                                  <  2:303024/150‑1 (MQ=255)
          gCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGCCGACCTGCaa                                                                                                                                  <  2:1395264/150‑1 (MQ=255)
           cGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTCAGCGTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGCCGACCTGCaac                                                                                                                                 <  1:1842643/150‑2 (MQ=255)
                gAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACt                                                                                                                            >  2:1269757/1‑150 (MQ=255)
                gAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACt                                                                                                                            >  1:1754032/1‑150 (MQ=255)
                 aaCAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtg                                                                                                                           >  2:2431534/1‑150 (MQ=255)
                  aCAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtgg                                                                                                                          <  1:784077/150‑1 (MQ=255)
                   cAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtggt                                                                                                                         <  1:1275632/150‑1 (MQ=255)
                     aGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtggtgg                                                                                                                       <  2:2233617/150‑1 (MQ=255)
                      ggTGAAACTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCCACTGTATGTCAGTC‑GTGTGCTCTTTTGCTGCAACCCCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGTCCTGCAAGGTACTGGTGGc                                                                                                                      <  1:280508/150‑1 (MQ=255)
                       gTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGc                                                                                                                      <  1:2364784/149‑1 (MQ=255)
                       gTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGc                                                                                                                      >  2:2364784/1‑149 (MQ=255)
                       gTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCt                                                                                                                     <  1:2348734/150‑1 (MQ=255)
                           aGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACa                                                                                                                 <  2:1899874/150‑1 (MQ=255)
                           aGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACa                                                                                                                 <  2:1849033/150‑1 (MQ=255)
                              tGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGt                                                                                                              <  1:1266184/150‑1 (MQ=255)
                               gTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTc                                                                                                             <  1:2328993/150‑1 (MQ=255)
                                  gCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAg                                                                                                          >  2:1196512/1‑150 (MQ=255)
                                   cTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAgg                                                                                                         <  1:1289031/150‑1 (MQ=255)
                                   cTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAgg                                                                                                         <  2:1413102/150‑1 (MQ=255)
                                     cACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCg                                                                                                       <  1:708011/150‑1 (MQ=255)
                                        cTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTc                                                                                                    >  1:2357778/1‑150 (MQ=255)
                                        cTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TATCCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTc                                                                                                    <  1:746190/150‑1 (MQ=255)
                                          ggTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGa                                                                                                  <  2:407305/150‑1 (MQ=255)
                                              ggtggtCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtgg                                                                                                                          >  1:2427308/1‑122 (MQ=255)
                                              ggtggtCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACtgg                                                                                                                          <  2:2427308/122‑1 (MQ=255)
                                              ggtggtCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTc                                                                                              >  1:2026633/1‑150 (MQ=255)
                                                     ttGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGc                                                                                       <  1:2431534/150‑1 (MQ=255)
                                                     ttGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGc                                                                                       >  1:2447376/1‑150 (MQ=255)
                                                     ttGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGGTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGc                                                                                       >  2:185957/1‑150 (MQ=255)
                                                        gCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAg                                                                                                          <  2:1527489/128‑1 (MQ=255)
                                                        gCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAg                                                                                                          >  1:1527489/1‑128 (MQ=255)
                                                        gCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCggg                                                                                    >  1:188949/1‑150 (MQ=255)
                                                         cAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGTAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGc                                                                                   <  1:497391/150‑1 (MQ=255)
                                                           gTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGa                                                                                 <  1:593255/150‑1 (MQ=255)
                                                            tGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGCCGACCTGCaac                                                                                                                                 <  2:786025/101‑2 (MQ=255)
                                                            tGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGCCGACCTGCaac                                                                                                                                 >  1:786025/1‑100 (MQ=255)
                                                              agcagcCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTg                                                                              >  1:2417518/1‑150 (MQ=255)
                                                               gcagcCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGa                                                                             <  1:1735318/150‑1 (MQ=255)
                                                                cagcCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAg                                                                            >  1:316232/1‑150 (MQ=255)
                                                                      cTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGc                                                                      <  1:362691/150‑1 (MQ=255)
                                                                      cTCGCCGATCAGGCCGATC‑TGTGCGGAAAATTCAGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGc                                                                      <  2:188949/150‑1 (MQ=255)
                                                                         gCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGa                                                                                                  <  1:1151511/119‑1 (MQ=255)
                                                                         gCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGa                                                                                                  >  2:1151511/1‑119 (MQ=255)
                                                                            gATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAg                                                                >  1:139230/1‑150 (MQ=255)
                                                                            gATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGACTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAg                                                                >  1:1859657/1‑150 (MQ=255)
                                                                             aTCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAgg                                                               <  2:540033/150‑1 (MQ=255)
                                                                              tCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGc                                                              >  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                                                                                                                     >  1:185252/1‑80 (MQ=255)
                                                                                              cGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCt                                                                                                                     <  2:185252/80‑1 (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:1060400/1‑150 (MQ=255)
                                                                                                                cAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTaa                             <  2:316232/150‑1 (MQ=255)
                                                                                                                 aaGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGAtt                                                                               >  2:2272273/1‑99 (MQ=255)
                                                                                                                 aaGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGAtt                                                                               <  1:2272273/99‑1 (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                                  <  1:425569/126‑1 (MQ=255)
                                                                                                                                   gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCGTGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCttgtt                                  >  2:425569/1‑126 (MQ=255)
                                                                                                                                   gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCCTGGGTGTCAGTTGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCAGTACGATGTTCgg          >  1:2113235/1‑150 (MQ=255)
                                                                                                                                   gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCgg          >  1:249914/1‑150 (MQ=255)
                                                                                                                                   gCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCgg          >  1:2349264/1‑150 (MQ=255)
                                                                                                                                     tAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCt        >  2:1767048/1‑150 (MQ=255)
                                                                                                                                     tAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCt        >  2:796308/1‑150 (MQ=255)
                                                                                                                                     tAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCt        >  1:956985/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)
                                                                                                                                                      |                                                                                                                                          
CATCGGCGGTGCGTTCGAACAAGGTGAAGCTGTAGCTCACCTGGTCGGTGGTCTTGGCAGTGAGCAGCCGCTCGCCGATCAGGCCGATC‑TGTGCGGCAAACTCCGGGTCGGCAAGGTTTTTCAGGTACTGGCTAACCTGCAGTTGCAGGTCGGCCTGCAAGGTACTGGTGGCTGACAGGTCGAGGCGGTCGAGGTCGTACAGCGGGCGATTGAGCATGGCCGACAGGCGGTTGCGCGCCAGGCTGATGCCCTTGTTGGTAATGATCGGTACGATGGTCGGCTGCGCGA  >  NC_002947/1207539‑1207826

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: ‑