Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 370,985 G→C 5.2% 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.2% 398.3 / 8.8 127A155P (GCT→CCT) dgcAdimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base G (59/61);  new base C (2/5);  total (61/66)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 128 total aligned reads displayed.

CCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTT  >  NC_002947/370840‑371134
                                                                                                                                                 |                                                                                                                                                     
ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTgc                                                                                                                                                   >  2:701126/1‑150 (MQ=255)
ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTgc                                                                                                                                                   >  2:186/1‑150 (MQ=255)
ccACATGGGCAGTAGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTgc                                                                                                                                                   >  2:392802/1‑150 (MQ=255)
  acaTGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCTCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTgcgc                                                                                                                                                 >  1:1635751/1‑150 (MQ=255)
  acaTGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCAGCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTgcgc                                                                                                                                                 >  2:88482/1‑150 (MQ=255)
   caTGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCAGCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGcgcg                                                                                                                                                >  2:1888252/1‑150 (MQ=255)
     tGGGCCGTCGCTCGCGCTTGTACGGCTTCAACTGGCCGACGCTGATGTCGCCCTCGGGCATTCGCGAACCCGTGCACCGCCCCACCTGCAAACCCCTCGAGGTGGAAGAGGTCGGTCGGGTAATCGGCATCTACGCGCAGCCTGCGCGta                                                                                                                                              <  2:313730/150‑2 (MQ=255)
     tGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGCCTGCGCGta                                                                                                                                              <  1:1814060/150‑2 (MQ=255)
      gggCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGCCTGCGCGtag                                                                                                                                             <  2:2278470/150‑3 (MQ=255)
               acggcgcTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAaga                                                                                                                                    >  1:82839/4‑150 (MQ=255)
                tcgcGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAagag                                                                                                                                   >  1:1010873/1‑150 (MQ=255)
                    gcTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggcg                                                                                                                               <  1:2014789/150‑1 (MQ=255)
                       gggACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCTGCGCGCCAAAGAGGGCGGcc                                                                                                                            <  1:1255710/150‑1 (MQ=255)
                        ggACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGCGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCt                                                                                                                           <  2:2337294/150‑1 (MQ=255)
                        ggACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCt                                                                                                                           >  1:419248/1‑150 (MQ=255)
                        ggACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCt                                                                                                                           >  2:1894970/1‑150 (MQ=255)
                               ttCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGc                                                                                                                    >  2:1626727/1‑150 (MQ=255)
                                       gCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAtt                                                                                                            <  2:1962760/150‑1 (MQ=255)
                                       gCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCCCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAtt                                                                                                            >  1:259343/1‑150 (MQ=255)
                                          gACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTc                                                                                                         <  2:579549/150‑1 (MQ=255)
                                              cTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCa                                                                                                     <  1:1161888/150‑1 (MQ=255)
                                              cTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCa                                                                                                     <  2:419248/150‑1 (MQ=255)
                                               tGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGATGGCGTGGAATTGTCGGCAg                                                                                                    >  1:674556/1‑150 (MQ=255)
                                               tGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAg                                                                                                    >  2:1349689/1‑150 (MQ=255)
                                                 aTGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTg                                                                                                  <  1:849734/150‑1 (MQ=255)
                                                    tCGCCATCGGGCATTCGCGAACCTGTGGACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAc                                                                                               <  1:186/150‑1 (MQ=255)
                                                    tCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGCTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAACGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAc                                                                                               >  2:1071607/1‑150 (MQ=255)
                                                    tCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAc                                                                                               <  1:695059/150‑1 (MQ=255)
                                                    tCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAc                                                                                               <  2:2008244/150‑1 (MQ=255)
                                                      gCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGGGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCa                                                                                             <  2:1280449/150‑1 (MQ=255)
                                                         aTCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGca                                                                                          <  2:1169294/150‑1 (MQ=255)
                                                           cGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGcaca                                                                                        <  1:1349689/150‑1 (MQ=255)
                                                               cATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACANGAt                                                                                    <  2:1507910/150‑1 (MQ=255)
                                                                aTTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATc                                                                                   >  1:2040020/1‑150 (MQ=255)
                                                                aTTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATc                                                                                   <  1:392802/150‑1 (MQ=255)
                                                                      gAACCTGTGCACCGCACCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGaa                                                                                                              >  1:199075/1‑117 (MQ=255)
                                                                      gAACCTGTGCACCGCACCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGaa                                                                                                              <  2:199075/117‑1 (MQ=255)
                                                                        aCCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCGACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                                >  1:1658537/1‑145 (MQ=255)
                                                                          ctatgCAACGCGCCAACTGCAAAACCATCGAGGTGGAAGAGATCTGGAGGGTAATCGGCAACTACGCGCAGGCTGCTAGTCACGCCAAAGAGGGCGGCCTGGACGGCGTTGAATTGTCGGCAGTGTACCAGCACATGATCGACCAGTTCt                                                                         <  1:558644/147‑1 (MQ=255)
                                                                          cTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGt                                                                                                                  >  1:809031/1‑109 (MQ=255)
                                                                          cTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTgg                                                                                                                <  1:264445/111‑1 (MQ=255)
                                                                          cTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTgg                                                                                                                >  2:264445/1‑111 (MQ=255)
                                                                          cTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCt                                                                         <  2:1573196/150‑1 (MQ=255)
                                                                             tgCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGa                                                                      <  1:1894970/150‑1 (MQ=255)
                                                                                aCCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAtt                                                                                                            >  1:2118737/1‑109 (MQ=255)
                                                                                aCCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAtt                                                                                                            <  2:2118737/109‑1 (MQ=255)
                                                                                     gcCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACgg                                                              <  2:2251841/150‑1 (MQ=255)
                                                                                     gcCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGAACGCCTTGGAATTGTCGTCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACgg                                                              <  2:82839/150‑1 (MQ=255)
                                                                                      ccACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTg                                                                                                                          <  2:1560400/89‑1 (MQ=255)
                                                                                      ccACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTg                                                                                                                          >  1:1560400/1‑89 (MQ=255)
                                                                                      ccACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCt                                                                         <  2:1692612/138‑1 (MQ=255)
                                                                                      ccACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCt                                                                         >  1:1692612/1‑138 (MQ=255)
                                                                                        aCCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAc                                                                                               <  1:1944867/114‑1 (MQ=255)
                                                                                        aCCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACgg                                                              >  2:659995/1‑147 (MQ=255)
                                                                                        aCCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACgg                                                              <  1:659995/147‑1 (MQ=255)
                                                                                        aCCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTc                                                           <  2:1649857/150‑1 (MQ=255)
                                                                                          cTGCCAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGCGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGATCAGTTCTGGAGCCCACGGGTcaa                                                         <  2:1188123/150‑1 (MQ=255)
                                                                                                   cATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGCCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGATCCCTTGGGTCGACAAGCGTAc                                                >  2:723862/1‑150 (MQ=255)
                                                                                                    aTCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGATGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTAcc                                               >  2:1660667/1‑150 (MQ=255)
                                                                                                    aTCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                                >  2:1198118/1‑117 (MQ=255)
                                                                                                    aTCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc                                                                                <  1:1198118/117‑1 (MQ=255)
                                                                                                    aTCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTAcc                                               >  1:1858938/1‑150 (MQ=255)
                                                                                                     tCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGt                                                                                                          >  2:1708761/1‑90 (MQ=255)
                                                                                                     tCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGt                                                                                                          <  1:1708761/90‑1 (MQ=255)
                                                                                                     tCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACcg                                              >  1:138737/1‑150 (MQ=255)
                                                                                                      cGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACcga                                             <  1:1093661/150‑1 (MQ=255)
                                                                                                        aGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTcaacaa                                                      <  2:1886457/139‑1 (MQ=255)
                                                                                                           tGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAAt                                        >  2:819222/1‑150 (MQ=255)
                                                                                                            ggAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATg                                       <  2:1635751/150‑1 (MQ=255)
                                                                                                             gAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATgg                                      >  1:1037448/1‑150 (MQ=255)
                                                                                                              aaGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATggg                                     >  2:1419242/1‑150 (MQ=255)
                                                                                                               agagATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTcaacaa                                                      <  2:1368580/132‑1 (MQ=255)
                                                                                                               agagATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATgggg                                    >  1:1622826/1‑150 (MQ=255)
                                                                                                                gagaTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGggc                                   <  2:576310/150‑1 (MQ=255)
                                                                                                                  gaTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGggcgg                                 >  1:1491769/1‑150 (MQ=255)
                                                                                                                   aTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGa                                                                                                                        <  1:34600/62‑1 (MQ=255)
                                                                                                                   aTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGa                                                                                                                        >  2:34600/1‑62 (MQ=255)
                                                                                                                   aTCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGggcggc                                <  1:2203232/150‑1 (MQ=255)
                                                                                                                        gCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCtt                           >  2:747385/1‑150 (MQ=255)
                                                                                                                          gggTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATggg                                     <  1:2130518/138‑1 (MQ=255)
                                                                                                                          gggTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATggg                                     >  2:2130518/1‑138 (MQ=255)
                                                                                                                           ggTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGa                        >  2:734774/1‑150 (MQ=255)
                                                                                                                            gTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTAc                                                >  1:2210493/1‑125 (MQ=255)
                                                                                                                            gTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTAc                                                <  2:2210493/125‑1 (MQ=255)
                                                                                                                            gTAATCGGAACCTACGCGTAGGCTGCGCGGCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCTCACGGGTCAACAAGCGTACCAACGAATGGGGCGGCACCTTTGaa                       <  2:1098555/150‑1 (MQ=255)
                                                                                                                             tAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAg                      <  2:2417317/150‑1 (MQ=255)
                                                                                                                             tAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAg                      <  1:1419242/150‑1 (MQ=255)
                                                                                                                               aTCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGc                    >  1:1809037/1‑150 (MQ=255)
                                                                                                                                   gCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAg                                                                                            <  1:1281867/74‑1 (MQ=255)
                                                                                                                                   gCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCgaa                >  2:1449450/1‑148 (MQ=255)
                                                                                                                                     aaCTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATg              <  2:1622826/150‑1 (MQ=255)
                                                                                                                                        tACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAg                      >  1:2043788/1‑139 (MQ=255)
                                                                                                                                        tACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAg                      <  2:2043788/139‑1 (MQ=255)
                                                                                                                                        tACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGaaa           >  1:950251/1‑150 (MQ=255)
                                                                                                                                             gcAGCCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAcc                                                                                              <  2:2449187/62‑1 (MQ=255)
                                                                                                                                             gcAGCCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCAcc                                                                                              >  1:2449187/1‑62 (MQ=255)
                                                                                                                                              cAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGGCGGCGTGGAGTTGTCGGCAGTGCACCAGCACATGATCGACCACTTCTGGAGCCCACGGGTCAACTAGCGTACCGACGAATGGGGCCGCACCTTTGAAGGCCGCATGAAATTCGGc     <  2:2118599/150‑1 (MQ=255)
                                                                                                                                              cAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGc     <  2:1037448/150‑1 (MQ=255)
                                                                                                                                              cAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGGGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGc     <  1:1888252/150‑1 (MQ=255)
                                                                                                                                                 gCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCCCCTTTGAAGGCCGCATGAAATTCGGCCtt  <  1:1626727/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGACCAGTTCTGGAGCCCACGGGTCAACAAGCGTACCGACGAATGGGGCGGCACCTTTGAAGGCCGCATGAAATTCGGCCTT  >  NC_002947/370840‑371134

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