Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,984,342 G→T 57.8% intergenic (‑216/‑171) hemC ← / → cyaA hydroxymethylbilane synthase/adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,984,3420GT57.8% 37.4 / inf 344intergenic (‑216/‑171)hemC/cyaAhydroxymethylbilane synthase/adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (55/90);  new base T (93/106);  total (148/196)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 344 total aligned reads displayed.

AAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAC  >  CP009273/3984197‑3984488
                                                                                                                                                 |                                                                                                                                                  
 aaCGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGTATCATGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTTCTACACTTGTATGTAGCGCATATTTATTTACGGTAAATCAGCAATGTGtt                                                                                                                                               >  2:2296447/1‑150 (MQ=255)
 aaCGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGtt                                                                                                                                               >  1:21058/1‑150 (MQ=255)
   cGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTaa                                                                                                                                             <  1:4170291/150‑1 (MQ=255)
     tAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAAt                                                                                                                                           <  2:2665850/150‑1 (MQ=255)
       gACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAAAGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTGAATCAGCAATGTGTTAAATTg                                                                                                                                         <  2:2558763/150‑1 (MQ=255)
         cGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGAt                                                                                                                                       >  1:1691435/1‑150 (MQ=255)
          gCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATc                                                                                                                                      <  1:3334608/150‑1 (MQ=255)
            cTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCAc                                                                                                                                    <  2:3618112/150‑1 (MQ=255)
               cgcAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACATAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGtt                                                                                                                                 <  2:2328471/150‑1 (MQ=255)
                   aaCCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAg                                                                                                                             <  1:2838646/150‑1 (MQ=255)
                    aCCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGa                                                                                                                            >  2:3010008/1‑150 (MQ=255)
                     ccGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTTCACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGAc                                                                                                                           <  1:50382/150‑1 (MQ=255)
                      cGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGAcc                                                                                                                          >  1:1843938/1‑150 (MQ=255)
                         aaaaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCAtt                                                                                                                       <  1:948118/150‑1 (MQ=255)
                         aaaaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCAtt                                                                                                                       <  1:2968130/150‑1 (MQ=255)
                          aaaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAttt                                                                                                                      <  2:1869569/150‑1 (MQ=255)
                            aTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAttttt                                                                                                                    >  1:1576079/1‑150 (MQ=255)
                               aGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTtcg                                                                                                                 <  2:600332/150‑1 (MQ=255)
                               aGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTtcg                                                                                                                 <  2:3929622/150‑1 (MQ=255)
                                  tCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTtcgtcg                                                                                                              <  1:550551/150‑1 (MQ=255)
                                  tCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTtcgtcg                                                                                                              <  1:2459620/150‑1 (MQ=255)
                                    tGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTg                                                                                                            <  1:1020843/150‑1 (MQ=255)
                                        aGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAAc                                                                                                        >  2:1269536/1‑150 (MQ=255)
                                        aGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAAc                                                                                                        <  1:3569324/150‑1 (MQ=255)
                                           ggATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaa                                                                                                     >  1:4337116/1‑150 (MQ=255)
                                           ggATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaa                                                                                                     <  2:713239/150‑1 (MQ=255)
                                            gATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaaa                                                                                                    <  2:1650157/150‑1 (MQ=255)
                                             aTCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaaaa                                                                                                   >  1:980933/1‑150 (MQ=255)
                                                 tGACAGGCGTTTCTCTCCGTTGTAATAAGGAATTTACAGTGAATAATCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTTAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaaaaattc                                                                                               <  2:3580612/150‑4 (MQ=255)
                                                 tGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAAc                                                                                               <  1:1881088/150‑1 (MQ=255)
                                                 tGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAAc                                                                                               >  1:4115963/1‑150 (MQ=255)
                                                 tGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAAc                                                                                               <  2:3743985/150‑1 (MQ=255)
                                                  gACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAAcc                                                                                              <  1:1078832/150‑1 (MQ=255)
                                                      ggCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGgc                                                                                          <  2:1788661/150‑1 (MQ=255)
                                                       gCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGgcg                                                                                         <  2:4216105/150‑1 (MQ=255)
                                                         gTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGcgcg                                                                                       >  2:1770378/1‑150 (MQ=255)
                                                            tCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAAACAGGCGCGaaa                                                                                    >  2:4420421/1‑150 (MQ=255)
                                                            tCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGaaa                                                                                    >  2:3405856/1‑150 (MQ=255)
                                                               cGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGt                                                                                 >  1:2269178/1‑150 (MQ=255)
                                                                    ttGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTaa                                                                            >  2:74221/1‑150 (MQ=255)
                                                                     tGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAAACAGGCGCGAAAAGTGGTAAc                                                                           >  1:467329/1‑150 (MQ=255)
                                                                     tGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGATTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAAc                                                                           <  1:2156735/150‑1 (MQ=255)
                                                                        aataaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGt                                                                        >  2:2010398/1‑150 (MQ=255)
                                                                         ataaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTCGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGtt                                                                       >  1:648187/1‑150 (MQ=255)
                                                                         ataaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGtt                                                                       <  1:127495/150‑1 (MQ=255)
                                                                          taaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTa                                                                      >  1:3423908/1‑150 (MQ=255)
                                                                          taaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTa                                                                      <  1:906144/150‑1 (MQ=255)
                                                                            aGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTAcc                                                                    <  2:1988168/150‑1 (MQ=255)
                                                                                aTTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTg                                                                >  2:3428083/1‑150 (MQ=255)
                                                                                aTTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTg                                                                <  2:2120044/150‑1 (MQ=255)
                                                                                  ttACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGAc                                                              >  1:1496632/1‑150 (MQ=255)
                                                                                  ttACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGAc                                                              >  2:4512975/1‑150 (MQ=255)
                                                                                   tACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACa                                                             >  1:799011/1‑150 (MQ=255)
                                                                                     cAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATa                                                           >  2:1185497/1‑150 (MQ=255)
                                                                                      agagAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATAc                                                          <  2:3859545/150‑1 (MQ=255)
                                                                                      agagAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATAc                                                          >  2:4017345/1‑150 (MQ=255)
                                                                                      agagAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATAc                                                          <  1:2578433/150‑1 (MQ=255)
                                                                                      agagAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATAc                                                          >  1:2742306/1‑150 (MQ=255)
                                                                                      agagAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTAAGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATAc                                                          <  2:1350693/150‑1 (MQ=255)
                                                                                       gagaATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACg                                                         <  1:3019860/150‑1 (MQ=255)
                                                                                       gagaATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACg                                                         <  1:3130381/150‑1 (MQ=255)
                                                                                           aTAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAat                                                     >  1:1507945/1‑150 (MQ=255)
                                                                                            tAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAata                                                    <  1:2942497/150‑1 (MQ=255)
                                                                                            tAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAata                                                    <  1:1414541/150‑1 (MQ=255)
                                                                                             aaaCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAatat                                                   <  2:869319/150‑1 (MQ=255)
                                                                                             aaaCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAatat                                                   <  2:3209205/150‑1 (MQ=255)
                                                                                               aCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATcc                                                 <  2:3320983/150‑1 (MQ=255)
                                                                                                cGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATccc                                                >  2:1582062/1‑150 (MQ=255)
                                                                                                 ggTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCg                                               >  2:564148/1‑150 (MQ=255)
                                                                                                 ggTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCg                                               <  1:4024136/150‑1 (MQ=255)
                                                                                                    gCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAAt                                            >  1:4621543/1‑150 (MQ=255)
                                                                                                     cTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATg                                           <  1:1909793/150‑1 (MQ=255)
                                                                                                        cacTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcg                                        <  2:1931270/150‑1 (MQ=255)
                                                                                                        cacTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcg                                        >  1:3783422/1‑150 (MQ=255)
                                                                                                           ttGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCgt                                     >  1:1061528/1‑150 (MQ=255)
                                                                                                           ttGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCgt                                     >  1:166430/1‑150 (MQ=255)
                                                                                                           ttGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCgt                                     >  1:2790888/1‑150 (MQ=255)
                                                                                                           ttGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCgt                                     <  2:21058/150‑1 (MQ=255)
                                                                                                            tgtatgtaGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCgtg                                    <  1:4176314/150‑1 (MQ=255)
                                                                                                                 gtaGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTAcc                               <  1:3901765/150‑1 (MQ=255)
                                                                                                                 gtaGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTAcc                               <  2:2505054/150‑1 (MQ=255)
                                                                                                                 gtaGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTAcc                               <  2:1148381/150‑1 (MQ=255)
                                                                                                                    gcgcATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGtt                            >  1:2393783/1‑150 (MQ=255)
                                                                                                                    gcgcATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGtt                            <  2:1104324/150‑1 (MQ=255)
                                                                                                                      gcATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGa                          <  2:189489/150‑1 (MQ=255)
                                                                                                                       cATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGAt                         <  1:565683/150‑1 (MQ=255)
                                                                                                                       cATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGAt                         >  1:3653508/1‑150 (MQ=255)
                                                                                                                        aTCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATg                        >  2:1438988/1‑150 (MQ=255)
                                                                                                                         tctttcttTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGt                       >  1:4405968/1‑150 (MQ=255)
                                                                                                                         tctttcttTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGt                       >  2:2710011/1‑150 (MQ=255)
                                                                                                                          ctttctttACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGtt                      >  1:1446151/1‑150 (MQ=255)
                                                                                                                           tttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTg                     >  1:1803937/1‑150 (MQ=255)
                                                                                                                                ttACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGa                >  2:283823/1‑150 (MQ=255)
                                                                                                                                   cGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATc             >  2:4360883/1‑150 (MQ=255)
                                                                                                                                   cGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATc             >  2:4248919/1‑150 (MQ=255)
                                                                                                                                   cGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATc             >  2:4085387/1‑150 (MQ=255)
                                                                                                                                   cGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATc             >  2:831973/1‑150 (MQ=255)
                                                                                                                                    ggTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATca            >  2:4123501/1‑150 (MQ=255)
                                                                                                                                     gTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATcac           >  1:3243329/1‑150 (MQ=255)
                                                                                                                                              cAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAc  <  1:121380/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAC  >  CP009273/3984197‑3984488

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: