Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 591,082 G→A 5.8% R109H (CGC→CAC)  PP_0501 → NAD‑dependent epimerase/dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947591,0820GA 5.8% 421.1 / 10.4 138R109H (CGC→CAC) PP_0501NAD‑dependent epimerase/dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (66/64);  new base A (4/4);  total (70/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 140 total aligned reads displayed.

CGATGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGAC  >  NC_002947/590933‑591230
                                                                                                                                                     |                                                                                                                                                    
   tGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGt                                                                                                                                                   <  1:231936/150‑1 (MQ=255)
       gggCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCCCCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCACGCGCAt                                                                                                                                               <  2:1762858/150‑1 (MQ=255)
        ggCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATg                                                                                                                                              >  1:739554/1‑150 (MQ=255)
         gCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATgg                                                                                                                                             >  2:1992992/1‑150 (MQ=255)
           tGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATg                                                                                                                                              <  1:2008418/147‑1 (MQ=255)
           tGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATg                                                                                                                                              >  2:2008418/1‑147 (MQ=255)
             gTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTg                                                                                                                                         >  2:1206915/1‑150 (MQ=255)
              tCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTgc                                                                                                                                        <  2:284808/150‑1 (MQ=255)
                  tCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCg                                                                                                                                    >  2:1606149/1‑150 (MQ=255)
                      gCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTg                                                                                                                                >  2:1854239/1‑150 (MQ=255)
                      gCTACGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGACTCCGTTCTGGCATCCGTCGAAGACCCAGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTg                                                                                                                                <  1:1967994/150‑1 (MQ=255)
                         gCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTATCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGCAGTTGGCCGTGTGctc                                                                                                                             <  2:711963/150‑1 (MQ=255)
                                    ccGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcag                                                                                                                  >  2:1308623/1‑150 (MQ=255)
                                    ccGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcag                                                                                                                  <  2:1032718/150‑1 (MQ=255)
                                     cGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagc                                                                                                                 >  1:1453901/1‑150 (MQ=255)
                                     cGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagc                                                                                                                 <  2:725802/150‑1 (MQ=255)
                                       tGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTgcgc                                                                                                                                      >  1:2059674/1‑127 (MQ=255)
                                             ggtgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgt                                                                                                         <  1:1095571/150‑1 (MQ=255)
                                               tgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgtgt                                                                                                       >  1:2257129/1‑150 (MQ=255)
                                                   tCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACgg                                                                                                   <  1:741152/150‑1 (MQ=255)
                                                           cggTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                           <  2:2368088/150‑1 (MQ=255)
                                                           cggTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                           <  1:1992992/150‑1 (MQ=255)
                                                            ggTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATgg                                                                                                                                             <  2:275178/99‑1 (MQ=255)
                                                             gTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGt                                                                                         >  1:2264787/1‑150 (MQ=255)
                                                             gTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGt                                                                                         <  1:1308623/150‑1 (MQ=255)
                                                              tAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAAt                                                                                            <  1:77065/146‑1 (MQ=255)
                                                              tAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAAt                                                                                            >  2:77065/1‑146 (MQ=255)
                                                               aGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGTGTGCTCTTCTCCTCTAGCGCAGCTGTGTACGGCAACAATTGTGa                                                                                       <  2:499355/150‑1 (MQ=255)
                                                                gCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGaa                                                                                      >  2:1360865/1‑150 (MQ=255)
                                                                acGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGaa                                                                                      <  2:808674/149‑1 (MQ=255)
                                                                 cGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAg                                                                                     <  1:2113357/150‑1 (MQ=255)
                                                                  gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagcgcagcg                                                                                                           <  1:559863/127‑1 (MQ=255)
                                                                  gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCcagcgcagcg                                                                                                           >  2:559863/1‑127 (MQ=255)
                                                                  gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAgg                                                                                    <  1:1540168/150‑1 (MQ=255)
                                                                   tCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAggg                                                                                   <  2:1973463/150‑1 (MQ=255)
                                                                       ttCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAAt                                                                               <  2:1793884/150‑1 (MQ=255)
                                                                            gCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTc                                                                                                                          <  2:1935872/102‑1 (MQ=255)
                                                                            gCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTc                                                                                                                          >  1:1935872/1‑102 (MQ=255)
                                                                              aTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTc                                                                        <  1:996946/150‑1 (MQ=255)
                                                                              aTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTc                                                                        >  1:830953/1‑150 (MQ=255)
                                                                              aTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTc                                                                        >  2:919433/1‑150 (MQ=255)
                                                                               tCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTcc                                                                       <  2:2257129/150‑1 (MQ=255)
                                                                               tCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGAGTACATGGCGTGCGCCGGGTGCTCTTCGCCTCCGGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGAtttac                                                                       >  1:1364831/1‑148 (MQ=255)
                                                                                  gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                                    <  1:1051896/150‑1 (MQ=255)
                                                                                  gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCACGCGCAt                                                                                                                                               <  1:448175/75‑1 (MQ=255)
                                                                                  gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCACGCGCAt                                                                                                                                               >  2:448175/1‑75 (MQ=255)
                                                                                   tCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                                    <  1:273979/149‑1 (MQ=255)
                                                                                     gAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGac                                                                 >  1:460053/1‑150 (MQ=255)
                                                                                     gAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGac                                                                 >  1:1802705/1‑150 (MQ=255)
                                                                                     gAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGac                                                                 <  1:248496/150‑1 (MQ=255)
                                                                                      aaGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAAAAATGGTGAAGGGGAATCGATTTCCGAAGaca                                                                >  1:2254584/1‑150 (MQ=255)
                                                                                      aaGACCCGGTGCGTACCCACCAGAGTAACTTAATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGaca                                                                >  1:1170493/1‑150 (MQ=255)
                                                                                           ccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCa                                                           >  1:1680943/1‑150 (MQ=255)
                                                                                                   tACCCACCAGAGTAAGTGTATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCAGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGcc                                                   <  2:1061404/150‑1 (MQ=255)
                                                                                                   tACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgt                                                                                                         <  2:1192771/96‑1 (MQ=255)
                                                                                                   tACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGgt                                                                                                         >  1:1192771/1‑96 (MQ=255)
                                                                                                      ccaccaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGac                                                                 >  1:103050/1‑133 (MQ=255)
                                                                                                        accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCAt                                                                                                                                               >  2:653528/1‑53 (MQ=255)
                                                                                                        accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGa                                              >  2:1128580/1‑150 (MQ=255)
                                                                                                        accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGa                                              >  1:838632/1‑150 (MQ=255)
                                                                                                        accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGa                                              >  1:1713117/1‑150 (MQ=255)
                                                                                                           agagTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGAccc                                           <  2:1453901/150‑1 (MQ=255)
                                                                                                              gTACCTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGGGCGCCGGGTGCTCTTCGCCCCCAGCGCAGCGGTGTGCGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCtt                                        >  2:400159/1‑150 (MQ=255)
                                                                                                               tAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGa                                                                            <  1:1213540/113‑1 (MQ=255)
                                                                                                                aaCTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAAt                                                                               >  1:1615159/1‑109 (MQ=255)
                                                                                                                aaCTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAAt                                                                               <  2:1615159/109‑1 (MQ=255)
                                                                                                                   ttCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGgc                                                      >  2:580767/1‑131 (MQ=255)
                                                                                                                   ttCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGgc                                                      <  1:580767/131‑1 (MQ=255)
                                                                                                                   ttCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGcc                                   >  2:361663/1‑150 (MQ=255)
                                                                                                                    tCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCg                                  >  2:2121421/1‑150 (MQ=255)
                                                                                                                      aTCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                <  2:1170493/150‑1 (MQ=255)
                                                                                                                      aTCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                <  1:361663/150‑1 (MQ=255)
                                                                                                                          gTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGt                                                                                         >  2:1531281/1‑89 (MQ=255)
                                                                                                                          gTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGt                                                                                         <  1:1531281/89‑1 (MQ=255)
                                                                                                                          gTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAggg                                                                                   >  2:1778632/1‑95 (MQ=255)
                                                                                                                          gTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAggg                                                                                   <  1:1778632/95‑1 (MQ=255)
                                                                                                                            aCCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAg                          <  2:2041554/150‑1 (MQ=255)
                                                                                                                               cTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACg                                                              <  2:1521384/111‑1 (MQ=255)
                                                                                                                               cTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTg                       >  1:2190394/1‑150 (MQ=255)
                                                                                                                                tCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCca                                                                                                                   >  1:1098533/1‑57 (MQ=255)
                                                                                                                                tCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCca                                                                                                                   <  2:1098533/57‑1 (MQ=255)
                                                                                                                                  aaCGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAg                                                        <  1:274884/114‑1 (MQ=255)
                                                                                                                                  aaCGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAg                                                        >  2:274884/1‑114 (MQ=255)
                                                                                                                                     gTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGt                 <  2:1802705/150‑1 (MQ=255)
                                                                                                                                     gTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGt                 <  1:1074086/150‑1 (MQ=255)
                                                                                                                                       ttGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGa               <  2:2190394/150‑1 (MQ=255)
                                                                                                                                       ttGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGa               <  2:1680943/150‑1 (MQ=255)
                                                                                                                                         gCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGc             >  2:249470/1‑150 (MQ=255)
                                                                                                                                          gcacgCCATGCACGCGCATGGCTTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCa            >  1:2211765/5‑150 (MQ=255)
                                                                                                                                          cGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCa            <  2:838632/150‑1 (MQ=255)
                                                                                                                                           gAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACaa                           >  2:972108/1‑134 (MQ=255)
                                                                                                                                           gAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACaa                           <  1:972108/134‑1 (MQ=255)
                                                                                                                                           gAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAg           >  2:100625/1‑150 (MQ=255)
                                                                                                                                            aaGCCATGCGCGTGCATGGCTTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCa                   >  1:279969/1‑141 (MQ=255)
                                                                                                                                              gCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTAc        <  1:2121421/150‑1 (MQ=255)
                                                                                                                                              gCCATGCACGCGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTAc        >  1:1159326/1‑150 (MQ=255)
                                                                                                                                                cATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCt      >  1:233971/1‑150 (MQ=255)
                                                                                                                                                cATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCt      <  1:1786251/150‑1 (MQ=255)
                                                                                                                                                  tGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTgg    >  1:54230/1‑150 (MQ=255)
                                                                                                                                                    cgcgTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGAc  <  1:1461289/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGATGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGCCGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGAC  >  NC_002947/590933‑591230

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: