Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,597,173 T→G 100% L666R (CTG→CGG)  hyfB → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,597,1730TG97.7% 853.5 / 1.3 214L666R (CTG→CGG) hyfBhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (99/110);  total (102/112)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 219 total aligned reads displayed.

GGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGC  >  CP009273/2597027‑2597317
                                                                                                                                                  |                                                                                                                                                
ggCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGtgctg                                                                                                                                                 >  2:1197525/1‑148 (MQ=255)
  cTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTgcggct                                                                                                                                               >  2:789492/1‑148 (MQ=255)
     agagCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAAt                                                                                                                                            >  1:766067/1‑148 (MQ=255)
         cacaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCt                                                                                                                                        >  2:2506932/1‑148 (MQ=255)
          acaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTa                                                                                                                                       <  2:3351340/148‑1 (MQ=255)
          acaGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTa                                                                                                                                       <  1:1782985/148‑1 (MQ=255)
          acaGAACCCTTCTGGGATGAGCGGGCGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTa                                                                                                                                       <  1:2425221/148‑1 (MQ=255)
           caGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTAt                                                                                                                                      <  2:799146/148‑1 (MQ=255)
           caGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTAt                                                                                                                                      <  1:256707/148‑1 (MQ=255)
           caGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTAt                                                                                                                                      <  2:2863554/148‑1 (MQ=255)
             gAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTg                                                                                                                                    <  1:1494846/148‑1 (MQ=255)
              aaCCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGc                                                                                                                                   >  1:2442858/1‑148 (MQ=255)
               aCCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTCATCGCTATTGcc                                                                                                                                  >  2:1863156/1‑148 (MQ=255)
                cccTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGACg                                                                                                                                 >  1:2638261/1‑148 (MQ=255)
                 ccTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGt                                                                                                                                <  1:3353653/148‑1 (MQ=255)
                 ccTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGt                                                                                                                                >  1:1290284/1‑148 (MQ=255)
                 ccTTCTGGGAAGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGt                                                                                                                                <  1:1666817/148‑1 (MQ=255)
                           tGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAAt                                                                                                                      <  1:3400688/148‑1 (MQ=255)
                            gAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATc                                                                                                                     <  2:184812/148‑1 (MQ=255)
                             aGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCa                                                                                                                    >  1:908481/1‑148 (MQ=255)
                                   tGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAACTCACCATga                                                                                                              >  2:1712108/1‑148 (MQ=255)
                                   tGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATga                                                                                                              <  1:2777736/148‑1 (MQ=255)
                                   tGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATga                                                                                                              >  2:3476142/1‑148 (MQ=255)
                                     aTCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATgaga                                                                                                            <  1:145777/148‑1 (MQ=255)
                                      tCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGAc                                                                                                           <  2:1594621/148‑1 (MQ=255)
                                      tCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGAc                                                                                                           >  1:1947520/1‑148 (MQ=255)
                                        cGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACaa                                                                                                         <  1:995695/148‑1 (MQ=255)
                                         gCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACaaa                                                                                                        <  2:2336350/148‑1 (MQ=255)
                                         gCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACaaa                                                                                                        <  2:1891273/148‑1 (MQ=255)
                                          ccccATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAAc                                                                                                       >  1:707949/1‑148 (MQ=255)
                                            ccATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAActc                                                                                                     >  2:1980279/1‑148 (MQ=255)
                                            ccATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAActc                                                                                                     >  2:910918/1‑148 (MQ=255)
                                             cATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAActct                                                                                                    <  2:2734073/148‑1 (MQ=255)
                                             cATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAActct                                                                                                    >  1:554901/1‑148 (MQ=255)
                                              aTCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCtt                                                                                                   <  1:2201687/148‑1 (MQ=255)
                                              aTCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCtt                                                                                                   <  2:2413658/148‑1 (MQ=255)
                                                cGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTg                                                                                                 >  1:1995606/1‑148 (MQ=255)
                                                  tGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCg                                                                                               <  2:3427622/148‑1 (MQ=255)
                                                   gAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGa                                                                                              >  2:2664317/1‑148 (MQ=255)
                                                    aGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGAc                                                                                             <  1:327936/148‑1 (MQ=255)
                                                    aGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGAc                                                                                             <  1:2878126/148‑1 (MQ=255)
                                                     gcgcCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACg                                                                                            >  2:408880/1‑148 (MQ=255)
                                                     gcgcCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACg                                                                                            >  2:3185087/1‑148 (MQ=255)
                                                       gcCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGa                                                                                          <  1:1724869/148‑1 (MQ=255)
                                                        ccACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGat                                                                                         >  1:567622/1‑148 (MQ=255)
                                                           cccAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATc                                                                                      <  2:503666/148‑1 (MQ=255)
                                                            ccAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCt                                                                                     <  1:677767/148‑1 (MQ=255)
                                                               aCGGCTGGCCAAAGACATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGt                                                                                  <  2:2201578/148‑1 (MQ=255)
                                                               aCGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGt                                                                                  <  2:1816522/148‑1 (MQ=255)
                                                                cGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTc                                                                                 >  2:869398/1‑148 (MQ=255)
                                                                cGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTc                                                                                 >  2:1090391/1‑148 (MQ=255)
                                                                  gCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAt                                                                               <  1:93217/148‑1 (MQ=255)
                                                                   cTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAtt                                                                              <  1:408420/148‑1 (MQ=255)
                                                                   cTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAtt                                                                              <  1:3549142/148‑1 (MQ=255)
                                                                     ggCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGACTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAtttt                                                                            <  2:2106169/148‑1 (MQ=255)
                                                                      gCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCAttttt                                                                           <  2:1468146/148‑1 (MQ=255)
                                                                           aGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGtt                                                                      >  1:1060233/1‑148 (MQ=255)
                                                                           aGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGtt                                                                      >  1:882726/1‑148 (MQ=255)
                                                                              aaTACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGc                                                                   <  2:3148704/148‑1 (MQ=255)
                                                                                tACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGcac                                                                 >  2:662471/1‑148 (MQ=255)
                                                                                    gCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGc                                                             <  1:112495/148‑1 (MQ=255)
                                                                                    gCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGc                                                             <  1:2104592/148‑1 (MQ=255)
                                                                                     cATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGcc                                                            >  2:2790783/1‑148 (MQ=255)
                                                                                      aTCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCg                                                           <  1:2418706/148‑1 (MQ=255)
                                                                                          gCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGAt                                                       <  2:139098/148‑1 (MQ=255)
                                                                                           cAAAGCGGCGACTTTCGTCTCTATTGCCTGTAAGTGGTAGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGaat                                                      >  2:2954277/1‑146 (MQ=255)
                                                                                            aaaGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTc                                                     <  2:2306089/148‑1 (MQ=255)
                                                                                              aGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTctg                                                   >  1:3260916/1‑148 (MQ=255)
                                                                                               gcggcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTctgc                                                  <  2:2529074/148‑1 (MQ=255)
                                                                                               gcggcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTctgc                                                  <  2:2047515/148‑1 (MQ=255)
                                                                                               gcggcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTctgc                                                  >  1:3464461/1‑148 (MQ=255)
                                                                                                  gcgACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGa                                               <  2:3066941/148‑1 (MQ=255)
                                                                                                    gACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATg                                             <  1:2506932/148‑1 (MQ=255)
                                                                                                     aCTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGc                                            <  2:832935/148‑1 (MQ=255)
                                                                                                        ttCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTaa                                         >  1:3524519/1‑148 (MQ=255)
                                                                                                        ttCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATAGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATATGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTaa                                         >  1:1339698/1‑148 (MQ=255)
                                                                                                            tctcTATTGCCTGTATGTGGTCGCCTCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAAcccc                                     >  1:1269777/1‑148 (MQ=255)
                                                                                                            tctcTATTGCCTGTATGTGGTCGCCTCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAAcccc                                     >  1:3048758/1‑148 (MQ=255)
                                                                                                             ctctATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCa                                    <  2:1947520/148‑1 (MQ=255)
                                                                                                                 aTTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACttt                                >  1:2722624/1‑148 (MQ=255)
                                                                                                                 aTTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACttt                                <  1:2786004/148‑1 (MQ=255)
                                                                                                                 aTTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACttt                                <  1:2854770/148‑1 (MQ=255)
                                                                                                                 aTTGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCGATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGAGTCTGCTGATGCTAACCCCACttt                                >  1:3103780/1‑148 (MQ=255)
                                                                                                                  ttGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACtttt                               <  2:3518103/148‑1 (MQ=255)
                                                                                                                   tGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACttttt                              >  2:1354898/1‑148 (MQ=255)
                                                                                                                   tGCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACttttt                              >  2:3129027/1‑148 (MQ=255)
                                                                                                                    gCCTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTa                             <  1:2653198/148‑1 (MQ=255)
                                                                                                                     ccTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTAc                            >  2:2252802/1‑148 (MQ=255)
                                                                                                                      cTGTATGTGGTCGCCGCACTGGTTGTGCGGCTAAACGCTATTGCCGTCTAAGGAAATCAACATGATACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTATTCTGATGCTAACCCCACATTTTACg                           >  1:1596368/1‑148 (MQ=255)
                                                                                                                        gTATGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACggg                         <  2:1507442/148‑1 (MQ=255)
                                                                                                                          aTGTGGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTa                       >  1:2204289/1‑148 (MQ=255)
                                                                                                                           tgtgGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTAt                      <  1:1516839/148‑1 (MQ=255)
                                                                                                                             tgGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTCttt                    >  2:1643632/1‑148 (MQ=255)
                                                                                                                             tgGTCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTAttt                    <  1:1532856/148‑1 (MQ=255)
                                                                                                                                tCGCCGCACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTccc                 >  1:186625/1‑148 (MQ=255)
                                                                                                                                  gccgcACTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCgg               <  2:1246329/148‑1 (MQ=255)
                                                                                                                                       aCTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGAt          <  2:707949/148‑1 (MQ=255)
                                                                                                                                       aCTGGTTGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGAt          <  1:1020456/148‑1 (MQ=255)
                                                                                                                                            ttGTGCGGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATGTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTCTTTCCCGGCAGATAcgcg     >  2:2160925/1‑148 (MQ=255)
                                                                                                                                               tgcggcTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATAcgcgcgc  >  2:329413/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGCTCAGAGCACAGAACCCTTCTGGGATGAGCGGGTGATCCGCCCCATCGTGAGCGCCACCCAACGGCTGGCCAAAGAAATACAGCATCTGCAAAGCGGCGACTTTCGTCTCTATTGCCTGTATGTGGTCGCCGCACTGGTTGTGCTGCTAATCGCTATTGCCGTCTAAGGAAATCACCATGAGACAAACTCTTTGCGACGGATATCTGGTCATTTTTGCGTTAGCACAGGCCGTGATTCTGCTGATGCTAACCCCACTTTTTACGGGTATTTCCCGGCAGATACGCGCGC  >  CP009273/2597027‑2597317

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: