Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,962,962 T→A 91.6% 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,9620TA91.6% 961.7 / 58.9 298intergenic (‑626/‑59)rppH/mutHRNA pyrophosphohydrolase/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (16/9);  new base A (149/124);  total (165/133)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

CCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTGG  >  CP009273/2962814‑2963108
                                                                                                                                                    |                                                                                                                                                  
 ccAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATa                                                                                                                                                  >  2:2203787/1‑150 (MQ=255)
   aGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATaca                                                                                                                                                >  2:381835/1‑150 (MQ=255)
    gCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGTGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATacac                                                                                                                                               <  1:2643948/150‑1 (MQ=255)
    gCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATacac                                                                                                                                               >  2:393621/1‑150 (MQ=255)
         cTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGa                                                                                                                                          >  2:2427499/1‑150 (MQ=255)
           gtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAat                                                                                                                                        <  1:1561732/150‑1 (MQ=255)
           gtgGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAat                                                                                                                                        <  2:3477089/150‑1 (MQ=255)
              gATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATAtt                                                                                                                                     <  2:4408306/150‑1 (MQ=255)
                tAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCg                                                                                                                                   <  1:786662/150‑1 (MQ=255)
                    aTGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGcac                                                                                                                               <  2:2499702/150‑1 (MQ=255)
                      ggTGTGAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACAt                                                                                                                             <  1:1710045/150‑1 (MQ=255)
                           aaGCTCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGATTTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTg                                                                                                                        <  2:134720/150‑1 (MQ=255)
                           aaGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTg                                                                                                                        >  1:979238/1‑150 (MQ=255)
                           aaGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTg                                                                                                                        >  2:144879/1‑150 (MQ=255)
                            aGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGc                                                                                                                       >  1:937957/1‑150 (MQ=255)
                             gATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCt                                                                                                                      <  2:229668/150‑1 (MQ=255)
                              aTCCTGTTTATTTTCAGTGACCTGATTTGGAAAACCAGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtg                                                                                                                     <  2:4168074/150‑1 (MQ=255)
                                ccTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtt                                                                                                                   <  2:1891722/150‑1 (MQ=255)
                                 cTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgttt                                                                                                                  <  1:1230113/150‑1 (MQ=255)
                                  tGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCtgtttg                                                                                                                 <  1:863255/150‑1 (MQ=255)
                                   gTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCtgtttgt                                                                                                                >  2:3649686/1‑150 (MQ=255)
                                     ttATTTTCAGTGACCTGATTTGGTAAACCCGTTGCAGTGTTGCGCTACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACTACGCATGTGGTAAAAGATACTCTGCGAATATTCGGCACATAATTGCTgtttgttt                                                                                                              <  2:2666465/150‑1 (MQ=255)
                                     ttATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTgtttgttt                                                                                                              >  1:4551562/1‑150 (MQ=255)
                                          ttCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGACTTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTaa                                                                                                         >  2:2726189/1‑150 (MQ=255)
                                             aGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTACTGTTTGTTTTTTAATCa                                                                                                      <  2:2762813/150‑1 (MQ=255)
                                                gACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAgg                                                                                                   >  1:48380/1‑150 (MQ=255)
                                                gACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAgg                                                                                                   >  1:1602718/1‑150 (MQ=255)
                                                  ccAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCAATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTa                                                                                                 <  2:4461283/150‑1 (MQ=255)
                                                   cAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTAt                                                                                                <  1:3040554/150‑1 (MQ=255)
                                                   cAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTAt                                                                                                >  2:1048241/1‑150 (MQ=255)
                                                     gATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATACATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCa                                                                                              <  1:1139686/150‑1 (MQ=255)
                                                      aTTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCAt                                                                                             <  1:2124915/150‑1 (MQ=255)
                                                        ttGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGa                                                                                           >  2:3870037/1‑150 (MQ=255)
                                                         tGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGAc                                                                                          <  1:3919949/150‑1 (MQ=255)
                                                          ggAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACa                                                                                         >  1:596841/1‑150 (MQ=255)
                                                           gAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACAt                                                                                        >  1:509547/1‑150 (MQ=255)
                                                            aaaaCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCCTGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATg                                                                                       >  2:1958613/1‑150 (MQ=255)
                                                              aaCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTc                                                                                     >  1:1848197/1‑150 (MQ=255)
                                                              aaCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTc                                                                                     <  2:1486940/150‑1 (MQ=255)
                                                               aCCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTcc                                                                                    >  2:3354700/1‑150 (MQ=255)
                                                                 ccGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCa                                                                                  >  1:3991025/1‑150 (MQ=255)
                                                                   gTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAAc                                                                                <  2:2841700/150‑1 (MQ=255)
                                                                    ttGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAAcc                                                                               <  1:1148544/150‑1 (MQ=255)
                                                                                 cAACTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctc                                                                  >  2:3279727/1‑150 (MQ=255)
                                                                                  aaCTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctct                                                                 <  2:1171782/150‑1 (MQ=255)
                                                                                  aaCTCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGctct                                                                 >  1:126532/1‑150 (MQ=255)
                                                                                     tCGATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCtctctc                                                              >  2:2147168/1‑150 (MQ=255)
                                                                                       gATTACCTGCAACCTTAAAAGCAATATAAATCAGAGAATTGATCATCGCATGTGGAAAAAGATACACTGCGTATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcct                                                            <  2:4635165/150‑1 (MQ=255)
                                                                                       gATTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcct                                                            >  1:3818238/1‑150 (MQ=255)
                                                                                        aTTACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTctcctc                                                           >  1:2187697/1‑150 (MQ=255)
                                                                                          tACCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTccc                                                         >  2:3244999/1‑150 (MQ=255)
                                                                                           aCCGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCg                                                        >  1:3764926/1‑150 (MQ=255)
                                                                                            ccGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGa                                                       >  1:3871680/1‑150 (MQ=255)
                                                                                            ccGGCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGa                                                       <  1:2862236/150‑1 (MQ=255)
                                                                                              ggCAACCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGaaa                                                     <  2:1656963/150‑1 (MQ=255)
                                                                                                 aaCCTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTg                                                  <  2:3078024/150‑1 (MQ=255)
                                                                                                   ccTAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaa                                                <  2:2646503/150‑1 (MQ=255)
                                                                                                     tAAAAAGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTgaaga                                              >  2:605754/1‑150 (MQ=255)
                                                                                                        aaaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACa                                           >  1:404029/1‑150 (MQ=255)
                                                                                                        aaaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACa                                           >  1:2324693/1‑150 (MQ=255)
                                                                                                         aaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAg                                          <  1:1349050/150‑1 (MQ=255)
                                                                                                         aaGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAg                                          <  1:2129133/150‑1 (MQ=255)
                                                                                                          aGCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgt                                         >  1:707879/1‑150 (MQ=255)
                                                                                                           gCAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgtt                                        >  2:1295388/1‑150 (MQ=255)
                                                                                                            cAATATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttg                                       >  2:2229334/1‑150 (MQ=255)
                                                                                                             aaTATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgt                                      <  1:2245017/150‑1 (MQ=255)
                                                                                                             aaTATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgt                                      <  1:1908895/150‑1 (MQ=255)
                                                                                                             aaTATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgt                                      >  2:2577028/1‑150 (MQ=255)
                                                                                                             aaTATAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgt                                      <  1:2519528/150‑1 (MQ=255)
                                                                                                              atatAAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAgttgtt                                     <  1:3549108/150‑1 (MQ=255)
                                                                                                               tataAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTa                                    >  1:195953/1‑150 (MQ=255)
                                                                                                               tataAATCAGAGAATTGAACAACGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTa                                    >  1:4453069/1‑150 (MQ=255)
                                                                                                               tataAATCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTa                                    >  2:2897327/1‑150 (MQ=255)
                                                                                                                   aaTCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAgcgc                                >  1:603157/1‑150 (MQ=255)
                                                                                                                    aTCAGAGAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCa                               <  1:2770472/150‑1 (MQ=255)
                                                                                                                       agagAATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGc                            >  1:4574347/1‑150 (MQ=255)
                                                                                                                        gagaATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGca                           >  1:247788/1‑150 (MQ=255)
                                                                                                                        gagaATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGca                           >  1:3442121/1‑150 (MQ=255)
                                                                                                                          gaATTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGcaca                         <  1:3115483/150‑1 (MQ=255)
                                                                                                                            aTTGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGc                       >  2:3744234/1‑150 (MQ=255)
                                                                                                                              tGAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCaa                     <  2:4314250/150‑1 (MQ=255)
                                                                                                                               gAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAAc                    >  2:512278/1‑150 (MQ=255)
                                                                                                                               gAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAAc                    >  2:3389900/1‑150 (MQ=255)
                                                                                                                               gAACAACGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAAc                    <  2:4220383/150‑1 (MQ=255)
                                                                                                                                 acaacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACtt                  <  2:2687971/150‑1 (MQ=255)
                                                                                                                                 acaacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACtt                  <  2:1770813/150‑1 (MQ=255)
                                                                                                                                  caacGCATGTGGAAAAAGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACttt                 >  2:2075182/1‑150 (MQ=255)
                                                                                                                                  caacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACttt                 >  2:861834/1‑150 (MQ=255)
                                                                                                                                   aacGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTc                <  1:2203787/150‑1 (MQ=255)
                                                                                                                                     cGCATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTg              >  2:1556425/1‑150 (MQ=255)
                                                                                                                                       cATGTGGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGt            <  2:2265849/150‑1 (MQ=255)
                                                                                                                                         tgtgGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGttc          >  1:831523/1‑149 (MQ=255)
                                                                                                                                           tgGAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTtata        >  1:242172/1‑150 (MQ=255)
                                                                                                                                            ggAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATAc       >  1:4540841/1‑150 (MQ=255)
                                                                                                                                             gAAAAAGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACa      <  1:3322869/150‑1 (MQ=255)
                                                                                                                                               aaaaGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACAtt    <  2:4249018/150‑1 (MQ=255)
                                                                                                                                                 aaGTTACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTgg  <  1:1947630/150‑1 (MQ=255)
                                                                                                                                                 aaGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTgg  <  1:88222/150‑1 (MQ=255)
                                                                                                                                                 aaGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTgg  >  2:4527491/1‑150 (MQ=255)
                                                                                                                                                 aaGATACACTGCGAATATTCGGCACATAATTGCTGTTTGTTTTTTAATCAAGGTATCATGACATGTCCCAACCTCGCCCACTGCTCTCTCCTCCCGAAACTGAAGAACAGTTGTTAGCGCAAGCACAGCAACTTTCTGGTTATACATTgg  >  2:2395468/1‑150 (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: