Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,047,336 C→G 5.2% L16V (CTG→GTG)  PP_2660 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,047,3360CG 5.2% 311.4 / 5.1 95L16V (CTG→GTG) PP_2660hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (47/43);  new base G (3/2);  total (50/45)
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.81e-01

CAGCTAAACCATGCAAGCTCAAGGTGGGCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGACCTGTTCGAT  >  NC_002947/3047190‑3047482
                                                                                                                                                  |                                                                                                                                                  
cAGCTAAACCATGCAAGCTCAAGGTGGGCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTgc                                                                                                                                                 <  1:1075836/150‑1 (MQ=255)
    tAAACCATGCAAGCTCAAGGTGGGCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTCACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGa                                                                                                                                             <  2:1464689/150‑1 (MQ=255)
      aaCCATGCAAGCTCAAGGTGGGCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAgg                                                                                                                                           >  1:1505787/1‑150 (MQ=255)
                  tCAAGGTGGGCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGtt                                                                                                                               <  2:331707/150‑1 (MQ=255)
                    aaGGTGGGCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGtttt                                                                                                                             <  2:1505787/150‑1 (MQ=255)
                          ggCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGTTGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGaa                                                                                                                       <  1:128229/150‑1 (MQ=255)
                           gCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAg                                                                                                                                            <  2:447907/128‑1 (MQ=255)
                           gCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAg                                                                                                                                            >  1:447907/1‑128 (MQ=255)
                           gCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAAc                                                                                                                      >  1:2158377/1‑150 (MQ=255)
                           gCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGCTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAAc                                                                                                                      >  1:1870121/1‑150 (MQ=255)
                            cGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACt                                                                                                                     >  2:1508510/1‑150 (MQ=255)
                                  ttGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTcagcca                                                                                                               >  2:1964204/1‑150 (MQ=255)
                                       gCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCAtt                                                                                                          <  1:1508510/150‑1 (MQ=255)
                                         cAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTc                                                                                                        <  2:1697728/150‑1 (MQ=255)
                                                 ttGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAAcc                                                                                                >  2:1249188/1‑150 (MQ=255)
                                                  tGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACcg                                                                                               >  2:1291271/1‑150 (MQ=255)
                                                   gCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACcgc                                                                                              <  2:1494575/150‑1 (MQ=255)
                                                    cAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACcgcg                                                                                             <  2:1957948/150‑1 (MQ=255)
                                                     aaCACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCGGGTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCgg                                                                                            >  1:2308656/1‑150 (MQ=255)
                                                         cccGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGAcc                                                                                        >  2:1325179/1‑150 (MQ=255)
                                                          ccGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTcagc                                                                                                                 >  2:853703/1‑124 (MQ=255)
                                                          ccGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTcagc                                                                                                                 <  1:853703/124‑1 (MQ=255)
                                                            ggCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGt                                                                                     >  1:1638720/1‑150 (MQ=255)
                                                             gCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGc                                                                                                                                        >  2:2143698/1‑98 (MQ=255)
                                                             gCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGc                                                                                                                                        <  1:2143698/98‑1 (MQ=255)
                                                             gCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGtt                                                                                    >  1:2124150/1‑150 (MQ=255)
                                                                   ttACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACcga                                                                              <  2:2045720/150‑1 (MQ=255)
                                                                     aCCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACcgacg                                                                            >  2:1305889/1‑150 (MQ=255)
                                                                      ccTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACcgacga                                                                           <  2:2124150/150‑1 (MQ=255)
                                                                      ccTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACcgacga                                                                           <  1:1964204/150‑1 (MQ=255)
                                                                       cTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAg                                                                                                                                            >  1:1008650/1‑84 (MQ=255)
                                                                       cTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAg                                                                                                                                            <  2:1008650/84‑1 (MQ=255)
                                                                       cTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAg                                                                          >  1:1596338/1‑150 (MQ=255)
                                                                        ttCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAgg                                                                         >  1:275371/1‑150 (MQ=255)
                                                                            gTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACcgcg                                                                                             >  2:1027529/1‑126 (MQ=255)
                                                                            gTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACcgcg                                                                                             <  1:1027529/126‑1 (MQ=255)
                                                                              tAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGtgctgc                                                                   <  1:1305889/150‑1 (MQ=255)
                                                                                        cctaGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGc                                                                                                                           >  1:679877/4‑84 (MQ=255)
                                                                                        cctaGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGc                                                                                                                           <  2:679877/81‑1 (MQ=255)
                                                                                        cccAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGtt                                                                                    >  2:1535228/1‑123 (MQ=255)
                                                                                        cccAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGtt                                                                                    <  1:1535228/123‑1 (MQ=255)
                                                                                        cccAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGt                                                         >  2:2193967/1‑150 (MQ=255)
                                                                                        cccAGGATTAACAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGGATGTGGt                                                         >  1:1878426/1‑150 (MQ=255)
                                                                                         ccAGGGTTAAAAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGtt                                                        >  2:2332346/1‑150 (MQ=255)
                                                                                         ccAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGtt                                                        >  1:1921748/1‑150 (MQ=255)
                                                                                         ccAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGAGGTAGAGCGCGCACTGCGCGAGGTGCAg                                                                                                                                      <  2:59304/72‑1 (MQ=255)
                                                                                         ccAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGAGGTAGAGCGCGCACTGCGCGAGGTGCAg                                                                                                                                      >  1:59304/1‑72 (MQ=255)
                                                                                           aGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGa                                                      >  1:938975/1‑150 (MQ=255)
                                                                                            ggATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGaa                                                     <  1:1599510/150‑1 (MQ=255)
                                                                                             gATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAg                                                    <  1:1249188/150‑1 (MQ=255)
                                                                                              aTTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAgg                                                   <  1:252168/150‑1 (MQ=255)
                                                                                                  ccAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTAt                                               >  1:1657952/1‑150 (MQ=255)
                                                                                                  ccAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTAt                                               >  1:822270/1‑150 (MQ=255)
                                                                                                  ccAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTAt                                               <  1:2039992/150‑1 (MQ=255)
                                                                                                  ccAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACCGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTAt                                               >  1:1855956/1‑150 (MQ=255)
                                                                                                   cAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATg                                              >  1:1757174/1‑150 (MQ=255)
                                                                                                        cgagcgaCGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTc                                                                                                                    <  2:1246803/75‑1 (MQ=255)
                                                                                                        cgagcgaCGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTc                                                                                                                    >  1:1246803/1‑75 (MQ=255)
                                                                                                              aCGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTgac                                   <  2:822270/150‑1 (MQ=255)
                                                                                                                gcgTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTgacga                                 <  2:1657952/150‑1 (MQ=255)
                                                                                                                  gtcgtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTgctgct                                                                  >  2:2138131/1‑115 (MQ=255)
                                                                                                                  gtcgtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTgctgct                                                                  <  1:2138131/115‑1 (MQ=255)
                                                                                                                   tcgtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCa                                                                                                            >  2:1277120/1‑72 (MQ=255)
                                                                                                                   tcgtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCa                                                                                                            <  1:1277120/72‑1 (MQ=255)
                                                                                                                   tcgtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACcgac                                                                             <  1:42881/103‑1 (MQ=255)
                                                                                                                   tcgtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACcgac                                                                             >  2:42881/1‑103 (MQ=255)
                                                                                                                    cgtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTAt                             <  2:2158377/150‑1 (MQ=255)
                                                                                                                     gtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTgg                                                          <  1:1237832/120‑1 (MQ=255)
                                                                                                                     gtcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTgg                                                          >  2:1237832/1‑120 (MQ=255)
                                                                                                                      tcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAg                                                    <  2:638829/125‑1 (MQ=255)
                                                                                                                      tcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAg                                                    >  1:638829/1‑125 (MQ=255)
                                                                                                                      tcTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGt                           >  2:1325971/1‑150 (MQ=255)
                                                                                                                       cTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAg                                                                          <  2:1746358/102‑1 (MQ=255)
                                                                                                                       cTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAg                                                                          >  1:1746358/1‑102 (MQ=255)
                                                                                                                          gATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAg                                                                          <  1:121994/99‑1 (MQ=255)
                                                                                                                          gATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAg                                                                          >  2:121994/1‑99 (MQ=255)
                                                                                                                             cTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGtcgctcgc                    <  1:1325179/150‑1 (MQ=255)
                                                                                                                                tCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGCTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGt                 >  2:915756/1‑150 (MQ=255)
                                                                                                                                  aGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGtgctgc                                                                   <  1:1418112/98‑1 (MQ=255)
                                                                                                                                  aGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGtgctgc                                                                   >  2:1418112/1‑98 (MQ=255)
                                                                                                                                   gCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTgg                                                          <  2:270139/106‑1 (MQ=255)
                                                                                                                                   gCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTgg                                                          >  1:270139/1‑106 (MQ=255)
                                                                                                                                      gTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAAt                                                              >  2:949436/1‑99 (MQ=255)
                                                                                                                                      gTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAAt                                                              <  1:949436/99‑1 (MQ=255)
                                                                                                                                        acctcgcgcAGTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGACCt         >  2:991910/5‑150 (MQ=255)
                                                                                                                                         cctcgcgcAGTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGAc           <  2:1090000/144‑1 (MQ=255)
                                                                                                                                         cctcgcgcAGTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGAc           >  1:1090000/4‑147 (MQ=255)
                                                                                                                                          ctcgcgcAGTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGaa                                                     <  2:2235015/102‑1 (MQ=255)
                                                                                                                                          ctcgcgcAGTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGaa                                                     >  1:2235015/3‑104 (MQ=255)
                                                                                                                                          agCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGACCTGt       <  2:1757174/150‑1 (MQ=255)
                                                                                                                                           gcgcgcACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGACCTGtt      >  1:1882811/1‑150 (MQ=255)
                                                                                                                                            cgcgcACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCTTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGACCTGTTc     <  2:1882811/150‑1 (MQ=255)
                                                                                                                                               gcACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGtcg                         <  2:1864440/127‑1 (MQ=255)
                                                                                                                                               gcACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGtcg                         >  1:1864440/1‑127 (MQ=255)
                                                                                                                                               gcACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGACCTGTTCGAt  >  2:27457/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGCTAAACCATGCAAGCTCAAGGTGGGCGGGTGTTGCAGCCAGGAAGGTTGCAACACCCGGCCGCGTTACCTTCCGTTAAAGTGTTACCCAGGATTACCAATGCGAGCGACGCGTCGTCTTGATCTGTCAGCGGTAGAGCGCGCACTGCGCGAGGTGCAGGGCCGTTTTGCCGAACTCAGCCAGCATTTCACCGAACCGCGGGACCCGTTTACCGACGAGGTGCTGCTGAATGTGGTTGAAGGCTATGCGCTGATTGACGACTATGTCGCTCGCGGTATCGACCTGTTCGAT  >  NC_002947/3047190‑3047482

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