Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,161,034 G→C 5.3% intergenic (+16/‑193) PP_3661 → / → PP_3662 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,161,0340GC 5.3% 525.9 / 14.8 152intergenic (+16/‑193)PP_3661/PP_3662membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/100);  new base C (2/6);  total (46/106)
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.93e-01

Only 100 of 165 total aligned reads displayed.

TGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGG  >  NC_002947/4160885‑4161183
                                                                                                                                                     |                                                                                                                                                      
tGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCATGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCgg                                                                                                                                                        <  2:1653205/150‑1 (MQ=255)
 gggCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCggg                                                                                                                                                       >  2:3528579/1‑149 (MQ=255)
 gggCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGgc                                                                                                                                                       >  2:1996508/1‑150 (MQ=255)
  ggCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGgct                                                                                                                                                      <  1:3632869/150‑1 (MQ=255)
       ggAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTc                                                                                                                                                 >  2:2035990/1‑150 (MQ=255)
         agaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTcgc                                                                                                                                               <  1:2378678/150‑1 (MQ=255)
         agaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTcgc                                                                                                                                               <  1:2378673/150‑1 (MQ=255)
         agaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTcgc                                                                                                                                               <  1:238256/150‑1 (MQ=255)
         agaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTcgc                                                                                                                                               <  1:2495165/150‑1 (MQ=255)
         agaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTcgc                                                                                                                                               <  2:502965/150‑1 (MQ=255)
          gaCCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTcgcg                                                                                                                                              <  1:3040644/150‑1 (MQ=255)
          gaCCCGCTTCAGCCAGATTCGACAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTagcg                                                                                                                                              >  2:4486226/1‑150 (MQ=255)
          gaCCCGCGTCAGCCCCAGGCGCCAGGCAGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTcgcg                                                                                                                                              <  2:1003323/150‑1 (MQ=255)
             ccGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGt                                                                                                                                           >  2:2173318/1‑150 (MQ=255)
               gCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGGAA                                                                                                                                         >  2:4024546‑M1/1‑150 (MQ=255)
                       cAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTC                                                                                                                                 <  1:1153071‑M1/150‑1 (MQ=255)
                       cAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTC                                                                                                                                 <  1:895102‑M1/150‑1 (MQ=255)
                             cGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAG                                                                                                                           <  1:2757807‑M1/150‑1 (MQ=255)
                                cAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCGTCGCGGGTAAAACCCCTCCGACAggtg                                                                                                                        >  2:2328033/1‑150 (MQ=255)
                                 aGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACAGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgg                                                                                                                       <  2:2074203/150‑1 (MQ=255)
                                    cTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAggtgggtg                                                                                                                    <  1:115517/150‑1 (MQ=255)
                                     tGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGtgc                                                                                                                   <  1:1279164/150‑1 (MQ=255)
                                           aCGGGCGCTGGCCACCGGTGCGATCCTCACCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCt                                                                                                             <  1:2470535/150‑1 (MQ=255)
                                            cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTg                                                                                                            <  2:72833/150‑1 (MQ=255)
                                            cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTg                                                                                                            <  2:1277212/150‑1 (MQ=255)
                                            cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTg                                                                                                            <  2:72828/150‑1 (MQ=255)
                                             gggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGa                                                                                                           <  1:4358675/150‑1 (MQ=255)
                                              ggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                                          <  2:1566321/150‑1 (MQ=255)
                                              ggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                                          <  2:3540629/150‑1 (MQ=255)
                                              ggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGaa                                                                                                          <  2:244521/150‑1 (MQ=255)
                                               gcgcTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAg                                                                                                         <  2:3385944/150‑1 (MQ=255)
                                                   tGGCCACCGGTGCGATCCTCACCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACgg                                                                                                     <  1:3419738/150‑1 (MQ=255)
                                                   tGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGGTGGGTCGCGGGGGAAAACCCACCGACAGGGGGGGGGGGGCCGGAAGACgg                                                                                                     >  2:200119/1‑150 (MQ=255)
                                                      ccACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                  <  2:1803499/150‑1 (MQ=255)
                                                        aCCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATGACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAgg                                                                                                <  1:653233/150‑1 (MQ=255)
                                                         ccGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGc                                                                                               <  1:1178921/150‑1 (MQ=255)
                                                           ggTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCc                                                                                             >  2:1926721‑M2/1‑149 (MQ=255)
                                                             tGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcct                                                                                           >  2:1087762‑M2/1‑147 (MQ=255)
                                                             tGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcct                                                                                           >  2:2996259‑M2/1‑147 (MQ=255)
                                                             ggcGATCCTCAACCTGGGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcct                                                                                           >  2:4430589‑M2/2‑147 (MQ=255)
                                                               cGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgt                                                                                         >  2:2072510‑M2/1‑145 (MQ=255)
                                                                 aTCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtag                                                                                       <  1:1975436‑M1/150‑36 (MQ=255)
                                                                 aTCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtag                                                                                       <  1:2337236‑M1/150‑36 (MQ=255)
                                                                   ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtagga                                                                                     <  1:2345066‑M1/150‑38 (MQ=255)
                                                                   ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtagga                                                                                     <  1:2345091‑M1/150‑38 (MQ=255)
                                                                    cTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTGAACCCCCACCAACAGGggccgtgcccgcataagggcagcacccaccaggaggat                                                                                    >  1:2606957‑M1/1‑112 (MQ=255)
                                                                     tCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTGAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtaggagc                                                                                   <  1:4514243‑M1/150‑40 (MQ=255)
                                                                     tCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagc                                                                                   >  2:3761945‑M2/1‑139 (MQ=255)
                                                                     tCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagc                                                                                   >  2:1875541‑M2/1‑139 (MQ=255)
                                                                       aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcgg                                                                                 >  1:3216348‑M2/1‑137 (MQ=255)
                                                                       aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtaggagcgg                                                                                 <  2:3663496‑M1/150‑42 (MQ=255)
                                                                       aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtaggagcgg                                                                                 <  1:3374378‑M1/150‑42 (MQ=255)
                                                                       aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtaggagcgg                                                                                 <  1:3379844‑M1/150‑42 (MQ=255)
                                                                       aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGCTGCCTGCTTCGCGGGTGAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtaggagcgg                                                                                 <  2:2580365‑M1/150‑42 (MQ=255)
                                                                       aaCCGGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcgg                                                                                 >  1:1508961‑M2/1‑137 (MQ=255)
                                                                        aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtaggagcggg                                                                                <  2:4089389‑M1/150‑43 (MQ=255)
                                                                        aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcggg                                                                                >  1:741888‑M2/1‑136 (MQ=255)
                                                                        aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcggg                                                                                >  1:4421705‑M2/1‑136 (MQ=255)
                                                                        aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcggg                                                                                >  1:2196329‑M2/1‑136 (MQ=255)
                                                                        aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcggg                                                                                >  1:3569623‑M2/1‑136 (MQ=255)
                                                                        aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcggg                                                                                >  2:4253360‑M2/1‑136 (MQ=255)
                                                                           tgtgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTc                                                                             <  1:2173318/150‑1 (MQ=255)
                                                                             tgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAc                                                                           <  1:2955000/150‑1 (MQ=255)
                                                                               gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTGAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtagga                                                                                     >  2:3038951‑M1/1‑101 (MQ=255)
                                                                               gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                         <  2:3714637/150‑1 (MQ=255)
                                                                               gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                         <  1:81351/150‑1 (MQ=255)
                                                                                cTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGggccttgccgtcttcagggcagcacccacctgtaggagcgg                                                                                 >  1:1348280‑M1/1‑100 (MQ=255)
                                                                                 tggtggGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcgggtttacccgc                                                                       >  2:2807635‑M2/1‑127 (MQ=255)
                                                                                 tGGTTGGCGGGGGGCTGGCGATTACCCGGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                       <  1:226271/150‑1 (MQ=255)
                                                                                 tGGTTGGCGGGGGACTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                       <  2:3973157/150‑1 (MQ=255)
                                                                                 tGGTTGGCGGGGGACTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                       <  2:2673565/150‑1 (MQ=255)
                                                                                    ttGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGaa                                                                    <  2:755858/150‑1 (MQ=255)
                                                                                       gCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgca                                                                 <  1:1787647/150‑1 (MQ=255)
                                                                                               cTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTc                                                         <  2:3003243/150‑1 (MQ=255)
                                                                                               cTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTc                                                         <  2:4129881/150‑1 (MQ=255)
                                                                                                    gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtct                                                    <  1:1011391/150‑1 (MQ=255)
                                                                                                    gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtct                                                    <  1:1011399/150‑1 (MQ=255)
                                                                                                      ttACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCa                                                  <  1:3360216/150‑1 (MQ=255)
                                                                                                       tACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCaa                                                 <  2:307598/150‑1 (MQ=255)
                                                                                                          ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCgg                                                                                                                                             <  2:3425002/55‑1 (MQ=255)
                                                                                                           cTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCGGCTTCGCCGGTAAAACCCCCACCACAGGTGGGTGGGGCCCGGAAGAACGCAAGCCCCCCGGAGGAGCGGGGTGACCCCCCCAGCACGCACCGCGGGCGCAAAGtt                                            >  2:681805/1‑150 (MQ=255)
                                                                                                              ggggTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGttttt                                         <  2:3712799/150‑1 (MQ=255)
                                                                                                                     aGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGcc                                  <  1:4081163/150‑1 (MQ=255)
                                                                                                                     aGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGcc                                  <  2:1979396/150‑1 (MQ=255)
                                                                                                                     aGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGcc                                  <  2:1361127/150‑1 (MQ=255)
                                                                                                                     aGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGcc                                  <  1:4081214/150‑1 (MQ=255)
                                                                                                                      gAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCGGCGTACCCGGTAAAAACCCACCAAAACGGGGGGGCTGCCCCGAAGAACGCCCGGCCCCCGGAGGGGCGGGGTCACCCCCCCAGCACGCCCCACGGGCTCAACGGtttttaaaaccc                                  >  2:3467921/1‑143 (MQ=255)
                                                                                                                          cTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTTAAGACGGCAAGGCCCCTGTAGGAGGGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCt                             <  1:3528579/150‑1 (MQ=255)
                                                                                                                                     aCTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTcc                  <  2:3130445/150‑1 (MQ=255)
                                                                                                                                      cTGGAGACCGTGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACCGGTGGGTGCTGCCCTGAAGACGGCACGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCca                 <  1:1633028/150‑1 (MQ=255)
                                                                                                                                              cGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACTGGTGGGTGCGGCCCTGAAGACGGCACGGCCCCTGTAGGAGCGGGTTCACCCGCGACGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc         <  1:200119/150‑1 (MQ=255)
                                                                                                                                              cGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc         <  1:285445/150‑1 (MQ=255)
                                                                                                                                              cGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc         <  2:2907090/150‑1 (MQ=255)
                                                                                                                                               gAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg        <  1:1477546/150‑1 (MQ=255)
                                                                                                                                                aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCATGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCgg       <  1:1779282/150‑1 (MQ=255)
                                                                                                                                                    ggCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCACACGTTGCCGTCCACATTTCTCGGTAtg   >  1:1418049/1‑150 (MQ=255)
                                                                                                                                                     gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCCGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg  >  2:1169201/1‑150 (MQ=255)
                                                                                                                                                     gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccc                                 >  1:1866133/1‑119 (MQ=255)
                                                                                                                                                     gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccc                                 <  2:1866133/119‑1 (MQ=255)
                                                                                                                                                     gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg  <  2:1145692/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGG  >  NC_002947/4160885‑4161183

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