Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 370,985 G→C 5.8% A155P (GCT→CCT)  dgcA → dimethylglycine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947370,9850GC 5.8% 382.9 / 10.6 121A155P (GCT→CCT) dgcAdimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base G (57/56);  new base C (3/4);  total (61/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 123 total aligned reads displayed.

CCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGG  >  NC_002947/370840‑371130
                                                                                                                                                 |                                                                                                                                                 
ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTgc                                                                                                                                               >  1:555469/1‑150 (MQ=255)
ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTgc                                                                                                                                               >  1:1605132/1‑150 (MQ=255)
ccACATGGGCCGTCGCTCGCGCGGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTgc                                                                                                                                               >  2:34290/1‑150 (MQ=255)
     tGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGCCTGCGCGta                                                                                                                                          <  1:2386881/150‑2 (MQ=255)
      gggCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTcg                                                                                                                                         >  2:87019/1‑150 (MQ=255)
      gggCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGCCTGCGCGtag                                                                                                                                         <  1:2263622/150‑3 (MQ=255)
                   cgcTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggc                                                                                                                            <  1:539516/150‑1 (MQ=255)
                    gcTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggcg                                                                                                                           <  2:39285/150‑1 (MQ=255)
                          aCGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGCGGTGGAAGAGATCTGGCGGGTAATCGGCACCTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGcc                                                                                                                        >  1:51228/1‑147 (MQ=255)
                          aCGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGcc                                                                                                                        <  2:51228/147‑1 (MQ=255)
                           cGGCTTCAACTGGCCGACGCTGATGTCGGCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGa                                                                                                                    <  2:18473/150‑1 (MQ=255)
                            ggCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGAc                                                                                                                   >  1:1001938/1‑150 (MQ=255)
                               ttCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGc                                                                                                                >  1:1321818/1‑150 (MQ=255)
                               ttCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGc                                                                                                                >  2:2036754/1‑150 (MQ=255)
                                 cAACTGGCCGACTCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCTCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCTAAAGAGGGCGGCCTGGACGGCGt                                                                                                              <  1:1056473/150‑1 (MQ=255)
                                 cAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGt                                                                                                              <  2:1001938/150‑1 (MQ=255)
                                  aaCTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg                                                                                                             <  1:2337047/150‑1 (MQ=255)
                                  aaCTGGCCGACGCTGATGTCGCCATAGGGCATTCGCGAACCTGTGCACCGCGCCACCTGAAAAACCATCGAGGTGGAAGAGATCTAGCGTGAAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCGAAGAGGGGGACCTGGACGGGGTg                                                                                                             >  2:1056473/1‑150 (MQ=255)
                                       gCCGACGCTGATGTCGCCATCGGGCAAGTTCCAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAtt                                                                                                        <  1:1210636/150‑1 (MQ=255)
                                         cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggcggc                                                                                                                         >  2:465840/1‑131 (MQ=255)
                                         cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggcggc                                                                                                                         <  1:465840/131‑1 (MQ=255)
                                         cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGt                                                                                                      >  2:1989158/1‑150 (MQ=255)
                                          gACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTc                                                                                                     >  2:2202415/1‑150 (MQ=255)
                                           aCGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCg                                                                                                    <  2:885175/150‑1 (MQ=255)
                                             gCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAgg                                                                                                                                                   >  2:1393269/1‑101 (MQ=255)
                                             gCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAgg                                                                                                                                                   <  1:1393269/101‑1 (MQ=255)
                                             gCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGc                                                                                                  >  1:1523351/1‑150 (MQ=255)
                                              cTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCa                                                                                                 <  1:87019/150‑1 (MQ=255)
                                               tGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAg                                                                                                >  1:1712528/1‑150 (MQ=255)
                                                 aTGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggcgg                                                                                                                          >  1:799143/1‑122 (MQ=255)
                                                 aTGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggcgg                                                                                                                          <  2:799143/122‑1 (MQ=255)
                                                     cGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAcc                                                                                          >  1:1264586/1‑150 (MQ=255)
                                                      gCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCt                                                                                                                       <  1:933551/120‑1 (MQ=255)
                                                      gCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCt                                                                                                                       >  2:933551/1‑120 (MQ=255)
                                                        cATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCACATCTCGGTACGGAGGACGTCTGACTTAAGAGGTCGGCCTGGACGGCGTGGAATCGTTGGCAGGGCACCAGc                                                                                       >  1:891897/1‑150 (MQ=255)
                                                        cATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCgcgc                                                                                                                                      <  1:272275/103‑1 (MQ=255)
                                                        cATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCgcgc                                                                                                                                      >  2:272275/1‑103 (MQ=255)
                                                        cATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAc                                                                                           <  2:64538/146‑1 (MQ=255)
                                                        cATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAc                                                                                           >  1:64538/1‑146 (MQ=255)
                                                         aTCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggc                                                                                                                            <  2:1036035/112‑1 (MQ=255)
                                                         aTCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggc                                                                                                                            >  1:1036035/1‑112 (MQ=255)
                                                          tCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGcac                                                                                     >  2:2185195/1‑150 (MQ=255)
                                                            gggCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACAt                                                                                   <  2:1712528/150‑1 (MQ=255)
                                                              gCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTg                                                                                              <  1:724285/137‑1 (MQ=255)
                                                              gCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTg                                                                                              >  2:724285/1‑137 (MQ=255)
                                                                aTTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATc                                                                               <  1:2036754/150‑1 (MQ=255)
                                                                aTTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATc                                                                               <  2:1264586/150‑1 (MQ=255)
                                                                   cgcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                            <  1:483488/150‑1 (MQ=255)
                                                                   cgcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATATGGCGGGTAATCGGCAAATACGCGCACGTTCCTTTTCGCGCCAAAGATGGCGGCCTGGACGGTGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                            <  2:275368/150‑1 (MQ=255)
                                                                        aCCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGtt                                                                       <  1:34290/150‑1 (MQ=255)
                                                                         ccTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTc                                                                      >  2:1768049/1‑150 (MQ=255)
                                                                          cTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCt                                                                     <  2:1605132/150‑1 (MQ=255)
                                                                           tgtgCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTg                                                                    >  1:1580705/1‑150 (MQ=255)
                                                                             tgCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGa                                                                  <  1:1768049/150‑1 (MQ=255)
                                                                             tgCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGa                                                                  <  2:587533/150‑1 (MQ=255)
                                                                                aCCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGCCTGCGCGta                                                                                                                                          >  2:1217916/1‑74 (MQ=255)
                                                                                aCCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGCCTGCGCGta                                                                                                                                          <  1:1217916/75‑2 (MQ=255)
                                                                                 ccGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                            >  2:1074890/1‑136 (MQ=255)
                                                                                 ccGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                            <  1:1074890/136‑1 (MQ=255)
                                                                                     gcCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACgg                                                          <  2:256681/150‑1 (MQ=255)
                                                                                       cACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGCCTgc                                                                                                                                               >  2:1695665/1‑63 (MQ=255)
                                                                                       cACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGCCTgc                                                                                                                                               <  1:1695665/63‑1 (MQ=255)
                                                                                          cTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAg                                                                         >  1:351899/1‑130 (MQ=255)
                                                                                          cTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTcaa                                                     >  1:993397/1‑150 (MQ=255)
                                                                                           tGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGACGGCGAGTGGCGGCAAAGAGGGCGGCCTGGACGGCGGGGAATTGTCGGGAGTGCACCAGCACATGGTCGACCATTTCTGGAGCCAACGGGTcaac                                                    >  2:68898/1‑150 (MQ=255)
                                                                                            gCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTcaaca                                                   >  1:2149794/1‑150 (MQ=255)
                                                                                             cAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCATCACATGATCGACCAGTTCTTGAGCCCACGGGTcaaaaa                                                  >  2:2309159/1‑150 (MQ=255)
                                                                                              aaaaCCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGGCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAg                                                 <  2:591068/150‑1 (MQ=255)
                                                                                               aaaCCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGc                                                <  1:2309159/150‑1 (MQ=255)
                                                                                                aaCCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTcaa                                                     <  1:1359006/144‑1 (MQ=255)
                                                                                                aaCCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTcaa                                                     >  2:1359006/1‑144 (MQ=255)
                                                                                                  ccATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAcc                                                                                          >  2:1331197/1‑105 (MQ=255)
                                                                                                   cATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGGGCGGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTAc                                            <  2:891897/150‑1 (MQ=255)
                                                                                                     tCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGc                                                                                             >  2:666022/1‑99 (MQ=255)
                                                                                                     tCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGc                                                                                             <  1:666022/99‑1 (MQ=255)
                                                                                                          gTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCa                                                                                         <  2:412780/98‑1 (MQ=255)
                                                                                                           tGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                            >  2:611379/1‑110 (MQ=255)
                                                                                                           tGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                            <  1:611379/110‑1 (MQ=255)
                                                                                                             gAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGc                                                <  1:2199669/136‑1 (MQ=255)
                                                                                                             gAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATgg                                  >  2:2215948/1‑150 (MQ=255)
                                                                                                             gAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATgg                                  >  2:1313512/1‑150 (MQ=255)
                                                                                                             gAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATgg                                  <  2:1597866/150‑1 (MQ=255)
                                                                                                               agagATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCt                                                                     >  1:2337245/1‑113 (MQ=255)
                                                                                                               agagATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCt                                                                     <  2:2337245/113‑1 (MQ=255)
                                                                                                               agagATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGAGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCAGATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATgggg                                >  1:2296898/1‑150 (MQ=255)
                                                                                                                gagaTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGggc                               <  2:1580705/150‑1 (MQ=255)
                                                                                                                gagaTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGggc                               <  1:2185195/150‑1 (MQ=255)
                                                                                                                   aTCTGGCGGGTAATCGGCAACTAGGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                            <  2:1890957/102‑1 (MQ=255)
                                                                                                                   aTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAc                                                                                           <  1:67576/87‑1 (MQ=255)
                                                                                                                   aTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGggcggc                            >  1:1070650/1‑150 (MQ=255)
                                                                                                                   aTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGggcggc                            >  1:1539786/1‑150 (MQ=255)
                                                                                                                   aTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGggcggc                            >  2:2272346/1‑150 (MQ=255)
                                                                                                                        gCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGTCCGTGCACCAGCAGATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCtt                       >  1:782203/1‑150 (MQ=255)
                                                                                                                            gTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGaa                   >  2:1301659/1‑150 (MQ=255)
                                                                                                                                   gCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCa            >  2:2196580/1‑150 (MQ=255)
                                                                                                                                   gCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCa            <  2:2149794/150‑1 (MQ=255)
                                                                                                                                      gcgaCGCGCAGCCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGa         >  1:825484/4‑150 (MQ=255)
                                                                                                                                           gcgcAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCttt                      >  1:960524/1‑132 (MQ=255)
                                                                                                                                            cgcAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCg   <  2:700836/150‑1 (MQ=255)
                                                                                                                                             gcAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCgg  >  2:2097440/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGG  >  NC_002947/370840‑371130

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