Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,848,892 C→A 17.2% A399E (GCG→GAG)  fhlA → formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,848,8920CA17.2% 506.6 / 62.7 174A399E (GCG→GAG) fhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (78/66);  new base A (16/14);  total (94/80)
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.97e-01

Only 100 of 181 total aligned reads displayed.

ATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGATTGCTGGAAAG  >  CP009273/2848745‑2849039
                                                                                                                                                   |                                                                                                                                                   
aTCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGACTTTGGCGAGATTCTTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGgc                                                                                                                                                     >  1:1861928/1‑148 (MQ=255)
  cAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGgcgc                                                                                                                                                   >  1:1504718/1‑148 (MQ=255)
  aaGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCACTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGTCCGCAGCGAAGACATGTACAGCGTGCTTAAACAAGTGGAAATGGTGgcgc                                                                                                                                                   >  1:3144656/2‑148 (MQ=255)
         tCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGa                                                                                                                                            >  1:2290813/1‑148 (MQ=255)
           cATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCAATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACa                                                                                                                                          <  1:788838/148‑1 (MQ=255)
             tCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGt                                                                                                                                        <  2:3201386/148‑1 (MQ=255)
             tCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGt                                                                                                                                        >  1:3240528/1‑148 (MQ=255)
              cGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTa                                                                                                                                       >  1:2314127/1‑148 (MQ=255)
              cGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTa                                                                                                                                       <  2:2668028/148‑1 (MQ=255)
                tCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTAcc                                                                                                                                     >  1:2245865/1‑148 (MQ=255)
                 cTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCg                                                                                                                                    <  2:1785778/148‑1 (MQ=255)
                 cTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCg                                                                                                                                    >  2:2804226/1‑148 (MQ=255)
                      aGATCGGCTGTTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTg                                                                                                                               <  2:240655/148‑1 (MQ=255)
                        aaCGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGAt                                                                                                                             >  1:3101208/1‑148 (MQ=255)
                        aaCGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGAt                                                                                                                             >  2:2623690/1‑148 (MQ=255)
                        aaCGGCTGGTGGAGGAAAACCTCGCCCTGCCCGAGCAGCTCAACAAGGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGAt                                                                                                                             <  2:2654124/148‑1 (MQ=255)
                          cGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATcc                                                                                                                           >  1:1167159/1‑148 (MQ=255)
                           ggCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCt                                                                                                                          <  2:2319982/148‑1 (MQ=255)
                            gCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTc                                                                                                                         <  1:1258180/148‑1 (MQ=255)
                            gCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTc                                                                                                                         >  1:2747136/1‑148 (MQ=255)
                             cTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCg                                                                                                                        <  2:3229561/148‑1 (MQ=255)
                             cTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCg                                                                                                                        <  2:2972093/148‑1 (MQ=255)
                              tGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCgg                                                                                                                       <  2:2013145/148‑1 (MQ=255)
                               ggTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGt                                                                                                                      >  1:1941934/1‑148 (MQ=255)
                                  tgatgaAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGaa                                                                                                                   >  2:2463641/1‑148 (MQ=255)
                                    atgaAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAAc                                                                                                                 >  2:882052/1‑148 (MQ=255)
                                    atgaAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAAc                                                                                                                 >  2:882042/1‑148 (MQ=255)
                                    atgaAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAAc                                                                                                                 >  2:3054288/1‑148 (MQ=255)
                                     tgaAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACt                                                                                                                >  2:2514156/1‑148 (MQ=255)
                                       aaaaCCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTgg                                                                                                              >  1:48363/1‑148 (MQ=255)
                                             tCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACggg                                                                                                        <  2:1000411/148‑1 (MQ=255)
                                             tCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACggg                                                                                                        <  2:793766/148‑1 (MQ=255)
                                             tCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACggg                                                                                                        >  1:190081/1‑148 (MQ=255)
                                                     aCCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGc                                                                                                <  2:2498594/148‑1 (MQ=255)
                                                       cGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTg                                                                                              <  2:2506625/148‑1 (MQ=255)
                                                         agcagcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGAt                                                                                            >  2:282808/1‑148 (MQ=255)
                                                         agcagcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGAt                                                                                            <  1:1674041/148‑1 (MQ=255)
                                                         agcagcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGAt                                                                                            <  1:1670723/148‑1 (MQ=255)
                                                          gcagcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGAtt                                                                                           >  2:359036/1‑148 (MQ=255)
                                                          gcagcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGAtt                                                                                           >  1:2059027/1‑148 (MQ=255)
                                                           cagcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTg                                                                                          >  2:1551810/1‑148 (MQ=255)
                                                           cagcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTg                                                                                          <  1:2466311/148‑1 (MQ=255)
                                                           cagcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTg                                                                                          <  1:2466317/148‑1 (MQ=255)
                                                            agcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGc                                                                                         <  1:1950508/148‑1 (MQ=255)
                                                             gcTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGcc                                                                                        >  1:115968/1‑148 (MQ=255)
                                                                caacaaTGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTCAAGAGCTGATTGCCCGt                                                                                     >  2:3249071/1‑148 (MQ=255)
                                                                caacaaTGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGt                                                                                     >  1:1038937/1‑148 (MQ=255)
                                                                caacaaTGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGt                                                                                     >  1:2227910/1‑148 (MQ=255)
                                                                    aaTGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGa                                                                                 <  2:1465988/148‑1 (MQ=255)
                                                                    aaTGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGa                                                                                 <  1:683682/148‑1 (MQ=255)
                                                                      tGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATc                                                                               >  1:2610183/1‑148 (MQ=255)
                                                                       gTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATcc                                                                              <  1:2018811/148‑1 (MQ=255)
                                                                        ttGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCa                                                                             <  1:2865617/148‑1 (MQ=255)
                                                                           aTAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATaa                                                                          >  2:3299322/1‑148 (MQ=255)
                                                                            tAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAAt                                                                         >  2:1313794/1‑148 (MQ=255)
                                                                             aGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAAtc                                                                        >  1:1895848/1‑148 (MQ=255)
                                                                               tGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAAtctc                                                                      <  2:3260825/148‑1 (MQ=255)
                                                                                 aaTTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAg                                                                    >  2:3217499/1‑148 (MQ=255)
                                                                                 aaTTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAg                                                                    >  2:2937013/1‑148 (MQ=255)
                                                                                 aaTTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAg                                                                    <  1:2688385/148‑1 (MQ=255)
                                                                                  aTTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGt                                                                   <  2:696186/148‑1 (MQ=255)
                                                                                  aTTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGt                                                                   <  1:2632547/148‑1 (MQ=255)
                                                                                   tttGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTg                                                                  >  2:3137056/1‑148 (MQ=255)
                                                                                    ttGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTgg                                                                 <  1:897073/148‑1 (MQ=255)
                                                                                    ttGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTgg                                                                 <  1:897027/148‑1 (MQ=255)
                                                                                     tGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTggg                                                                <  2:2695176/148‑1 (MQ=255)
                                                                                       gCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCg                                                              <  1:882042/148‑1 (MQ=255)
                                                                                       gCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCg                                                              <  1:882052/148‑1 (MQ=255)
                                                                                        cGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGt                                                             <  2:1218611/148‑1 (MQ=255)
                                                                                            attattGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGtaata                                                         >  1:1451962/1‑148 (MQ=255)
                                                                                            attattGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGtaata                                                         >  2:718306/1‑148 (MQ=255)
                                                                                            attattGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGCGCGATCCATAATCTCAGTGGGCGtaata                                                         >  2:3310879/1‑148 (MQ=255)
                                                                                             ttattGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGtaataa                                                        >  1:2549130/1‑148 (MQ=255)
                                                                                             ttattGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGtaataa                                                        <  1:2228887/148‑1 (MQ=255)
                                                                                              tattGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTaataat                                                       <  1:1860396/148‑1 (MQ=255)
                                                                                              tattGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTaataat                                                       <  2:1324071/148‑1 (MQ=255)
                                                                                               attGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATc                                                      >  2:18414/1‑148 (MQ=255)
                                                                                                ttGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATcg                                                     <  1:461042/148‑1 (MQ=255)
                                                                                                ttGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATcg                                                     >  2:1020254/1‑148 (MQ=255)
                                                                                                 tGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATcgc                                                    <  2:374368/148‑1 (MQ=255)
                                                                                                     cGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCa                                                <  2:217233/148‑1 (MQ=255)
                                                                                                         gCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGt                                            <  2:1600338/148‑1 (MQ=255)
                                                                                                         gCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGt                                            >  2:1375055/1‑148 (MQ=255)
                                                                                                             aGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCaaa                                        >  2:1974770/1‑148 (MQ=255)
                                                                                                                cATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATg                                     >  2:1288621/1‑148 (MQ=255)
                                                                                                                  tGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGaa                                   >  1:97484/1‑148 (MQ=255)
                                                                                                                  tGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGaa                                   <  1:2800054/148‑1 (MQ=255)
                                                                                                                     aCAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTg                                >  1:1184446/1‑148 (MQ=255)
                                                                                                                     aCAGCGTGCTTAAACAAGTTGAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTg                                >  1:1028742/1‑148 (MQ=255)
                                                                                                                          gTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCgcgg                           <  1:821012/148‑1 (MQ=255)
                                                                                                                             cTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGa                        >  1:1730263/1‑148 (MQ=255)
                                                                                                                                 aaCAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGAtgcc                    >  1:313105/1‑148 (MQ=255)
                                                                                                                                   cAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGAtgcctg                  >  2:2910280/1‑148 (MQ=255)
                                                                                                                                        tGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGa             >  2:255389/1‑148 (MQ=255)
                                                                                                                                        tGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGa             >  2:1313313/1‑148 (MQ=255)
                                                                                                                                         gAAATGGTGGAGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGAt            <  2:661286/148‑1 (MQ=255)
                                                                                                                                           aagggtggCGCACAGTGACAGTTCCGTGCTGATCCTCGGTGAAACTGGCACGGGTACAGAGCTGATTGCGCGTGCGATCCCTAATCTCCGTGGGCGGAATCCTCGCCGCCTGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGCTTg          >  2:1763993/4‑148 (MQ=255)
                                                                                                                                             tggtggCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGATTGCt        <  1:1560581/148‑1 (MQ=255)
                                                                                                                                                tggCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGATTGCTGGa     <  2:48363/148‑1 (MQ=255)
                                                                                                                                                   cgcAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGATTGCTGGAAAg  >  1:1268992/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCAGGAAATCCATCGTCTGAAAGAACGGCTGGTTGATGAAAACCTCGCCCTGACCGAGCAGCTCAACAATGTTGATAGTGAATTTGGCGAGATTATTGGCCGCAGCGAAGCCATGTACAGCGTGCTTAAACAAGTTGAAATGGTGGCGCAAAGTGACAGTACCGTGCTGATCCTCGGTGAAACTGGCACGGGTAAAGAGCTGATTGCCCGTGCGATCCATAATCTCAGTGGGCGTAATAATCGCCGCATGGTCAAAATGAACTGCGCGGCGATGCCTGCCGGATTGCTGGAAAG  >  CP009273/2848745‑2849039

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