Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,765,286 G→T 5.4% G137C (GGC→TGC)  PP_1574 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,765,2860GT 5.4% 692.9 / 20.7 202G137C (GGC→TGC) PP_1574hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (107/84);  new base T (7/4);  total (114/88)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 208 total aligned reads displayed.

CGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGCC  >  NC_002947/1765137‑1765435
                                                                                                                                                     |                                                                                                                                                     
cGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTg                                                                                                                                                       <  2:1381879/150‑1 (MQ=255)
cGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTg                                                                                                                                                       <  1:4051991/150‑1 (MQ=255)
cGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTg                                                                                                                                                       <  2:276866/150‑1 (MQ=255)
cGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTg                                                                                                                                                       >  1:2764305/1‑150 (MQ=255)
  cTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGc                                                                                                                                                     >  1:3825657/1‑150 (MQ=255)
     gACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACg                                                                                                                                                  <  1:1955052/150‑1 (MQ=255)
      aCCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGa                                                                                                                                                 <  1:3622181/150‑1 (MQ=255)
       cccGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAg                                                                                                                                                <  2:4001667/150‑1 (MQ=255)
          ggTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTc                                                                                                                                             <  2:1337397/150‑1 (MQ=255)
             tGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTc                                                                                                                                          <  2:1172935/150‑1 (MQ=255)
             tGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCGTGCCAGAGCTCGTc                                                                                                                                          <  2:4450094/150‑1 (MQ=255)
               aaCTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCGTGCCAGAGCTCGTcga                                                                                                                                        <  1:2829296/150‑3 (MQ=255)
                 cTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCctct                                                                                                                                      <  2:3006183/150‑1 (MQ=255)
                   ggAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTcc                                                                                                                                    >  2:3807865/1‑150 (MQ=255)
                    gagaAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCg                                                                                                                                   >  2:2408589/1‑150 (MQ=255)
                         ggACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATc                                                                                                                              >  2:1476939/1‑150 (MQ=255)
                         ggACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCGATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATc                                                                                                                              >  2:4301422/1‑150 (MQ=255)
                              cGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCGGGGACAGACGAAGAGCTCACCGAGCTCGCCGAGCTCTGGCACGAGCTCGTGCTCTCCGGCCTcggcgg                                                                                                                         >  2:4037065/1‑150 (MQ=255)
                                gagaGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCt                                                                                                                       <  2:2185285/150‑1 (MQ=255)
                                   agCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcga                                                                                                                    <  1:2894465/150‑1 (MQ=255)
                                     cACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcgacg                                                                                                                  >  1:2991687/1‑150 (MQ=255)
                                      aCCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcgacga                                                                                                                 <  2:3691717/150‑1 (MQ=255)
                                        cAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATc                                                                                                               <  2:3482344/150‑1 (MQ=255)
                                        cAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATc                                                                                                               >  1:4034229/1‑150 (MQ=255)
                                        cAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATc                                                                                                               <  1:3836659/150‑1 (MQ=255)
                                         aaGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATcg                                                                                                              >  1:1555827/1‑150 (MQ=255)
                                           gcgcCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATcgcg                                                                                                            <  2:3772364/150‑1 (MQ=255)
                                              ccTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAg                                                                                                         <  2:1601373/150‑1 (MQ=255)
                                                tGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGc                                                                                                       >  2:3545383/1‑150 (MQ=255)
                                                    aTCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCaaa                                                                                                   <  1:1712754/150‑1 (MQ=255)
                                                    aTCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCaaa                                                                                                   >  2:334207/1‑150 (MQ=255)
                                                        ctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGTCGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGaaa                                                                                               >  2:3020503/1‑150 (MQ=255)
                                                         ttGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGaaa                                                                                               <  2:1366522/149‑1 (MQ=255)
                                                          tGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAAc                                                                                             >  1:2908442/1‑150 (MQ=255)
                                                           gATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAAcc                                                                                            <  2:2288749/150‑1 (MQ=255)
                                                             tGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTg                                                                                          >  1:342100/1‑150 (MQ=255)
                                                              ggCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                         <  2:3069888/150‑1 (MQ=255)
                                                              ggCGACCTGTCCGTCGCGCTGCTGGTAGCTAGTCGGGAGGGCAACGACCTGGGAAAGACGACGAGCTCACCGAGCTCGACGAGCTCGGGCACGAGCTCGTCCTCTCCGTCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                         <  2:2146664/150‑1 (MQ=255)
                                                                cGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCG‑‑‑‑‑‑‑‑TGCCA‑GAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAg                                                                              <  1:3939120/150‑1 (MQ=255)
                                                                     tGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCg                                                                                  <  1:2122805/150‑1 (MQ=255)
                                                                     tGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCg                                                                                  >  1:1494294/1‑150 (MQ=255)
                                                                       tCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTg                                                                                >  2:4086545/1‑150 (MQ=255)
                                                                        ccGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGa                                                                               <  2:2838614/150‑1 (MQ=255)
                                                                          gTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGt                                                                             <  2:1273579/150‑1 (MQ=255)
                                                                           tCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAg                                                                                                         >  1:1824339/1‑121 (MQ=255)
                                                                             gcgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCa                                                                          <  2:2648802/150‑1 (MQ=255)
                                                                             gcgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCa                                                                          <  2:2012653/150‑1 (MQ=255)
                                                                              cgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAg                                                                         <  2:2319308/150‑1 (MQ=255)
                                                                              cgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAg                                                                         <  2:179139/150‑1 (MQ=255)
                                                                                ctgctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGt                                                                       >  1:94132/1‑150 (MQ=255)
                                                                                 tgctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTc                                                                      >  1:1016748/1‑150 (MQ=255)
                                                                                    tgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTg                                                                   <  1:718085/150‑1 (MQ=255)
                                                                                     ggTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTgg                                                                  <  1:3811017/150‑1 (MQ=255)
                                                                                      gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                 >  2:1039701/1‑150 (MQ=255)
                                                                                      gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                 >  2:300676/1‑150 (MQ=255)
                                                                                      gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                 >  2:1039782/1‑150 (MQ=255)
                                                                                      gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                 <  1:1916237/150‑1 (MQ=255)
                                                                                         gCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGgaga                                                              >  2:779888/1‑150 (MQ=255)
                                                                                         gCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGgaga                                                              >  1:3161873/1‑150 (MQ=255)
                                                                                          cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                         <  1:2596694/122‑1 (MQ=255)
                                                                                          cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                         >  2:2596694/1‑122 (MQ=255)
                                                                                          cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                             <  1:134753/150‑1 (MQ=255)
                                                                                          cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                             <  1:4225329/150‑1 (MQ=255)
                                                                                          cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                             <  1:2253566/150‑1 (MQ=255)
                                                                                             ttCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTa                                                          <  1:820312/150‑1 (MQ=255)
                                                                                             ttCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTa                                                          >  2:455295/1‑150 (MQ=255)
                                                                                             ttCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTa                                                          >  1:2554953/1‑150 (MQ=255)
                                                                                                gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAg                                                       >  1:1263837/1‑150 (MQ=255)
                                                                                                gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAg                                                       >  2:2134946/1‑150 (MQ=255)
                                                                                                gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAg                                                       <  1:2961237/150‑1 (MQ=255)
                                                                                                 ggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAgg                                                      <  2:601419/150‑1 (MQ=255)
                                                                                                        aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                               >  2:1660766/1‑150 (MQ=255)
                                                                                                        aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                               >  2:3439796/1‑150 (MQ=255)
                                                                                                        aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                               >  1:4387105/1‑150 (MQ=255)
                                                                                                        aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                               >  1:52243/1‑150 (MQ=255)
                                                                                                              cTGGGAAAGACGAAGAGCTAACCGAGCTCGACGAGCTCTGGCACGAGCTCGACCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCt                                         >  2:3270422/1‑150 (MQ=255)
                                                                                                               tGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTc                                        >  2:2874951/1‑150 (MQ=255)
                                                                                                                  gAAAGACGAAGAGCTCACCGAGCTCGTCGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCt                                     <  1:3020503/150‑1 (MQ=255)
                                                                                                                  gAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCt                                     >  1:2377524/1‑150 (MQ=255)
                                                                                                                  gAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCt                                     >  2:3952302/1‑150 (MQ=255)
                                                                                                                    aaGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCCTCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGGGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGc                                   >  2:1094368/1‑150 (MQ=255)
                                                                                                                      gACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCAt                                 >  2:1030441/1‑150 (MQ=255)
                                                                                                                         gAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTc                              >  1:739485/1‑150 (MQ=255)
                                                                                                                         gAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTc                              >  1:3008949/1‑150 (MQ=255)
                                                                                                                           agagCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCgg                            >  2:2626645/1‑150 (MQ=255)
                                                                                                                           agagCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCgg                            >  1:949224/1‑150 (MQ=255)
                                                                                                                            gagCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCggg                           <  2:1605535/150‑1 (MQ=255)
                                                                                                                             agCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGa                          <  1:4314553/150‑1 (MQ=255)
                                                                                                                                tCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATgc                       <  2:1057195/150‑1 (MQ=255)
                                                                                                                                  aCCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATgcgc                     >  1:4426799/1‑150 (MQ=255)
                                                                                                                                    cGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGt                   >  1:162878/1‑150 (MQ=255)
                                                                                                                                         tCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgc              >  2:522039/1‑150 (MQ=255)
                                                                                                                                           gacgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc            >  1:1160359/1‑150 (MQ=255)
                                                                                                                                           gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc            >  2:2513599/1‑150 (MQ=255)
                                                                                                                                           gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc            >  2:573153/1‑150 (MQ=255)
                                                                                                                                           gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc            >  1:412064/1‑150 (MQ=255)
                                                                                                                                           gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc            >  2:1912136/1‑150 (MQ=255)
                                                                                                                                                   ctGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACggg    <  1:802615/150‑1 (MQ=255)
                                                                                                                                                    tGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGc   <  1:1078173/150‑1 (MQ=255)
                                                                                                                                                     ggCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGcc  <  2:3900895/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGCC  >  NC_002947/1765137‑1765435

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