Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,523,001 T→A 60.5% intergenic (‑1/+55) yjhP ← / ← yjhQ putative methyltransferase/putative acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,523,0010TA60.5% 260.4 / inf 309intergenic (‑1/+55)yjhP/yjhQputative methyltransferase/putative acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (73/49);  new base A (90/97);  total (163/146)
Fisher's exact test for biased strand distribution p-value = 4.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

CTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAG  >  CP009273/4522854‑4523146
                                                                                                                                                   |                                                                                                                                                 
 tcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCattatt                                                                                                                                                >  1:1275828/1‑150 (MQ=255)
 tcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCattatt                                                                                                                                                <  1:1948285/150‑1 (MQ=255)
       gAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTc                                                                                                                                          >  2:909774/1‑150 (MQ=255)
       gAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTc                                                                                                                                          <  2:3909412/150‑1 (MQ=255)
        aGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCGTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTcc                                                                                                                                         <  2:2342591/150‑1 (MQ=255)
         gCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCt                                                                                                                                        <  2:1348054/150‑1 (MQ=255)
         gCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGAATATCCATAATTTAATTCCt                                                                                                                                        <  1:3326098/150‑1 (MQ=255)
           cGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGt                                                                                                                                      >  1:2059701/1‑150 (MQ=255)
             cTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGTTGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTa                                                                                                                                    <  1:2921614/150‑1 (MQ=255)
                ccGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTg                                                                                                                                 >  1:3614314/1‑150 (MQ=255)
                    ggTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGtt                                                                                                                             <  1:3213299/150‑1 (MQ=255)
                       cAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAg                                                                                                                          >  2:2302689/1‑150 (MQ=255)
                         aGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGcc                                                                                                                        >  2:1951110/1‑150 (MQ=255)
                          gAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCt                                                                                                                       <  1:1546156/150‑1 (MQ=255)
                          gAATGCGGGTGCCCGGATTCATGCTCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGTTGGTAAAAATACGTGGGATTTCCATAATTTAATTCCTGTTATTGGGTTAAGCCt                                                                                                                       <  1:3467655/150‑1 (MQ=255)
                           aaTGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCtt                                                                                                                      >  1:1589/1‑150 (MQ=255)
                           aaTGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCtt                                                                                                                      >  1:886965/1‑150 (MQ=255)
                              gCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTaaa                                                                                                                   <  2:2686555/150‑1 (MQ=255)
                                  gTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAgg                                                                                                               >  2:1455393/1‑150 (MQ=255)
                                   tGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGTACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTACTGTTATTGGGTTAAGCCTTAAAAAGGa                                                                                                              <  1:1000139/150‑1 (MQ=255)
                                     cccGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGaca                                                                                                            <  1:1169048/150‑1 (MQ=255)
                                     cccGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGaca                                                                                                            >  1:825645/1‑150 (MQ=255)
                                      ccGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGacac                                                                                                           <  2:4331920/150‑1 (MQ=255)
                                       cGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGAcaca                                                                                                          >  2:2597244/1‑150 (MQ=255)
                                        ggCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaa                                                                                                         >  2:2195665/1‑150 (MQ=255)
                                        ggCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaa                                                                                                         >  1:581294/1‑150 (MQ=255)
                                        ggCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaa                                                                                                         >  1:3256514/1‑150 (MQ=255)
                                         gCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaac                                                                                                        >  2:1420738/1‑150 (MQ=255)
                                          cttcttGCGCAGCACGCGAAGCAGCGTGGCGTTCTTCTCTTCGGTGTTCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATTTCCATTATTTAATTCCTGTTATTGGTTTAAGCCTTAAATAGGAGacaaca                                                                                                       <  2:3013307/150‑1 (MQ=255)
                                          cTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaaca                                                                                                       <  1:4237834/150‑1 (MQ=255)
                                          cTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaaca                                                                                                       <  1:4237614/150‑1 (MQ=255)
                                            tCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAg                                                                                                     <  2:34498/150‑1 (MQ=255)
                                               tGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGAtt                                                                                                  <  2:3622191/150‑1 (MQ=255)
                                                  gcagcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGtt                                                                                               >  2:4240306/1‑150 (MQ=255)
                                                    agcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTAt                                                                                             <  2:4101457/150‑1 (MQ=255)
                                                       aCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCa                                                                                          >  1:1882154/1‑150 (MQ=255)
                                                           gACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                      <  2:1855281/150‑1 (MQ=255)
                                                               cAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTg                                                                                  >  1:3395939/1‑150 (MQ=255)
                                                               cAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTg                                                                                  <  2:1137664/150‑1 (MQ=255)
                                                                 gCGTGGGGTACTTCTCTTCGGTGAACGGGTTTTGGATGCGGTGTTCACTTTCGATGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAATAAGGACACAACAAGATTGTTTTTCACGTCAGTGCt                                                                                <  2:3793996/150‑1 (MQ=255)
                                                                       cGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGtt                                                                          <  2:2131203/150‑1 (MQ=255)
                                                                       cGTACTTCTCGTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGtt                                                                          <  2:3281833/150‑1 (MQ=255)
                                                                        gTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGttt                                                                         >  1:1486248/1‑150 (MQ=255)
                                                                          aCTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtca                                                                       >  2:3961771/1‑150 (MQ=255)
                                                                          aCTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtca                                                                       >  2:3962764/1‑150 (MQ=255)
                                                                             tctcTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTTTTATTAAAGTCAGTGCTCCGGTTtcatca                                                                    >  2:2315629/1‑150 (MQ=255)
                                                                              ctctTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAg                                                                   >  2:1096495/1‑150 (MQ=255)
                                                                               tctTCGGTGCACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAgg                                                                  <  2:1176793/150‑1 (MQ=255)
                                                                                ctTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAggg                                                                 >  1:2194768/1‑150 (MQ=255)
                                                                                   cGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCt                                                              <  2:1486248/150‑1 (MQ=255)
                                                                                      tGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGAGCTCCGGTTTCATCAGGGTCTgcg                                                           >  2:3409520/1‑150 (MQ=255)
                                                                                      tGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTgcg                                                           <  1:3905513/150‑1 (MQ=255)
                                                                                      tGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTgcg                                                           <  2:3590464/150‑1 (MQ=255)
                                                                                      tGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTgcg                                                           <  2:1235773/150‑1 (MQ=255)
                                                                                      tGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTgcg                                                           >  1:250407/1‑150 (MQ=255)
                                                                                        aaCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCa                                                         <  1:3774494/150‑1 (MQ=255)
                                                                                          cGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCAcc                                                       <  2:168646/150‑1 (MQ=255)
                                                                                           gggTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCg                                                      >  1:1289032/1‑150 (MQ=255)
                                                                                                aTGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtgt                                                 >  1:3035013/1‑150 (MQ=255)
                                                                                                 tGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtgtg                                                >  1:472817/1‑150 (MQ=255)
                                                                                                  ggATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTgg                                               <  1:4196614/150‑1 (MQ=255)
                                                                                                     tGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCg                                            >  1:2256043/1‑150 (MQ=255)
                                                                                                      gCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACTTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCgg                                           >  2:2593218/1‑150 (MQ=255)
                                                                                                            ttCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAggg                                     >  2:4267280/1‑150 (MQ=255)
                                                                                                            ttCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAggg                                     >  2:3359640/1‑150 (MQ=255)
                                                                                                                  ttCGCTGTTGTTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAATTTTGTTATTCACGTCAGTGCTCCGGTTTAATCAGGGTCTGCGCACCGTATGTGGCCGGTCTGGGTCAttg                               <  1:753231/150‑3 (MQ=255)
                                                                                                                  ttCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATgg                               <  2:1829973/150‑1 (MQ=255)
                                                                                                                   tCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATggg                              <  1:2315584/150‑1 (MQ=255)
                                                                                                                   tCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATggg                              <  1:2009218/150‑1 (MQ=255)
                                                                                                                   tCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATggg                              <  1:955532/150‑1 (MQ=255)
                                                                                                                   tCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATggg                              >  2:3149891/1‑150 (MQ=255)
                                                                                                                    cGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGc                             >  1:325061/1‑150 (MQ=255)
                                                                                                                      cTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTg                           >  1:1496696/1‑150 (MQ=255)
                                                                                                                       tGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGa                          <  2:814795/150‑1 (MQ=255)
                                                                                                                       tGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGTGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGa                          <  1:2813561/150‑1 (MQ=255)
                                                                                                                       tGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGa                          <  1:527391/150‑1 (MQ=255)
                                                                                                                         aTGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGc                        <  1:3716896/150‑1 (MQ=255)
                                                                                                                          tGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCg                       >  1:668351/1‑150 (MQ=255)
                                                                                                                          tGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCg                       <  2:886965/150‑1 (MQ=255)
                                                                                                                           ggTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCgg                      <  2:2259756/150‑1 (MQ=255)
                                                                                                                            gTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgt                     <  2:1876052/150‑1 (MQ=255)
                                                                                                                            gTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgt                     <  1:486503/150‑1 (MQ=255)
                                                                                                                              aaaaaTACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtga                   >  2:67966/1‑150 (MQ=255)
                                                                                                                              aaaaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtga                   >  1:2763286/1‑150 (MQ=255)
                                                                                                                                   tACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGAt              <  1:3438263/150‑1 (MQ=255)
                                                                                                                                   tACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGAt              >  2:3186516/1‑150 (MQ=255)
                                                                                                                                     cGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTg            >  1:4522029/1‑150 (MQ=255)
                                                                                                                                        gggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAATGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcat         >  2:3614376/1‑150 (MQ=255)
                                                                                                                                         ggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc        >  2:2532167/1‑150 (MQ=255)
                                                                                                                                         ggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc        >  2:2460724/1‑150 (MQ=255)
                                                                                                                                         ggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc        >  2:105629/1‑150 (MQ=255)
                                                                                                                                         ggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc        >  2:159181/1‑150 (MQ=255)
                                                                                                                                         ggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc        >  1:2585834/1‑150 (MQ=255)
                                                                                                                                         ggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc        <  1:4348414/150‑1 (MQ=255)
                                                                                                                                         ggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTACTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc        >  1:3870026/1‑150 (MQ=255)
                                                                                                                                           atatCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatcat      >  2:3746135/1‑150 (MQ=255)
                                                                                                                                              tCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCca   <  1:3098778/150‑1 (MQ=255)
                                                                                                                                               ccATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcag  >  1:3843710/1‑150 (MQ=255)
                                                                                                                                               ccATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcag  >  1:245404/1‑150 (MQ=255)
                                                                                                                                               ccATAATTTAATTCCTGTTATTGGGTTAAGCCTAAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcag  >  2:639341/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAG  >  CP009273/4522854‑4523146

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: