Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,213,666 C→A 6.4% intergenic (‑110/‑90) iclR ← / → metH transcriptional repressor/homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,213,6660CA 6.4% 829.7 / 27.0 239intergenic (‑110/‑90)iclR/metHtranscriptional repressor/homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base C (108/114);  new base A (12/4);  total (121/118)
Fisher's exact test for biased strand distribution p-value = 6.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGG  >  CP009273/4213521‑4213813
                                                                                                                                                 |                                                                                                                                                   
ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtgg                                                                                                                                                 <  1:4868211/150‑1 (MQ=255)
ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtgg                                                                                                                                                 <  1:2291775/150‑1 (MQ=255)
ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtgg                                                                                                                                                 <  1:2219041/150‑1 (MQ=255)
  ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGAgtggtg                                                                                                                                               >  2:1912929/1‑150 (MQ=255)
   tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGtggtgg                                                                                                                                              <  2:3285709/150‑1 (MQ=255)
   tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGtggtgg                                                                                                                                              <  1:51853/150‑1 (MQ=255)
     tGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTc                                                                                                                                            >  2:3890335/1‑150 (MQ=255)
      gCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCg                                                                                                                                           <  2:1490939/150‑1 (MQ=255)
          cgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCttt                                                                                                                                       <  1:1927780/150‑1 (MQ=255)
           gTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCtttt                                                                                                                                      >  1:944245/1‑150 (MQ=255)
              tCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTGTTAcc                                                                                                                                   >  1:1379124/1‑150 (MQ=255)
               cgcgGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACca                                                                                                                                  <  2:4799471/150‑1 (MQ=255)
                gcgGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcac                                                                                                                                 >  1:4717733/1‑150 (MQ=255)
                  gggAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAg                                                                                                                               >  2:759305/1‑150 (MQ=255)
                    gAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAt                                                                                                                             <  2:4273680/150‑1 (MQ=255)
                       tGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCg                                                                                                                          >  2:2067790/1‑150 (MQ=255)
                        gggTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGt                                                                                                                         <  1:3957146/150‑1 (MQ=255)
                          gTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGttt                                                                                                                       <  1:5126067/150‑1 (MQ=255)
                          gTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGttt                                                                                                                       <  2:560766/150‑1 (MQ=255)
                           tGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCAAAGATGCGTTTa                                                                                                                      >  1:440444/1‑150 (MQ=255)
                            gCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTAt                                                                                                                     <  2:3518613/150‑1 (MQ=255)
                            gCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTAt                                                                                                                     <  2:918868/150‑1 (MQ=255)
                            gCGACCATGACAGTATCTTTTTTCTGTATCGTGTAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTAt                                                                                                                     <  2:5134086/150‑1 (MQ=255)
                             cGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATg                                                                                                                    >  1:938579/1‑150 (MQ=255)
                             cGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATg                                                                                                                    <  1:324982/150‑1 (MQ=255)
                               aCCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGcc                                                                                                                  >  2:2098498/1‑150 (MQ=255)
                               aCCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGcc                                                                                                                  >  2:835269/1‑150 (MQ=255)
                                  aTGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGt                                                                                                               >  2:2607791/1‑150 (MQ=255)
                                    gACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATTCCAGTAt                                                                                                             >  1:3995416/1‑150 (MQ=255)
                                    gACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                             >  2:4410723/1‑150 (MQ=255)
                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATg                                                                                                            <  1:3823527/150‑1 (MQ=255)
                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATg                                                                                                            <  1:1226912/150‑1 (MQ=255)
                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATg                                                                                                            <  1:1057741/150‑1 (MQ=255)
                                     aCAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAAATGATCCGATGAGTTAATTTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATg                                                                                                            >  1:4248425/1‑150 (MQ=255)
                                      cAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATgg                                                                                                           <  2:4243276/150‑1 (MQ=255)
                                             ttttttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgttg                                                                                                    >  2:4423032/1‑150 (MQ=255)
                                                tttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAt                                                                                                 >  2:3700231/1‑150 (MQ=255)
                                                tttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAt                                                                                                 <  1:4895765/150‑1 (MQ=255)
                                                 ttCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAtt                                                                                                <  2:1559218/150‑1 (MQ=255)
                                                   cTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAtttt                                                                                              >  1:2786955/1‑150 (MQ=255)
                                                   cTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAtttt                                                                                              <  2:3871428/150‑1 (MQ=255)
                                                     gTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTtta                                                                                            <  2:1308658/150‑1 (MQ=255)
                                                      tATCGTGGAAATCATTTTCATTTTTATCGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttat                                                                                           <  2:2645764/150‑1 (MQ=255)
                                                      tATCGTGGAAATCATTTTCATTTTTATCGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttat                                                                                           <  2:27326/150‑1 (MQ=255)
                                                       aTCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttatt                                                                                          <  1:4451530/150‑1 (MQ=255)
                                                        tCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttatta                                                                                         >  1:972802/1‑150 (MQ=255)
                                                          gTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTaaa                                                                                       >  1:2612326/1‑150 (MQ=255)
                                                            ggAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGAGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATc                                                                                     >  2:1098597/1‑150 (MQ=255)
                                                             gAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCt                                                                                    <  2:4501127/150‑1 (MQ=255)
                                                             gAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCt                                                                                    <  1:2205831/150‑1 (MQ=255)
                                                               aaTCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTgg                                                                                  <  1:4097200/150‑1 (MQ=255)
                                                               aaTCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTgg                                                                                  >  1:4604653/1‑150 (MQ=255)
                                                                 tcattttcatttTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGt                                                                                <  1:1130663/150‑1 (MQ=255)
                                                                  cattttcattttTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGtt                                                                               <  1:2351175/150‑1 (MQ=255)
                                                                   attttcattttTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTTTGCCAGTATGGTTTGTTGAATTTTTATTAAATATGGGTTg                                                                              >  1:2200109/1‑150 (MQ=255)
                                                                    ttttcattttTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGa                                                                             <  1:4433484/150‑1 (MQ=255)
                                                                    ttttcattttTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGAGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGa                                                                             >  2:285575/1‑150 (MQ=255)
                                                                      ttcattttTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGc                                                                           <  2:4517782/150‑1 (MQ=255)
                                                                       tcattttTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCg                                                                          <  2:2255942/150‑1 (MQ=255)
                                                                          tttttATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGAGGTGGTCGCTTTTACCACAGATGCGGTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATATGGGTTGAGCgtgt                                                                       >  1:3052037/1‑150 (MQ=255)
                                                                           ttttATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTc                                                                      <  2:685318/150‑1 (MQ=255)
                                                                             ttATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCgg                                                                    <  2:4636839/150‑1 (MQ=255)
                                                                              tATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCggg                                                                   >  2:2511285/1‑150 (MQ=255)
                                                                               aTTTTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTTGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGa                                                                  <  1:736781/150‑1 (MQ=255)
                                                                               aTTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGa                                                                  >  2:1492245/1‑150 (MQ=255)
                                                                               aTTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGa                                                                  <  1:1850326/150‑1 (MQ=255)
                                                                               aTTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGa                                                                  >  1:4962410/1‑150 (MQ=255)
                                                                                  gTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCa                                                               >  1:5026002/1‑150 (MQ=255)
                                                                                    tAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAg                                                             >  1:773019/1‑150 (MQ=255)
                                                                                     aGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAgt                                                            >  2:331671/1‑150 (MQ=255)
                                                                                      gCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAgtg                                                           >  2:1607071/1‑150 (MQ=255)
                                                                                      gCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAgtg                                                           <  1:2926297/150‑1 (MQ=255)
                                                                                        tAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGtgtg                                                         <  2:123269/150‑1 (MQ=255)
                                                                                        tAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAATTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGtgtg                                                         >  1:2250952/1‑150 (MQ=255)
                                                                                        tAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTATGCCCATTTGGGTTGTTTAAATTTTAATAAATATGGGTTTAGCGTGTAGGGAGCCAGGGTGa                                                        >  2:3118996/1‑150 (MQ=255)
                                                                                            gCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGagca                                                     >  2:3291104/1‑150 (MQ=255)
                                                                                               aTAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCCCAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAgcagca                                                  >  1:3509983/1‑150 (MQ=255)
                                                                                                   ttACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                              <  2:3330229/150‑1 (MQ=255)
                                                                                                    tACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTg                                             <  1:3064431/150‑1 (MQ=255)
                                                                                                      ctgaactgaTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGa                                           <  2:2517456/150‑1 (MQ=255)
                                                                                                            tgaTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACt                                     >  1:4747934/1‑150 (MQ=255)
                                                                                                            tgaTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGGTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACt                                     >  1:2170685/1‑150 (MQ=255)
                                                                                                                  gATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtgc                               <  1:1669825/150‑1 (MQ=255)
                                                                                                                  gATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtgc                               <  1:920734/150‑1 (MQ=255)
                                                                                                                       gTTAATGTTGAACAAATCTCATGTTGAGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt                          <  2:4010495/150‑1 (MQ=255)
                                                                                                                         tAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                        <  1:1486349/150‑1 (MQ=255)
                                                                                                                         tAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGGATGGTGTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGttc                        >  2:305487/1‑149 (MQ=255)
                                                                                                                          aaTGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaa                       <  1:878346/150‑1 (MQ=255)
                                                                                                                             gTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATg                    >  2:435885/1‑150 (MQ=255)
                                                                                                                                 aaCAAATCTCATGTTGAGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCCGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg                >  2:649378/1‑150 (MQ=255)
                                                                                                                                  aCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGt               <  2:2786955/150‑1 (MQ=255)
                                                                                                                                   cAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTa              <  2:2173454/150‑1 (MQ=255)
                                                                                                                                   cAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTa              >  2:4679251/1‑150 (MQ=255)
                                                                                                                                   cAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCACTTAAATGAACGTa              <  1:41265/150‑1 (MQ=255)
                                                                                                                                     aaTCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt            >  2:3611080/1‑150 (MQ=255)
                                                                                                                                           aTGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTg      <  1:3039269/150‑1 (MQ=255)
                                                                                                                                            tGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc     <  1:2442694/150‑1 (MQ=255)
                                                                                                                                            tGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGc     >  1:4932502/1‑150 (MQ=255)
                                                                                                                                              ttGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTg   <  1:2883015/150‑1 (MQ=255)
                                                                                                                                               tGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTgg  <  2:3247847/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTGG  >  CP009273/4213521‑4213813

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: