Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,758 G→T 7.4% intergenic (‑342/‑28) yigE ← / → corA DUF2233 family protein/magnesium/nickel/cobalt transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,994,7580GT 7.4% 1134.8 / 40.7 310intergenic (‑342/‑28)yigE/corADUF2233 family protein/magnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base G (154/133);  new base T (12/11);  total (166/144)
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.59e-01

Only 100 of 316 total aligned reads displayed.

CCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGC  >  CP009273/3994610‑3994906
                                                                                                                                                    |                                                                                                                                                    
ccTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCg                                                                                                                                                      >  1:3875902/1‑149 (MQ=255)
 cTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGa                                                                                                                                                     <  2:3311271/149‑1 (MQ=255)
 cTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGa                                                                                                                                                     <  1:732649/149‑1 (MQ=255)
  tGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGat                                                                                                                                                    >  2:3570230/1‑149 (MQ=255)
  tGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGat                                                                                                                                                    <  1:2521795/149‑1 (MQ=255)
    tCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGatat                                                                                                                                                  <  2:2268186/149‑1 (MQ=255)
          aGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAg                                                                                                                                            >  2:3178189/1‑149 (MQ=255)
          aGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAg                                                                                                                                            <  1:319938/149‑1 (MQ=255)
          aGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAg                                                                                                                                            >  2:3657064/1‑149 (MQ=255)
          aGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACAGAACTGTCCTATATTTTAAg                                                                                                                                            >  2:4357807/1‑149 (MQ=255)
             aTCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCAt                                                                                                                                         >  2:3741726/1‑149 (MQ=255)
              tCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCAtt                                                                                                                                        >  1:1792730/1‑149 (MQ=255)
              tCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCAtt                                                                                                                                        <  2:4201390/149‑1 (MQ=255)
               cAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTg                                                                                                                                       <  2:1613006/149‑1 (MQ=255)
                 aCTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTggg                                                                                                                                     <  2:146357/149‑1 (MQ=255)
                 aCTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTggg                                                                                                                                     <  1:1010010/149‑1 (MQ=255)
                  cTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGa                                                                                                                                    <  2:4214179/149‑1 (MQ=255)
                   tCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAg                                                                                                                                   <  1:2681599/149‑1 (MQ=255)
                     cTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTc                                                                                                                                 <  2:4726883/149‑1 (MQ=255)
                        aGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCg                                                                                                                              >  2:2218418/1‑149 (MQ=255)
                         ggCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCgg                                                                                                                             >  2:3800260/1‑149 (MQ=255)
                               cccATTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTTTGTGCCCACCTAACTGTCAGATATTTTAAGCATTGGGAGTCCCGGTCATGc                                                                                                                       >  2:2430891/1‑149 (MQ=255)
                                  tttttGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGa                                                                                                                    <  1:4399200/149‑1 (MQ=255)
                                       gCTGATAGCCTGAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCTATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCa                                                                                                               <  1:1272729/149‑1 (MQ=255)
                                         tGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCAtt                                                                                                             <  2:1669991/149‑1 (MQ=255)
                                          gATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCAttt                                                                                                            >  2:2765240/1‑149 (MQ=255)
                                           aTAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTc                                                                                                           >  2:4334800/1‑149 (MQ=255)
                                            tAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCa                                                                                                          <  2:1212553/149‑1 (MQ=255)
                                            tAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCa                                                                                                          >  2:2701371/1‑149 (MQ=255)
                                                cTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTg                                                                                                      >  2:2239415/1‑149 (MQ=255)
                                                 ttAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTgg                                                                                                     >  1:2016670/1‑149 (MQ=255)
                                                    gCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGaaa                                                                                                  >  1:4242392/1‑149 (MQ=255)
                                                    gCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGaaa                                                                                                  >  2:3892804/1‑149 (MQ=255)
                                                    gCGGTTGTCAGCGACCTCAATTGTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGaaa                                                                                                  <  2:4624179/149‑1 (MQ=255)
                                                       gTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCCTTTCATCTGGaaactc                                                                                               >  2:2794574/1‑146 (MQ=255)
                                                        ttGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCTATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATaa                                                                                              >  2:3412018/1‑149 (MQ=255)
                                                         tGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAAc                                                                                             >  1:3281606/1‑149 (MQ=255)
                                                         tGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAAc                                                                                             <  1:2224094/149‑1 (MQ=255)
                                                            cAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGa                                                                                          <  2:789742/149‑1 (MQ=255)
                                                              gCGACCTCATTTTTTCCCGTCGCGCTGTGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACt                                                                                        <  1:4114490/149‑1 (MQ=255)
                                                                gACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGa                                                                                      >  2:1784605/1‑149 (MQ=255)
                                                                gACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGa                                                                                      >  1:4538254/1‑149 (MQ=255)
                                                                 aacTCAATTTTTCCCGTCGAGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAAATGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGAc                                                                                     <  2:4395724/147‑1 (MQ=255)
                                                                 aCCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGAc                                                                                     >  2:2956609/1‑149 (MQ=255)
                                                                   cTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGAccc                                                                                   <  2:2639627/149‑1 (MQ=255)
                                                                      aaTTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGc                                                                                >  2:4105739/1‑149 (MQ=255)
                                                                      aaTTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGc                                                                                >  1:1529733/1‑149 (MQ=255)
                                                                       aTTTGTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCt                                                                               <  2:2890896/149‑1 (MQ=255)
                                                                            tCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAg                                                                          >  1:4088015/1‑149 (MQ=255)
                                                                              ccGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTc                                                                        <  1:4732322/149‑1 (MQ=255)
                                                                              ccGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTc                                                                        <  1:3694683/149‑1 (MQ=255)
                                                                                   gcgcTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAaga                                                                   <  2:1343409/149‑1 (MQ=255)
                                                                                    cgcTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAagag                                                                  <  1:638657/149‑1 (MQ=255)
                                                                                      cTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTc                                                                <  2:2445695/149‑1 (MQ=255)
                                                                                       tGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTca                                                               >  1:3186055/1‑149 (MQ=255)
                                                                                       tGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTca                                                               <  2:4078584/149‑1 (MQ=255)
                                                                                       tGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTca                                                               <  2:4078674/149‑1 (MQ=255)
                                                                                        gAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTcac                                                              >  2:1445080/1‑149 (MQ=255)
                                                                                          gTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACaa                                                            >  2:1135341/1‑149 (MQ=255)
                                                                                           tCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAAc                                                           >  1:2904186/1‑149 (MQ=255)
                                                                                              ggCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCTATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAAcccc                                                        >  2:490888/1‑149 (MQ=255)
                                                                                              ggCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAAcccc                                                        >  1:496829/1‑149 (MQ=255)
                                                                                               gCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCt                                                       >  2:4097440/1‑149 (MQ=255)
                                                                                               gCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCt                                                       >  2:2554608/1‑149 (MQ=255)
                                                                                                 tGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCTATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTg                                                     >  2:206484/1‑149 (MQ=255)
                                                                                                 tGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTg                                                     <  1:1310000/149‑1 (MQ=255)
                                                                                                  gTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGt                                                    <  1:1257452/149‑1 (MQ=255)
                                                                                                       aTGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATg                                               <  2:3070064/149‑1 (MQ=255)
                                                                                                         ggTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCa                                             <  2:734689/149‑1 (MQ=255)
                                                                                                           tCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGt                                           >  2:190332/1‑149 (MQ=255)
                                                                                                           tCTGAAACCCAATTTGTTTCTGTGTGCCCACAGAACTGTCCTATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGt                                           <  1:4357608/149‑1 (MQ=255)
                                                                                                            cTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTa                                          <  1:1887281/149‑1 (MQ=255)
                                                                                                             tGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTAt                                         >  2:2475665/1‑149 (MQ=255)
                                                                                                              gAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATg                                        >  2:197005/1‑149 (MQ=255)
                                                                                                              gAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATg                                        <  1:376303/149‑1 (MQ=255)
                                                                                                               aaaCCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATgg                                       <  2:3240927/149‑1 (MQ=255)
                                                                                                               aaaCCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATgg                                       >  1:1456213/1‑149 (MQ=255)
                                                                                                                 aCCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGAt                                     >  1:854603/1‑149 (MQ=255)
                                                                                                                  cccAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGAtt                                    >  1:715644/1‑149 (MQ=255)
                                                                                                                  cccAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGAtt                                    >  1:2011064/1‑149 (MQ=255)
                                                                                                                    cAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGa                                  >  1:4062247/1‑149 (MQ=255)
                                                                                                                    cAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGa                                  <  2:570091/149‑1 (MQ=255)
                                                                                                                      aTTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATc                                <  2:4464732/149‑1 (MQ=255)
                                                                                                                       tttGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCt                               >  2:4499663/1‑149 (MQ=255)
                                                                                                                       tttGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCt                               >  1:4426428/1‑149 (MQ=255)
                                                                                                                        ttGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCtt                              <  2:1750410/149‑1 (MQ=255)
                                                                                                                          gTTTCTTTGTGACCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTTGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGt                            <  1:1767484/149‑1 (MQ=255)
                                                                                                                           tttCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTc                           >  1:1356032/1‑149 (MQ=255)
                                                                                                                               tgtgtgCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAAc                       >  2:2926565/1‑149 (MQ=255)
                                                                                                                               tgtgtgCCCACCGAACTGTCCGATATTTTAAGCAGTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAAc                       <  1:4614005/149‑1 (MQ=255)
                                                                                                                                   tgCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGga                   <  2:1826266/149‑1 (MQ=255)
                                                                                                                                    gCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGgac                  <  2:1356032/149‑1 (MQ=255)
                                                                                                                                      ccACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGgacga                >  2:843005/1‑149 (MQ=255)
                                                                                                                                        aCCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGgacgacg              >  1:4704377/1‑149 (MQ=255)
                                                                                                                                             aCTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAg         >  1:91357/1‑149 (MQ=255)
                                                                                                                                               tGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCg       <  1:2832139/149‑1 (MQ=255)
                                                                                                                                                gTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGa      <  1:4750827/149‑1 (MQ=255)
                                                                                                                                                 tCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGAc     <  1:558725/149‑1 (MQ=255)
                                                                                                                                                   cGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTg   >  1:3345112/1‑149 (MQ=255)
                                                                                                                                                    gATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTgc  >  1:3430489/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGC  >  CP009273/3994610‑3994906

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: