Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,962,962 T→A 43.8% intergenic (‑626/‑59) rppH ← / → mutH RNA pyrophosphohydrolase/methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,962,9620TA43.8% 100.7 / inf 233intergenic (‑626/‑59)rppH/mutHRNA pyrophosphohydrolase/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (56/75);  new base A (50/52);  total (106/127)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 234 total aligned reads displayed.

CCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTGG  >  CP009273/2962814‑2963108
                                                                                                                                                    |                                                                                                                                                  
cccAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGtt                                                                                                                                                   >  1:2298947/1‑150 (MQ=255)
 ccAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATa                                                                                                                                                  <  2:2394934/150‑1 (MQ=255)
    gCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTacac                                                                                                                                               >  2:77321/1‑150 (MQ=255)
           gtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAat                                                                                                                                        <  1:4747206/150‑1 (MQ=255)
             ggATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAatat                                                                                                                                      >  2:4650890/1‑150 (MQ=255)
              gATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATAtt                                                                                                                                     <  2:1046182/150‑1 (MQ=255)
                tAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCg                                                                                                                                   <  1:4520722/150‑1 (MQ=255)
                tAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCg                                                                                                                                   <  1:3513213/150‑1 (MQ=255)
                 aaCATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCgg                                                                                                                                  <  2:3833133/150‑1 (MQ=255)
                   cATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGca                                                                                                                                <  2:212017/150‑1 (MQ=255)
                   cATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGca                                                                                                                                <  2:2604534/150‑1 (MQ=255)
                    aTGTTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGcac                                                                                                                               <  2:2804837/150‑1 (MQ=255)
                    aTGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGcac                                                                                                                               >  1:4660545/1‑150 (MQ=255)
                    aTGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGcac                                                                                                                               <  2:1510044/150‑1 (MQ=255)
                     tGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGcaca                                                                                                                              <  1:917343/150‑1 (MQ=255)
                       gtgtAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATa                                                                                                                            >  1:3553833/1‑150 (MQ=255)
                       gtgtAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATa                                                                                                                            <  2:2668414/150‑1 (MQ=255)
                          tAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAAtt                                                                                                                         >  1:4688008/1‑150 (MQ=255)
                          tAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAAtt                                                                                                                         >  1:2260849/1‑150 (MQ=255)
                          tAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAAtt                                                                                                                         >  1:2797679/1‑150 (MQ=255)
                           aaGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTg                                                                                                                        <  2:3696472/150‑1 (MQ=255)
                             gATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCt                                                                                                                      <  2:4715004/150‑1 (MQ=255)
                             gATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCt                                                                                                                      <  2:4715769/150‑1 (MQ=255)
                                ccTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCtgtt                                                                                                                   <  2:3928096/150‑1 (MQ=255)
                                 gtGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCtgttt                                                                                                                  >  1:1155078/2‑150 (MQ=255)
                                 cTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgttt                                                                                                                  >  1:821091/1‑150 (MQ=255)
                                  tGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtttg                                                                                                                 <  2:4221016/150‑1 (MQ=255)
                                   gTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCtgtttgt                                                                                                                >  2:3448348/1‑150 (MQ=255)
                                   gTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtttgt                                                                                                                <  2:853799/150‑1 (MQ=255)
                                   gTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtttgt                                                                                                                <  2:1115833/150‑1 (MQ=255)
                                     ttATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCATAGAATTGAACAACGCTTGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTgtttgttt                                                                                                              <  2:2993695/150‑1 (MQ=255)
                                     ttATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTgtttgttt                                                                                                              <  1:2375964/150‑1 (MQ=255)
                                     ttATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTgtttgttt                                                                                                              <  1:2376210/150‑1 (MQ=255)
                                       aTTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGttttt                                                                                                            <  1:4327072/150‑1 (MQ=255)
                                       aTTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGttttt                                                                                                            <  1:3128081/150‑1 (MQ=255)
                                        ttttCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGtttttt                                                                                                           <  2:1178534/150‑1 (MQ=255)
                                          ttCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTaa                                                                                                         >  1:2131966/1‑150 (MQ=255)
                                            cAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATc                                                                                                       >  2:1998895/1‑150 (MQ=255)
                                             aGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCa                                                                                                      >  2:3550664/1‑150 (MQ=255)
                                             aGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCa                                                                                                      <  1:3292111/150‑1 (MQ=255)
                                               tGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAg                                                                                                    >  2:2493756/1‑150 (MQ=255)
                                                 aCCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGt                                                                                                  <  1:2528983/150‑1 (MQ=255)
                                                 aCCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGt                                                                                                  <  2:4709674/150‑1 (MQ=255)
                                                   cAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTAt                                                                                                >  1:2539457/1‑150 (MQ=255)
                                                    aGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATc                                                                                               >  1:3053524/1‑150 (MQ=255)
                                                    aGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATc                                                                                               <  1:3123003/150‑1 (MQ=255)
                                                         tGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCAGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGAc                                                                                          >  2:3356355/1‑150 (MQ=255)
                                                           gAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACAt                                                                                        <  2:4756167/150‑1 (MQ=255)
                                                            aaaaCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATg                                                                                       <  1:1706563/150‑1 (MQ=255)
                                                              aaCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTc                                                                                     >  1:2045066/1‑150 (MQ=255)
                                                              aaCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTc                                                                                     >  1:4442474/1‑150 (MQ=255)
                                                               aCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTcc                                                                                    <  1:2395648/150‑1 (MQ=255)
                                                                 ccGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCa                                                                                  <  2:228737/150‑1 (MQ=255)
                                                                    ttGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAAcc                                                                               >  2:2372045/1‑150 (MQ=255)
                                                                     tGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCACTATAAATCAGAGAATTGAACAACGCATGTGGACAAAGATACACTGCGACTCTTCGGCACCGAAGTGCTGTTTGTTTTGTAATCACGGTATCATGACATGTCCCAACCt                                                                              >  2:2970257/1‑150 (MQ=255)
                                                                     tGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCt                                                                              >  1:4581587/1‑150 (MQ=255)
                                                                             tGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACt                                                                      <  1:3923927/150‑1 (MQ=255)
                                                                               cgcAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGc                                                                    >  2:1629771/1‑150 (MQ=255)
                                                                                gcAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCGGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGct                                                                   >  1:2894260/1‑150 (MQ=255)
                                                                                   aCTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctctc                                                                <  2:1621526/150‑1 (MQ=255)
                                                                                    cTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctctct                                                               >  2:3210258/1‑150 (MQ=255)
                                                                                     tCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCtctctc                                                              >  1:1011263/1‑150 (MQ=255)
                                                                                      cGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcct                                                            <  2:4046044/150‑1 (MQ=255)
                                                                                        aTTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcctc                                                           <  2:1943265/150‑1 (MQ=255)
                                                                                         ttACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                                                          <  1:813508/150‑1 (MQ=255)
                                                                                         ttACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCtcctcc                                                          <  1:263206/150‑1 (MQ=255)
                                                                                            ccGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGa                                                       >  1:1429302/1‑150 (MQ=255)
                                                                                              ggCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaaa                                                     >  1:2321691/1‑150 (MQ=255)
                                                                                              ggCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaaa                                                     >  1:2320571/1‑150 (MQ=255)
                                                                                                cAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACt                                                   >  2:367589/1‑150 (MQ=255)
                                                                                                cAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACt                                                   >  1:3381199/1‑150 (MQ=255)
                                                                                                  aCCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTga                                                 >  1:2672227/1‑150 (MQ=255)
                                                                                                      aaaaaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaagaa                                             <  1:4163288/150‑1 (MQ=255)
                                                                                                      aaaaaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaagaa                                             <  1:2732286/150‑1 (MQ=255)
                                                                                                           gCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgtt                                        <  2:659531/150‑1 (MQ=255)
                                                                                                             aaTATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgt                                      >  1:2215100/1‑150 (MQ=255)
                                                                                                             aaTATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgt                                      >  1:3147255/1‑150 (MQ=255)
                                                                                                                ataAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAg                                   <  1:71043/150‑1 (MQ=255)
                                                                                                                ataAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAg                                   >  2:3753333/1‑150 (MQ=255)
                                                                                                                  aaaTCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAgcg                                 <  2:448281/150‑1 (MQ=255)
                                                                                                                   aaTCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAgcgc                                <  1:1804141/150‑1 (MQ=255)
                                                                                                                       agagAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGc                            >  1:76786/1‑150 (MQ=255)
                                                                                                                       agagAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGc                            <  1:3448440/150‑1 (MQ=255)
                                                                                                                        gagaATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGca                           >  2:146832/1‑150 (MQ=255)
                                                                                                                          gaATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGcaca                         <  1:3994092/150‑1 (MQ=255)
                                                                                                                              tGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCaa                     <  1:1560129/150‑1 (MQ=255)
                                                                                                                              tGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCaa                     >  2:4206307/1‑150 (MQ=255)
                                                                                                                              tGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCaa                     >  2:607447/1‑150 (MQ=255)
                                                                                                                               gAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAAc                    >  1:1364307/1‑150 (MQ=255)
                                                                                                                                aacaacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACt                   >  1:328046/1‑150 (MQ=255)
                                                                                                                                   aacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTc                <  2:1584702/150‑1 (MQ=255)
                                                                                                                                     cGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTg              >  2:3504519/1‑150 (MQ=255)
                                                                                                                                     cGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTg              >  2:1061937/1‑150 (MQ=255)
                                                                                                                                       cATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGt            <  1:4210389/150‑1 (MQ=255)
                                                                                                                                          gtgGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTtat         >  2:4368259/1‑150 (MQ=255)
                                                                                                                                            ggAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATAc       <  2:3108671/150‑1 (MQ=255)
                                                                                                                                             gAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACa      >  1:1054374/1‑150 (MQ=255)
                                                                                                                                                aaaGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTg   >  1:2958529/1‑150 (MQ=255)
                                                                                                                                                aaaGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTg   >  1:4118566/1‑150 (MQ=255)
                                                                                                                                                 aaGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTgg  <  2:3061327/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTGG  >  CP009273/2962814‑2963108

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: