Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,045,386 G→T 22.1% intergenic (+170/‑19) yihI → / → hemN activator of Der GTPase/coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,045,3860GT22.1% 625.6 / 117.6 245intergenic (+170/‑19)yihI/hemNactivator of Der GTPase/coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base G (77/113);  new base T (20/35);  total (97/148)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

TCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGTTTTTACAA  >  CP009273/4045243‑4045533
                                                                                                                                               |                                                                                                                                                   
tCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCAt                                                                                                                                                 >  2:1666993/1‑148 (MQ=255)
    ctgATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCTGa                                                                                                                                             <  2:2902675/148‑1 (MQ=255)
       aTAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAgg                                                                                                                                          >  1:2560677/1‑148 (MQ=255)
        gaaCCCTCCCGGGTATCTGTTATCTGGGGCGGGTAGTCGAGAAACTCCCGCTGTGGTCGCGGCGACAAAAGTGGCTGCGCACCCGGGTGATGACCCGAAACGGGCCTCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGa                                                                                                                                         <  2:3365734/147‑1 (MQ=255)
          aCCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGt                                                                                                                                       >  2:456733/1‑148 (MQ=255)
          aCCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGt                                                                                                                                       <  2:1281852/148‑1 (MQ=255)
          aCCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGt                                                                                                                                       <  1:373139/148‑1 (MQ=255)
           cccTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTg                                                                                                                                      <  2:1715381/148‑1 (MQ=255)
           cccTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTg                                                                                                                                      >  1:3672920/1‑148 (MQ=255)
              ttCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGc                                                                                                                                   >  1:1391014/1‑148 (MQ=255)
               tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCa                                                                                                                                  >  1:192357/1‑148 (MQ=255)
               tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCa                                                                                                                                  <  2:2280410/148‑1 (MQ=255)
               tCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCa                                                                                                                                  >  2:3010440/1‑148 (MQ=255)
                 cGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATg                                                                                                                                <  1:2239472/148‑1 (MQ=255)
                     tataTGTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTg                                                                                                                            >  2:949817/1‑148 (MQ=255)
                          gTTATCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAg                                                                                                                       >  2:557830/1‑148 (MQ=255)
                             aTCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                                                                                                                    <  2:2275128/148‑1 (MQ=255)
                              tCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaaa                                                                                                                   <  2:1430738/148‑1 (MQ=255)
                              tCTGGTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaaa                                                                                                                   <  1:2531590/148‑1 (MQ=255)
                                 ggTGTGGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCg                                                                                                                >  2:3062735/1‑148 (MQ=255)
                                   tgtgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGAc                                                                                                              >  2:1827090/1‑148 (MQ=255)
                                    gtgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACt                                                                                                             >  1:2326460/1‑148 (MQ=255)
                                    gtgGTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACt                                                                                                             >  1:2241610/1‑148 (MQ=255)
                                      ggTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTgg                                                                                                           >  1:963861/1‑148 (MQ=255)
                                      ggTTATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTgg                                                                                                           >  1:2875675/1‑148 (MQ=255)
                                        ttATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGa                                                                                                         >  1:1307587/1‑148 (MQ=255)
                                         tATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCGTCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGAt                                                                                                        >  2:177910/1‑148 (MQ=255)
                                         tATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGAt                                                                                                        >  1:2189148/1‑148 (MQ=255)
                                         tATTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGAt                                                                                                        >  1:3784402/1‑148 (MQ=255)
                                          aTTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATc                                                                                                       <  2:195554/148‑1 (MQ=255)
                                          aTTCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATc                                                                                                       <  2:1083490/148‑1 (MQ=255)
                                           ttCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCt                                                                                                      <  2:2103756/148‑1 (MQ=255)
                                           ttCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCt                                                                                                      <  2:2208116/148‑1 (MQ=255)
                                           ttCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                                                                                                                    <  1:3217625/134‑1 (MQ=255)
                                           ttCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCaa                                                                                                                    >  2:3217625/1‑134 (MQ=255)
                                           ttCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCt                                                                                                      >  2:341700/1‑148 (MQ=255)
                                            tCGTTAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTg                                                                                                     >  2:3767414/1‑148 (MQ=255)
                                               ttAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc                                                                                                  <  1:1274962/148‑1 (MQ=255)
                                               ttAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGcc                                                                                                  <  2:889870/148‑1 (MQ=255)
                                                tAAACTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGccc                                                                                                 >  2:1141135/1‑148 (MQ=255)
                                                  aaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTg                                                                                               <  2:3591513/148‑1 (MQ=255)
                                                  aaCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTg                                                                                               >  2:725862/1‑148 (MQ=255)
                                                   aCTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGa                                                                                              <  1:2412416/148‑1 (MQ=255)
                                                    cTACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGAt                                                                                             >  1:312410/1‑148 (MQ=255)
                                                     tACAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATc                                                                                            >  2:28116/1‑148 (MQ=255)
                                                      aCAGCTGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATcc                                                                                           >  1:678817/1‑148 (MQ=255)
                                                            gTTGTCGCGGCGACAAAAGTGGCTGCGCACCCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAat                                                                                     <  1:3523259/148‑1 (MQ=255)
                                                            gTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGTAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAat                                                                                     <  1:1633197/148‑1 (MQ=255)
                                                             ttGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAata                                                                                    <  2:3250799/148‑1 (MQ=255)
                                                              tGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAatat                                                                                   <  2:1493963/148‑1 (MQ=255)
                                                               gTCGGGGCGACAAATGTGGCTGCGAAACGGGGTGATGACCCGTAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGTGACGCCATCGGTAGGAGTTAGCTTGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAAtata                                                                                  <  1:2967041/148‑1 (MQ=255)
                                                               gTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAAtata                                                                                  <  1:275232/148‑1 (MQ=255)
                                                                tCGCGGCGACAAAAGTGGCTGCGCACCGGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGTGCATTTCTGTACAGCAAATCGACTGGTATCTGGCCCTGATCCAGAAATATaa                                                                                 <  1:3290795/148‑1 (MQ=255)
                                                                tCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATaa                                                                                 <  1:531985/148‑1 (MQ=255)
                                                                tCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATaa                                                                                 <  2:505510/148‑1 (MQ=255)
                                                                 cgcgGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAAc                                                                                <  2:39078/148‑1 (MQ=255)
                                                                  gcggcgACAAATGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACt                                                                               <  1:3706651/148‑1 (MQ=255)
                                                                   cggcgACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTa                                                                              <  1:3595697/148‑1 (MQ=255)
                                                                   cggcgACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTa                                                                              <  1:535338/148‑1 (MQ=255)
                                                                     gcgACAAAAGTGGCTGCGCAACCGGCTTATGACCCTAAACGGGCATCCGCCGGTACGCCGTAGCCGCCATAGACTCCATCGGAAGTAGTGAGCATGTCTGTACAGAAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTAtt                                                                            <  1:599249/148‑1 (MQ=255)
                                                                      cgACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTc                                                                           <  2:1553393/148‑1 (MQ=255)
                                                                      cgACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTc                                                                           <  2:3450958/148‑1 (MQ=255)
                                                                        aCAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCg                                                                         >  2:1518451/1‑148 (MQ=255)
                                                                               tGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTTCGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACg                                                                  <  2:3670737/148‑1 (MQ=255)
                                                                                      gcAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACAcc                                                           <  1:1318756/148‑1 (MQ=255)
                                                                                        aaCCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTc                                                         <  2:2875675/148‑1 (MQ=255)
                                                                                         aCCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCg                                                        >  1:1082153/1‑148 (MQ=255)
                                                                                         aCCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCg                                                        <  1:2804874/148‑1 (MQ=255)
                                                                                         aCCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCg                                                        <  1:891343/148‑1 (MQ=255)
                                                                                             gCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTAcc                                                    >  2:1653709/1‑148 (MQ=255)
                                                                                             gCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTAcc                                                    >  2:2573409/1‑148 (MQ=255)
                                                                                                gatgaCCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGa                                                 >  1:3744895/1‑148 (MQ=255)
                                                                                                gatgaCCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGa                                                 >  2:2768773/1‑148 (MQ=255)
                                                                                                     cccGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACcgcg                                            <  1:2841945/148‑1 (MQ=255)
                                                                                                      ccGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCgcgc                                           >  1:736776/1‑148 (MQ=255)
                                                                                                       cGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCt                                          >  1:450523/1‑148 (MQ=255)
                                                                                                          aaCGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGa                                       >  1:3038465/1‑148 (MQ=255)
                                                                                                           aCGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAg                                      >  2:1103671/1‑148 (MQ=255)
                                                                                                           aCGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAg                                      >  2:3707344/1‑148 (MQ=255)
                                                                                                           aCGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAg                                      <  2:3672047/148‑1 (MQ=255)
                                                                                                             gggCATCCGCCGGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGtt                                    <  1:2338385/148‑1 (MQ=255)
                                                                                                                     gccgGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCagaag                            <  2:3349803/148‑1 (MQ=255)
                                                                                                                     gccgGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCagaag                            <  2:3585879/148‑1 (MQ=255)
                                                                                                                      ccgGTACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCagaaga                           <  1:2082176/148‑1 (MQ=255)
                                                                                                                         gTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACtt                        <  1:1763071/148‑1 (MQ=255)
                                                                                                                         gTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACtt                        <  1:2836662/148‑1 (MQ=255)
                                                                                                                          tACGCCGTAGCCGCCAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTc                       >  1:698313/1‑148 (MQ=255)
                                                                                                                          tACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTc                       <  1:1742588/148‑1 (MQ=255)
                                                                                                                                 tAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAAc                <  1:3010440/148‑1 (MQ=255)
                                                                                                                                   gccgccAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAg              >  2:781720/1‑148 (MQ=255)
                                                                                                                                      gccAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCg           <  1:2185606/148‑1 (MQ=255)
                                                                                                                                      gccAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCg           <  1:2838881/148‑1 (MQ=255)
                                                                                                                                      gccAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCg           <  1:2116337/148‑1 (MQ=255)
                                                                                                                                       ccAGAGACTCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGt          <  1:3212602/148‑1 (MQ=255)
                                                                                                                                       ccAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGt          >  2:1820512/1‑148 (MQ=255)
                                                                                                                                          gagaCGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGtttt       >  2:1416698/1‑148 (MQ=255)
                                                                                                                                           agaCGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGttttt      <  1:1046967/148‑1 (MQ=255)
                                                                                                                                              cGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGTTTTTACa   <  2:533442/148‑1 (MQ=255)
                                                                                                                                              cGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGTTTTTACa   <  1:880383/148‑1 (MQ=255)
                                                                                                                                              cGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGTTTTTACa   <  2:2456331/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCTGCTGATAACCCTTCCGGTTATATGTTATCTGGTGTGGTTATTCGTTAAACTACAGCGGTTGTCGCGGCGACAAAAGTGGCTGCGCAACCGGCTGATGACCCGAAACGGGCATCCGCCGGTACGCCGTAGCCGCCAGAGACGCCATCGGAAGGAGTGAGCATGTCTGTACAGCAAATCGACTGGGATCTGGCCCTGATCCAGAAATATAACTATTCCGGGCCACGATACACCTCGTACCCGACCGCGCTGGAGTTTTCAGAAGACTTCGGCGAACAGGCGTTTTTACAA  >  CP009273/4045243‑4045533

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