Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,523,001 T→A 88.7% 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,0010TA88.7% 415.8 / 37.8 142intergenic (‑1/+55)yjhP/yjhQputative methyltransferase/putative acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/14);  new base A (43/83);  total (45/97)
Fisher's exact test for biased strand distribution p-value = 9.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 150 total aligned reads displayed.

CGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGC  >  CP009273/4522856‑4523147
                                                                                                                                                 |                                                                                                                                                  
cGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCataat                                                                                                                                                  <  1:3189612/148‑1 (MQ=255)
cGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCataat                                                                                                                                                  <  2:1820461/148‑1 (MQ=255)
cGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCataat                                                                                                                                                  >  1:1750307/1‑148 (MQ=255)
   ccGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTa                                                                                                                                               >  2:1324975/1‑148 (MQ=255)
    cGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTaa                                                                                                                                              <  2:2247287/148‑1 (MQ=255)
        ccGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTcc                                                                                                                                          >  1:3054789/1‑148 (MQ=255)
         cGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCt                                                                                                                                         <  1:144748/148‑1 (MQ=255)
          gCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTg                                                                                                                                        <  1:1320575/148‑1 (MQ=255)
           cTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGt                                                                                                                                       <  2:1968797/148‑1 (MQ=255)
           cTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGt                                                                                                                                       <  2:1968801/148‑1 (MQ=255)
            tGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGtt                                                                                                                                      <  2:2490059/148‑1 (MQ=255)
            tGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGtt                                                                                                                                      <  2:3181930/148‑1 (MQ=255)
            tGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGtt                                                                                                                                      <  2:3182400/148‑1 (MQ=255)
              ccGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTAt                                                                                                                                    <  1:972875/148‑1 (MQ=255)
                gAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTg                                                                                                                                  <  1:2694794/148‑1 (MQ=255)
                gAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTg                                                                                                                                  <  1:2694886/148‑1 (MQ=255)
                 aGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCAATAAGTTAATTCCTGTTATTgg                                                                                                                                 >  2:2012143/1‑148 (MQ=255)
                     cAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTa                                                                                                                             >  2:194598/1‑148 (MQ=255)
                     cAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTAGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTa                                                                                                                             <  2:404084/148‑1 (MQ=255)
                      aaGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTaa                                                                                                                            <  2:3122190/148‑1 (MQ=255)
                         aaTGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGcc                                                                                                                         >  1:282990/1‑148 (MQ=255)
                         aaTGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGcc                                                                                                                         <  2:2499520/148‑1 (MQ=255)
                          aTGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCt                                                                                                                        <  1:2011003/148‑1 (MQ=255)
                            gCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTa                                                                                                                      <  1:1314124/148‑1 (MQ=255)
                            gCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTa                                                                                                                      <  1:834895/148‑1 (MQ=255)
                            gCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTa                                                                                                                      <  1:834826/148‑1 (MQ=255)
                               ggTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTaaaa                                                                                                                   >  1:2474268/1‑148 (MQ=255)
                               ggTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTaaaa                                                                                                                   <  1:2797196/148‑1 (MQ=255)
                                 tGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTAATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAg                                                                                                                 <  2:600439/148‑1 (MQ=255)
                                   cccGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGa                                                                                                               <  1:2650050/148‑1 (MQ=255)
                                    ccGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTCTGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGac                                                                                                              <  2:1858093/148‑1 (MQ=255)
                                     cGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGaca                                                                                                             <  1:411250/148‑1 (MQ=255)
                                      ggCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGacac                                                                                                            <  1:1992785/148‑1 (MQ=255)
                                      ggCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGacac                                                                                                            <  1:325978/148‑1 (MQ=255)
                                        cTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaa                                                                                                          >  2:1608743/1‑148 (MQ=255)
                                          tCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaaca                                                                                                        >  1:462115/1‑148 (MQ=255)
                                          tCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaaca                                                                                                        <  2:8177/148‑1 (MQ=255)
                                          tCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCAGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaaca                                                                                                        <  1:2586263/148‑1 (MQ=255)
                                            aTGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAg                                                                                                      <  2:1448959/148‑1 (MQ=255)
                                              gcgcAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGAt                                                                                                    <  1:882056/148‑1 (MQ=255)
                                                gcagcaCGCGACCCAGCGTGGCGTACTACTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTg                                                                                                  <  1:3010146/148‑1 (MQ=255)
                                                   gcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTa                                                                                               <  2:1943285/148‑1 (MQ=255)
                                                   gcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTa                                                                                               <  2:1467990/148‑1 (MQ=255)
                                                   gcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTGATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTa                                                                                               >  2:2220106/1‑148 (MQ=255)
                                                    caCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTAt                                                                                              <  2:1284170/148‑1 (MQ=255)
                                                    caCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTAt                                                                                              <  1:1707883/148‑1 (MQ=255)
                                                        cgACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCAc                                                                                          >  2:1126767/1‑148 (MQ=255)
                                                           cccAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                       <  2:485345/148‑1 (MQ=255)
                                                           cccAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                       <  1:2570530/148‑1 (MQ=255)
                                                           cccAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                       <  1:1203553/148‑1 (MQ=255)
                                                           cccAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                       >  1:1876809/1‑148 (MQ=255)
                                                           cccAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                       >  1:3309896/1‑148 (MQ=255)
                                                                 gTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCt                                                                                 <  2:1599925/148‑1 (MQ=255)
                                                                    gCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCg                                                                              >  1:2205664/1‑148 (MQ=255)
                                                                    gCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCg                                                                              <  1:1561326/148‑1 (MQ=255)
                                                                     cGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCgg                                                                             >  2:788946/1‑148 (MQ=255)
                                                                     cGTACTTCTCTTCGGTGAAAGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCgg                                                                             >  1:2425535/1‑148 (MQ=255)
                                                                        aCTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGttt                                                                          >  2:3104418/1‑148 (MQ=255)
                                                                         cTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtc                                                                         <  1:892564/148‑1 (MQ=255)
                                                                         cTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtc                                                                         <  2:2536236/148‑1 (MQ=255)
                                                                         cTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtc                                                                         <  2:2960588/148‑1 (MQ=255)
                                                                           tctcTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtcat                                                                       <  2:1526769/148‑1 (MQ=255)
                                                                           tctcTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtcat                                                                       <  1:2625824/148‑1 (MQ=255)
                                                                           tctcTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtcat                                                                       <  1:1568545/148‑1 (MQ=255)
                                                                                tCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAggg                                                                  <  2:2262828/148‑1 (MQ=255)
                                                                                tCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAggg                                                                  >  2:2242582/1‑148 (MQ=255)
                                                                                tCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAggg                                                                  >  2:2242431/1‑148 (MQ=255)
                                                                                tCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAggg                                                                  >  1:3047317/1‑148 (MQ=255)
                                                                                 cGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGt                                                                 <  2:2806833/148‑1 (MQ=255)
                                                                                 cGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGt                                                                 <  1:1791936/148‑1 (MQ=255)
                                                                                  ggTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTc                                                                <  2:2674736/148‑1 (MQ=255)
                                                                                    tGACCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTg                                                              <  2:3011664/148‑1 (MQ=255)
                                                                                            ttATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCgg                                                      <  1:2953666/148‑1 (MQ=255)
                                                                                            ttATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCgg                                                      >  1:2619525/1‑148 (MQ=255)
                                                                                              aTGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAt                                                    >  2:1346076/1‑148 (MQ=255)
                                                                                                 gATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtgtg                                                 <  2:2315329/148‑1 (MQ=255)
                                                                                                  aTGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTgg                                                <  1:793980/148‑1 (MQ=255)
                                                                                                  aTGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTgg                                                <  1:793963/148‑1 (MQ=255)
                                                                                                       gtgtTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGt                                           >  1:2249211/1‑148 (MQ=255)
                                                                                                           tCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAAGTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAgg                                       <  1:1988439/148‑1 (MQ=255)
                                                                                                              cTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTc                                    <  2:2025137/148‑1 (MQ=255)
                                                                                                              cTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTc                                    >  1:17109/1‑148 (MQ=255)
                                                                                                                ttCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCAt                                  >  2:2684029/1‑148 (MQ=255)
                                                                                                                   gCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATggg                               <  1:236179/148‑1 (MQ=255)
                                                                                                                      gATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTg                            >  2:1318765/1‑148 (MQ=255)
                                                                                                                          gTAAAAATACGTGGGATATCCATAATTTAATTACTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATACACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCg                        >  2:2427128/1‑148 (MQ=255)
                                                                                                                            aaaaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgt                      <  2:2745607/148‑1 (MQ=255)
                                                                                                                                     tGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTg             <  2:2619525/148‑1 (MQ=255)
                                                                                                                                     tGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTg             >  1:1938205/1‑148 (MQ=255)
                                                                                                                                      gggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGc            <  1:2533089/148‑1 (MQ=255)
                                                                                                                                       ggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGca           <  2:1750302/148‑1 (MQ=255)
                                                                                                                                       ggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGca           <  2:1750307/148‑1 (MQ=255)
                                                                                                                                         atatCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc         >  1:1513721/1‑148 (MQ=255)
                                                                                                                                          tatCCATAATTTAATTCCTGTTTTTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatca        <  1:2012143/148‑1 (MQ=255)
                                                                                                                                           atCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatcat       >  2:2340490/1‑148 (MQ=255)
                                                                                                                                            tCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCatcatc      >  1:1181178/1‑148 (MQ=255)
                                                                                                                                             ccATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATcc     <  1:2107575/148‑1 (MQ=255)
                                                                                                                                             ccATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATcc     >  1:3179347/1‑148 (MQ=255)
                                                                                                                                                ttattTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcagc  >  2:3257110/1‑148 (MQ=255)
                                                                                                                                                ttattTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcagc  >  2:1619966/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGC  >  CP009273/4522856‑4523147

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: