Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,116,964 A→G 92.2% intergenic (‑21/+53) dinI ← / ← pyrC DNA damage‑inducible protein I/dihydro‑orotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,116,9640AG92.2% 764.8 / 39.5 230intergenic (‑21/+53)dinI/pyrCDNA damage‑inducible protein I/dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base A (8/10);  new base G (107/104);  total (116/114)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 233 total aligned reads displayed.

TACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCT  >  CP009273/1116818‑1117110
                                                                                                                                                  |                                                                                                                                                  
tACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTg                                                                                                                                                 >  1:276464/1‑150 (MQ=255)
  ccGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGt                                                                                                                                               <  1:1944093/150‑1 (MQ=255)
    gATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGtta                                                                                                                                             >  1:3674784/1‑150 (MQ=255)
     aTACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttat                                                                                                                                            >  1:3889766/1‑150 (MQ=255)
      tACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttatt                                                                                                                                           <  1:3817596/150‑1 (MQ=255)
      tACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttatt                                                                                                                                           <  1:3817216/150‑1 (MQ=255)
            gCCTTCATTATCAGGAAACGCATACTGAATCCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTAGGGTGACTTCAATTCGCATAAGAGCCCCCTGTTGACTATGCTGGTTATTTATAcc                                                                                                                                     <  1:784535/150‑2 (MQ=255)
            gCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACa                                                                                                                                     <  1:3486944/150‑1 (MQ=255)
                tCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGtc                                                                                                                                 >  2:2920391/1‑149 (MQ=255)
                   ttATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTaaaa                                                                                                                              <  1:2698376/150‑1 (MQ=255)
                     aTCAGGAAACGCATTCTGAATTCGGCGGGAATGTTCGCCAGCCAGGGCGTCAATAGCCCAAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATa                                                                                                                            <  2:258673/150‑1 (MQ=255)
                         ggAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCt                                                                                                                        >  1:1115804/1‑150 (MQ=255)
                         ggAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGATTATTTATACAGGTAAAATAACCt                                                                                                                        >  2:617921/1‑150 (MQ=255)
                          gAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTa                                                                                                                       <  1:2777734/150‑1 (MQ=255)
                          gAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTa                                                                                                                       >  1:1489923/1‑150 (MQ=255)
                           aaaCGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTaa                                                                                                                      >  2:300855/1‑150 (MQ=255)
                              cGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGa                                                                                                                   <  2:1651472/150‑1 (MQ=255)
                                  tACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGacaac                                                                                                               <  1:80531/150‑1 (MQ=255)
                                    cTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAg                                                                                                             <  1:2987681/150‑1 (MQ=255)
                                     tGAATACGGCGGGAAAGTTCGCAAGCCAGGGCGTCAATAGCCCCAGCAGGCAATGGAGAAGTTTTATCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGacaacatg                                                                                                            >  2:3627189/1‑148 (MQ=255)
                                      gAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGa                                                                                                           >  2:3906740/1‑150 (MQ=255)
                                        agaCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAg                                                                                                         <  1:1940007/148‑1 (MQ=255)
                                        aTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAg                                                                                                         <  1:2380671/150‑1 (MQ=255)
                                         tACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGATTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                        <  1:2901003/150‑1 (MQ=255)
                                          aCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACTACAGGAAGCt                                                                                                       >  2:296358/1‑150 (MQ=255)
                                                ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAtt                                                                                                 <  2:4289539/150‑1 (MQ=255)
                                                 gAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttt                                                                                                <  2:1394983/150‑1 (MQ=255)
                                                   aaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttttt                                                                                              <  1:457673/150‑1 (MQ=255)
                                                   aaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttttt                                                                                              <  1:457643/150‑1 (MQ=255)
                                                   aaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttttt                                                                                              >  2:2160407/1‑150 (MQ=255)
                                                   aaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttttt                                                                                              <  1:84639/150‑1 (MQ=255)
                                                    aGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTa                                                                                             >  1:1986201/1‑150 (MQ=255)
                                                      ttCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAtt                                                                                           <  2:1120165/150‑1 (MQ=255)
                                                         gccagccaGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGtt                                                                                        >  1:3897275/1‑150 (MQ=255)
                                                           cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                      <  2:1169602/150‑1 (MQ=255)
                                                            agccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGATAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTaa                                                                                     >  2:4032218/1‑150 (MQ=255)
                                                              ccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACg                                                                                   <  1:1296485/150‑1 (MQ=255)
                                                              ccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACg                                                                                   <  1:312370/150‑1 (MQ=255)
                                                              ccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATATTTATTGTTTAAct                                                                                   >  2:2672851/1‑149 (MQ=255)
                                                               cagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACgg                                                                                  <  2:1180962/150‑1 (MQ=255)
                                                                agGGCGTCAATAGCCCCAGCTGGCACTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGa                                                                                 >  2:619925/1‑150 (MQ=255)
                                                                 gggCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGAc                                                                                <  2:848772/150‑1 (MQ=255)
                                                                    cGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAg                                                                             <  1:1847369/150‑1 (MQ=255)
                                                                      tCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCg                                                                           >  1:2474312/1‑150 (MQ=255)
                                                                                  gCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcccc                                                               <  1:4324368/150‑1 (MQ=255)
                                                                                   cTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCg                                                              >  2:2881550/1‑150 (MQ=255)
                                                                                      gCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCg                                                           <  2:1438511/150‑1 (MQ=255)
                                                                                       cAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGa                                                          >  2:351070/1‑150 (MQ=255)
                                                                                        aaTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAg                                                         >  1:1970519/1‑150 (MQ=255)
                                                                                        aaTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAg                                                         <  1:2855532/150‑1 (MQ=255)
                                                                                         aTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                        >  2:2790624/1‑150 (MQ=255)
                                                                                           ggAGAAGTTTTCGCTATGGTGACTTCAATTAGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACGGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGaa                                                      >  1:882112/1‑150 (MQ=255)
                                                                                            gagaAGTTTTCGCTATGGTGACTTCAATTCGCATTATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAAt                                                     <  1:560850/150‑1 (MQ=255)
                                                                                             agaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATg                                                    >  2:3325020/1‑150 (MQ=255)
                                                                                               aaGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGc                                                  <  1:4148335/150‑1 (MQ=255)
                                                                                                aGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCa                                                 >  2:3589647/1‑150 (MQ=255)
                                                                                                aGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCa                                                 >  2:3589637/1‑150 (MQ=255)
                                                                                                  ttttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCAcc                                               >  1:1241526/1‑150 (MQ=255)
                                                                                                    ttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAg                                             <  2:3816847/150‑1 (MQ=255)
                                                                                                          aTGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtc                                       >  1:4323527/1‑150 (MQ=255)
                                                                                                            ggTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcat                                     <  1:2161538/150‑1 (MQ=255)
                                                                                                              tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatca                                   <  1:1050228/150‑1 (MQ=255)
                                                                                                              tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatca                                   <  1:752609/150‑1 (MQ=255)
                                                                                                               gACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcag                                  >  1:1359678/1‑150 (MQ=255)
                                                                                                                aCTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcagt                                 >  2:483936/1‑150 (MQ=255)
                                                                                                                  ttCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGACTGGCACCAGCGTGTCAtcagtca                               >  2:2362423/1‑150 (MQ=255)
                                                                                                                        tCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGa                         >  1:3314406/1‑150 (MQ=255)
                                                                                                                           cATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGc                      >  2:663906/1‑150 (MQ=255)
                                                                                                                           cATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGc                      >  2:3543286/1‑150 (MQ=255)
                                                                                                                           cATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGAGATGc                      >  2:3063610/1‑150 (MQ=255)
                                                                                                                           cATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGAGCGATGc                      >  2:4424559/1‑150 (MQ=255)
                                                                                                                            ataataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCt                     <  1:2526570/150‑1 (MQ=255)
                                                                                                                            ataataGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCt                     <  1:4263066/150‑1 (MQ=255)
                                                                                                                             taataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCtt                    <  2:1489923/150‑1 (MQ=255)
                                                                                                                             taataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCtt                    <  2:1311256/150‑1 (MQ=255)
                                                                                                                               ataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc                  >  1:1360985/1‑150 (MQ=255)
                                                                                                                                taGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCa                 >  1:2967697/1‑150 (MQ=255)
                                                                                                                                taGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCa                 <  1:52683/150‑1 (MQ=255)
                                                                                                                                 aGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAg                <  2:2641565/150‑1 (MQ=255)
                                                                                                                                 aGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAg                <  2:3889766/150‑1 (MQ=255)
                                                                                                                                 aGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAg                <  2:1840507/150‑1 (MQ=255)
                                                                                                                                  gCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGc               >  1:692898/1‑150 (MQ=255)
                                                                                                                                   cccccTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCa              >  2:629299/1‑150 (MQ=255)
                                                                                                                                    ccccTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCaa             <  2:410580/150‑1 (MQ=255)
                                                                                                                                    ccccTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCaa             <  1:618769/150‑1 (MQ=255)
                                                                                                                                     cccTGTTGAATATACTTGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAAc            <  2:4527887/150‑1 (MQ=255)
                                                                                                                                     cccTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAAc            <  2:4527441/150‑1 (MQ=255)
                                                                                                                                       cTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCt          <  1:2167825/150‑1 (MQ=255)
                                                                                                                                        tGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTg         <  1:1385947/150‑1 (MQ=255)
                                                                                                                                         gTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt        >  2:4400835/1‑150 (MQ=255)
                                                                                                                                          ttGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGtt       >  1:923668/1‑150 (MQ=255)
                                                                                                                                          ttGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGtt       >  1:976582/1‑150 (MQ=255)
                                                                                                                                          ttGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGtt       <  1:287381/150‑1 (MQ=255)
                                                                                                                                          ttGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGtt       <  2:2067451/150‑1 (MQ=255)
                                                                                                                                           tGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTg      >  1:1701980/1‑150 (MQ=255)
                                                                                                                                            gAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGc     >  2:4413467/1‑150 (MQ=255)
                                                                                                                                            gAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGc     >  2:2173627/1‑150 (MQ=255)
                                                                                                                                             aaTATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGct    >  2:698038/1‑150 (MQ=255)
                                                                                                                                             aaTATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGct    <  1:3515883/150‑1 (MQ=255)
                                                                                                                                               tatgCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGctct  >  2:1287425/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCT  >  CP009273/1116818‑1117110

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: