Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,601,829 T→G 6.1% S518S (TCT→TCG hyfF → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,601,8290TG 6.1% 823.7 / 40.7 265S518S (TCT→TCGhyfFhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (136/113);  new base G (10/6);  total (146/119)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

Only 100 of 270 total aligned reads displayed.

CTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAAATCAATCTGC  >  CP009273/2601680‑2601978
                                                                                                                                                     |                                                                                                                                                     
ctgGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCt                                                                                                                                                       <  1:2919379/150‑1 (MQ=255)
   gTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGtt                                                                                                                                                    >  2:368582/1‑150 (MQ=255)
     cATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTc                                                                                                                                                  <  2:3344938/150‑1 (MQ=255)
      atgatgCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCg                                                                                                                                                 <  1:324346/150‑1 (MQ=255)
       tgatgCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTCTAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCgg                                                                                                                                                >  2:950256/1‑150 (MQ=255)
         atgCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGgag                                                                                                                                              >  1:3060046/1‑150 (MQ=255)
          tgCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGgaga                                                                                                                                             <  2:2127221/150‑1 (MQ=255)
            cTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGaaa                                                                                                                                           >  1:3351056/1‑150 (MQ=255)
              ggCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCGGTTTCGGAGAAaca                                                                                                                                         >  1:3388863/1‑150 (MQ=255)
                 gATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAg                                                                                                                                      >  1:1118772/1‑150 (MQ=255)
                 gATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCGGTTTCGGAGAAACACAg                                                                                                                                      >  1:2216418/1‑150 (MQ=255)
                  aTGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGt                                                                                                                                     <  1:3379430/150‑1 (MQ=255)
                  aTGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGt                                                                                                                                     >  2:2525897/1‑150 (MQ=255)
                  aTGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCACGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGt                                                                                                                                     >  1:2122883/1‑150 (MQ=255)
                  aTGGGAACGCATATTCCACAAACTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCTGGCACCGCATCTGTTTCGGAGAAACACAGt                                                                                                                                     >  2:2462486/1‑150 (MQ=255)
                   tGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTg                                                                                                                                    >  1:2204688/1‑150 (MQ=255)
                    gggAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGa                                                                                                                                   >  1:665746/1‑150 (MQ=255)
                       aaCGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACg                                                                                                                                <  2:2909570/150‑1 (MQ=255)
                         cGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGtt                                                                                                                              <  2:4064792/150‑1 (MQ=255)
                         cGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGtt                                                                                                                              <  2:1387653/150‑1 (MQ=255)
                         cGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGtt                                                                                                                              >  2:125165/1‑150 (MQ=255)
                            atatTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAAt                                                                                                                           >  1:1854422/1‑150 (MQ=255)
                             tatTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAACACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAAtt                                                                                                                          >  2:1281558/1‑150 (MQ=255)
                              atTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTc                                                                                                                         <  2:4179050/150‑1 (MQ=255)
                                 ccACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATc                                                                                                                      <  2:2372966/150‑1 (MQ=255)
                                 ccACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGAACGATCTGCCTGCACAACGTAGCACATGGCATGATTTTTTTCCTTCAGGCACCGCATCGGTTTCGGAGAAACACAGTGAACGTTAATTCATc                                                                                                                      <  1:3597419/150‑1 (MQ=255)
                                   acaACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGt                                                                                                                    <  2:4013811/150‑1 (MQ=255)
                                    caACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTc                                                                                                                   >  2:3211937/1‑150 (MQ=255)
                                    caACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTc                                                                                                                   >  2:549017/1‑150 (MQ=255)
                                     aaCCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCa                                                                                                                  <  1:620596/150‑1 (MQ=255)
                                     aaCCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCa                                                                                                                  >  2:461882/1‑150 (MQ=255)
                                     aaCCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCa                                                                                                                  >  2:359808/1‑150 (MQ=255)
                                       ccTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCaaa                                                                                                                >  1:4419128/1‑150 (MQ=255)
                                       ccTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCaaa                                                                                                                <  1:3142043/150‑1 (MQ=255)
                                        cTGTCATCAGGATCCTGGCGGGCGCTTCCACAATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTAAGGCACCGCAACTGTTTCGGAGAAACACAGTGAAAGTTAAATAATAGTCAAAt                                                                                                               >  1:3671324/1‑150 (MQ=255)
                                          gTCATCAGGATCCTGGCGGGCGCTTCCACGATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGGTAATTCAAAGTCAAATCg                                                                                                             >  1:3791501/1‑150 (MQ=255)
                                             atcaGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTgg                                                                                                          <  2:3962894/150‑1 (MQ=255)
                                             atcaGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTgg                                                                                                          >  1:2645192/1‑150 (MQ=255)
                                               caGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCg                                                                                                        <  1:3224076/150‑1 (MQ=255)
                                                aGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGa                                                                                                       >  1:4046903/1‑150 (MQ=255)
                                                 ggATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGaa                                                                                                      <  2:44972/150‑1 (MQ=255)
                                                          cgggcgCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATtctc                                                                                             <  1:4138278/150‑1 (MQ=255)
                                                          cgggcgCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATtctc                                                                                             >  1:1700115/1‑150 (MQ=255)
                                                             gcgcTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTcgcc                                                                                          <  1:1628878/150‑1 (MQ=255)
                                                              cgcTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTcgccg                                                                                         >  2:1060015/1‑150 (MQ=255)
                                                              cgcTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTcgccg                                                                                         <  2:302969/150‑1 (MQ=255)
                                                              cgcTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTcgccg                                                                                         >  1:1425515/1‑150 (MQ=255)
                                                                 ttCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGccc                                                                                      >  2:4424754/1‑150 (MQ=255)
                                                                    cACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGa                                                                                   >  1:441390/1‑150 (MQ=255)
                                                                        ataGTCCGCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAAc                                                                               <  2:3382541/150‑1 (MQ=255)
                                                                         taGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACg                                                                              >  2:2818010/1‑150 (MQ=255)
                                                                          aGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGc                                                                             >  2:2414285/1‑150 (MQ=255)
                                                                             ccTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGt                                                                          >  1:2068982/1‑150 (MQ=255)
                                                                               tctcAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTc                                                                        <  1:2419392/150‑1 (MQ=255)
                                                                                ctcAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTcc                                                                       <  2:2801435/150‑1 (MQ=255)
                                                                                ctcAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTcc                                                                       <  1:2643289/150‑1 (MQ=255)
                                                                                   aGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCg                                                                    <  2:3564154/150‑1 (MQ=255)
                                                                                   aGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATAGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCg                                                                    >  1:4460859/1‑150 (MQ=255)
                                                                                    gggACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCgg                                                                   <  1:1045143/150‑1 (MQ=255)
                                                                                      gACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGgcg                                                                 <  1:582440/150‑1 (MQ=255)
                                                                                        cGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGcgcg                                                               <  2:875920/150‑1 (MQ=255)
                                                                                          cACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGt                                                             >  1:4427612/1‑150 (MQ=255)
                                                                                           aCGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTg                                                            <  1:847252/150‑1 (MQ=255)
                                                                                                 tgcctgcACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCGGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGAt                                                      <  1:1882284/150‑1 (MQ=255)
                                                                                                  gcctgcACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATg                                                     <  1:717624/150‑1 (MQ=255)
                                                                                                   cctgcACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGa                                                    >  2:666559/1‑150 (MQ=255)
                                                                                                    ctgcACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGaa                                                   >  1:2474316/1‑150 (MQ=255)
                                                                                                    ctgcACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGaa                                                   <  1:1230728/150‑1 (MQ=255)
                                                                                                      gcACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAaga                                                 <  1:3773476/150‑1 (MQ=255)
                                                                                                       cacaACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAagag                                                <  1:461882/150‑1 (MQ=255)
                                                                                                        acaACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGc                                               >  2:1938079/1‑150 (MQ=255)
                                                                                                         caACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCg                                              <  1:4277416/150‑1 (MQ=255)
                                                                                                             gTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCGGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACaa                                          >  2:4463317/1‑150 (MQ=255)
                                                                                                             gTAGCACCTGGCATGATTTTATGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCATGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACaa                                          >  1:1469597/1‑150 (MQ=255)
                                                                                                              tAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACaaa                                         >  2:2256317/1‑150 (MQ=255)
                                                                                                                gCACCTGGCATGATTTTTTGCCTTTAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACg                                       <  1:1342400/150‑1 (MQ=255)
                                                                                                                  aCCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGcc                                     <  1:1491789/150‑1 (MQ=255)
                                                                                                                  aCCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGcc                                     >  1:920556/1‑150 (MQ=255)
                                                                                                                  aCCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGcc                                     >  1:1830640/1‑150 (MQ=255)
                                                                                                                     tGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGa                                  <  1:4009788/150‑1 (MQ=255)
                                                                                                                         aTGATTTTTTGCCTTCAGGCACCGCATCGGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAg                              >  2:317869/1‑150 (MQ=255)
                                                                                                                            aTTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTc                           <  1:1275484/150‑1 (MQ=255)
                                                                                                                            aTTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTc                           <  2:1919750/150‑1 (MQ=255)
                                                                                                                            aTTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTc                           <  2:1576389/150‑1 (MQ=255)
                                                                                                                            aTTTTTGGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTc                           <  2:3883474/150‑1 (MQ=255)
                                                                                                                             ttttttGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCa                          <  1:278515/150‑1 (MQ=255)
                                                                                                                               ttttGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCAcc                        <  1:222262/150‑1 (MQ=255)
                                                                                                                                tttGCCTTCAGGCACCGCATCTGTGTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCa                       <  2:2674938/150‑1 (MQ=255)
                                                                                                                                 ttGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCAt                      <  2:3581205/150‑1 (MQ=255)
                                                                                                                                   gCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTCATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTa                    >  2:2792738/1‑150 (MQ=255)
                                                                                                                                     cTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACg                  <  2:3774786/150‑1 (MQ=255)
                                                                                                                                       tCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGt                <  2:4283747/150‑1 (MQ=255)
                                                                                                                                               cGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAaatca        >  2:1775142/1‑150 (MQ=255)
                                                                                                                                               cGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAaatca        >  1:4106492/1‑150 (MQ=255)
                                                                                                                                                 cATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAaatcaat      <  2:4265059/150‑1 (MQ=255)
                                                                                                                                                 cATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAaatcaat      >  2:779554/1‑150 (MQ=255)
                                                                                                                                                  aTCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAaatcaatc     <  1:2818022/150‑1 (MQ=255)
                                                                                                                                                   tCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAAATCAATct    >  1:3847216/1‑150 (MQ=255)
                                                                                                                                                    cTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAAATCAATctg   >  1:712166/1‑150 (MQ=255)
                                                                                                                                                     tGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTAGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAAATCAATctgc  >  1:3265457/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCATCTGTTTCGGAGAAACACAGTGAACGTTAATTCATCGTCAAATCGTGGCGAAGCGATTCTCGCCGCCCTGAAAACGCAGTTCCCCGGCGCGGTGCTGGATGAAGAGCGACAAACGCCTGAACAGGTCACCATTACGGTGAAAATCAATCTGC  >  CP009273/2601680‑2601978

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: