Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,993,566 A→G intergenic (+61/+86) uvrD → / ← yigE DNA‑dependent ATPase I and helicase II/DUF2233 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,993,5660AG99.4% 1202.9 / NA 306intergenic (+61/+86)uvrD/yigEDNA‑dependent ATPase I and helicase II/DUF2233 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (173/131);  minor base T (1/1);  total (174/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

GGTGAGCATAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCC  >  CP009273/3993419‑3993712
                                                                                                                                                   |                                                                                                                                                  
 gTGAGCATAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTg                                                                                                                                                 >  1:181900/1‑150 (MQ=255)
 gTGAGCATAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTg                                                                                                                                                 >  2:1373404/1‑150 (MQ=255)
 gTGAGCATAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTGTGGCCTAAAAAATCATTCAGATTCAATAAGTTg                                                                                                                                                 >  2:359320/1‑150 (MQ=255)
   gAGCATAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCa                                                                                                                                               <  2:459713/150‑1 (MQ=255)
    aGCATAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCaa                                                                                                                                              >  2:2848579/1‑150 (MQ=255)
    aGCATAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCaa                                                                                                                                              >  2:1005346/1‑150 (MQ=255)
       aTAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAATTTGCAACGt                                                                                                                                           >  1:2337923/1‑150 (MQ=255)
       aTAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGt                                                                                                                                           >  1:1326498/1‑150 (MQ=255)
        tAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTc                                                                                                                                          <  2:3536349/150‑1 (MQ=255)
          gCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCAt                                                                                                                                        <  1:1451962/150‑1 (MQ=255)
          gCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCAt                                                                                                                                        <  1:3736651/150‑1 (MQ=255)
           ccGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATg                                                                                                                                       <  2:2689330/150‑1 (MQ=255)
             gTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTa                                                                                                                                     >  2:1465990/1‑150 (MQ=255)
                tGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGc                                                                                                                                  <  2:2203934/150‑1 (MQ=255)
                  cAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCg                                                                                                                                <  2:1871646/150‑1 (MQ=255)
                         cATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAggg                                                                                                                         >  1:2098108/1‑150 (MQ=255)
                         cATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAggg                                                                                                                         >  2:1806518/1‑150 (MQ=255)
                             tCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGtt                                                                                                                     <  1:2779271/150‑1 (MQ=255)
                             tCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGtt                                                                                                                     >  2:2247527/1‑150 (MQ=255)
                              cAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTTACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGttt                                                                                                                    <  1:2645599/150‑1 (MQ=255)
                               aGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTa                                                                                                                   >  1:1653629/1‑150 (MQ=255)
                                   ccAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgcc                                                                                                               >  2:857607/1‑150 (MQ=255)
                                     aGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACTTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTAcgccgc                                                                                                             >  1:3228951/1‑150 (MQ=255)
                                        gTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATc                                                                                                          >  1:1589179/1‑150 (MQ=255)
                                          aTTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCg                                                                                                        >  2:2577979/1‑150 (MQ=255)
                                            tAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGc                                                                                                      >  1:3109374/1‑150 (MQ=255)
                                               aTGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATc                                                                                                   >  1:299151/1‑150 (MQ=255)
                                                tGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCt                                                                                                  >  1:2263118/1‑150 (MQ=255)
                                                tGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCt                                                                                                  >  1:3611157/1‑150 (MQ=255)
                                                   cTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgcg                                                                                               >  2:123439/1‑150 (MQ=255)
                                                   cTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgcg                                                                                               <  2:69210/150‑1 (MQ=255)
                                                   cTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgcg                                                                                               >  1:3620165/1‑150 (MQ=255)
                                                    tggtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgcgc                                                                                              >  1:3247592/1‑150 (MQ=255)
                                                    tggtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTgcgc                                                                                              >  1:2233750/1‑150 (MQ=255)
                                                      gtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCa                                                                                            >  1:3069517/1‑150 (MQ=255)
                                                      gtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCa                                                                                            <  2:1586872/150‑1 (MQ=255)
                                                           ggcATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTc                                                                                       >  2:1640363/1‑150 (MQ=255)
                                                              aTACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTGAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGt                                                                                    <  2:1296994/150‑1 (MQ=255)
                                                              aTACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGt                                                                                    >  2:3449026/1‑150 (MQ=255)
                                                              aTACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGt                                                                                    >  1:3668370/1‑150 (MQ=255)
                                                                  gCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCt                                                                                <  2:1976253/150‑1 (MQ=255)
                                                                     cGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGac                                                                             <  1:3465572/150‑1 (MQ=255)
                                                                      ggCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGaca                                                                            <  1:1508623/150‑1 (MQ=255)
                                                                         tGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACaa                                                                         >  2:813518/1‑150 (MQ=255)
                                                                            aGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAAc                                                                      >  1:879946/1‑150 (MQ=255)
                                                                             gTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACt                                                                     >  2:3429842/1‑150 (MQ=255)
                                                                             gTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACt                                                                     <  2:2171957/150‑1 (MQ=255)
                                                                              tCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTa                                                                    <  1:215553/150‑1 (MQ=255)
                                                                                ggTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATc                                                                  <  2:239624/150‑1 (MQ=255)
                                                                                  tgtAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGt                                                                <  2:2511530/150‑1 (MQ=255)
                                                                                    tAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGttt                                                              >  1:3718838/1‑150 (MQ=255)
                                                                                    tAACGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGttt                                                              >  2:59398/1‑150 (MQ=255)
                                                                                     aaCGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtttt                                                             <  2:1383921/150‑1 (MQ=255)
                                                                                     aaCGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtttt                                                             >  1:1162592/1‑150 (MQ=255)
                                                                                     aaCGTTGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGtttt                                                             <  1:2578614/150‑1 (MQ=255)
                                                                                          tGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCt                                                        <  1:813518/150‑1 (MQ=255)
                                                                                          tGCCGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCt                                                        >  1:682022/1‑150 (MQ=255)
                                                                                            ccGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCttt                                                      >  1:3571476/1‑150 (MQ=255)
                                                                                            ccGGATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCttt                                                      >  2:1708578/1‑150 (MQ=255)
                                                                                              ggATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTcc                                                    >  1:3682058/1‑150 (MQ=255)
                                                                                               gATGCGGTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCg                                                   >  2:2982798/1‑150 (MQ=255)
                                                                                                     gTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGAc                                             <  2:3232891/150‑1 (MQ=255)
                                                                                                     gTGCTGCGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGAc                                             <  1:1868270/150‑1 (MQ=255)
                                                                                                         tgcGCACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGaa                                         >  2:3500426/1‑150 (MQ=255)
                                                                                                           cgcACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAAt                                       <  1:2247527/150‑1 (MQ=255)
                                                                                                            gcACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGACATc                                      >  1:47340/1‑150 (MQ=255)
                                                                                                             cACCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCa                                     >  1:909719/1‑150 (MQ=255)
                                                                                                              aCCTTGTTTGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCAt                                    >  2:2142165/1‑150 (MQ=255)
                                                                                                                 ttgtttGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGt                                 <  1:123439/147‑1 (MQ=255)
                                                                                                                  tgtttGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTg                                >  2:413510/3‑150 (MQ=255)
                                                                                                                    tttGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGAc                              <  1:3285446/150‑1 (MQ=255)
                                                                                                                    tttGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGAc                              <  1:1385728/150‑1 (MQ=255)
                                                                                                                      tGGCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACaa                            >  1:42044/1‑150 (MQ=255)
                                                                                                                        gCCTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACaaaa                          >  1:599880/1‑150 (MQ=255)
                                                                                                                         ccTAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAg                         <  2:2070469/150‑1 (MQ=255)
                                                                                                                           tAAAAAATCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGa                       <  2:3069742/150‑1 (MQ=255)
                                                                                                                            aaaaaaTCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGAt                      >  1:2355878/1‑150 (MQ=255)
                                                                                                                             aaaaaTCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATa                     <  2:3611157/150‑1 (MQ=255)
                                                                                                                              aaaaTCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAg                    >  1:2843433/1‑150 (MQ=255)
                                                                                                                              aaaaTCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAg                    >  1:3262337/1‑150 (MQ=255)
                                                                                                                              aaaaTCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAg                    >  2:18824/1‑150 (MQ=255)
                                                                                                                               aaaTCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGc                   <  1:3504127/150‑1 (MQ=255)
                                                                                                                               aaaTCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGc                   >  1:1030691/1‑150 (MQ=255)
                                                                                                                                aaTCATTCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCg                  <  1:2493725/150‑1 (MQ=255)
                                                                                                                                  tCATTGAGATTAAATAATTTGCAACGTCATGTAGGACGGATAGGGCGTTTACGCCGCATCCGTCATCTGCTCCATATTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCTACGGATATCATGGTGACATATGGATAGCGtt                <  2:2337923/150‑1 (MQ=255)
                                                                                                                                     ttCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGcc             <  2:181900/150‑1 (MQ=255)
                                                                                                                                      tCAGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCa            <  1:2950868/150‑1 (MQ=255)
                                                                                                                                       taGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCAc           >  2:3714995/2‑150 (MQ=255)
                                                                                                                                        aGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACg          <  2:3620165/150‑1 (MQ=255)
                                                                                                                                        aGATTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACg          <  2:888905/150‑1 (MQ=255)
                                                                                                                                          aTTCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACggg        <  1:576630/150‑1 (MQ=255)
                                                                                                                                           ttCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGa       <  2:564060/150‑1 (MQ=255)
                                                                                                                                            tCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGAt      <  1:1708578/150‑1 (MQ=255)
                                                                                                                                            tCAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGAt      >  2:3141996/1‑150 (MQ=255)
                                                                                                                                             cAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGAtt     <  1:1936024/150‑1 (MQ=255)
                                                                                                                                             cAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGAtt     <  2:3247592/150‑1 (MQ=255)
                                                                                                                                             cAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGAtt     <  2:3201339/150‑1 (MQ=255)
                                                                                                                                             cAATAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGAtt     <  1:857607/150‑1 (MQ=255)
                                                                                                                                               aTAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGc   >  1:1752494/1‑150 (MQ=255)
                                                                                                                                               aTAAGTTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGc   >  2:154553/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGTGAGCATAGCCGTTTGCAGGTGGCATTTCAGGGCCAGGGTATTAAATGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCC  >  CP009273/3993419‑3993712

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: