Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 1,407,139 A→T C118* (TGT→TGA intR ← Rac prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,407,1390AT99.6% 835.0 / ‑6.5 250C118* (TGT→TGAintRRac prophage; integrase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (121/128);  minor base C (1/0);  total (122/128)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

CAGGGCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGT  >  CP009273/1406992‑1407285
                                                                                                                                                   |                                                                                                                                                  
cAGGGCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTctca                                                                                                                                                  >  2:2600609/1‑150 (MQ=255)
  gggCGATATAGTTATCAACTGTGCGCACGGTTCGTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCaaa                                                                                                                                                <  1:3076933/150‑1 (MQ=255)
   ggCGATATAGTTATCTACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAat                                                                                                                                               <  1:1232190/150‑1 (MQ=255)
   ggCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAat                                                                                                                                               <  2:2192031/150‑1 (MQ=255)
    gCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAatt                                                                                                                                              <  1:1971278/150‑1 (MQ=255)
     cGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAatta                                                                                                                                             >  1:1889324/1‑150 (MQ=255)
         ataGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAattattat                                                                                                                                         >  2:3797205/1‑150 (MQ=255)
         ataGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAattattat                                                                                                                                         >  2:152872/1‑150 (MQ=255)
           aGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTc                                                                                                                                       <  2:1170693/150‑1 (MQ=255)
            gTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCt                                                                                                                                      <  2:2158527/150‑1 (MQ=255)
              tATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGa                                                                                                                                    <  2:3028355/150‑1 (MQ=255)
                tCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAgt                                                                                                                                  <  1:3776176/150‑1 (MQ=255)
                tCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAgt                                                                                                                                  <  2:2001231/150‑1 (MQ=255)
                 cAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAgtg                                                                                                                                 >  2:3472946/1‑150 (MQ=255)
                 cAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAgtg                                                                                                                                 <  2:3350609/150‑1 (MQ=255)
                   aCTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGtt                                                                                                                               <  2:1894172/150‑1 (MQ=255)
                     tgtgCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGtttt                                                                                                                             >  1:473269/1‑150 (MQ=255)
                       tgCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTa                                                                                                                           >  2:556248/1‑150 (MQ=255)
                              ggTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGa                                                                                                                    >  2:2839375/1‑150 (MQ=255)
                                ttcttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGAtt                                                                                                                  <  1:1327936/150‑1 (MQ=255)
                                ttcttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGAtt                                                                                                                  <  1:3228112/150‑1 (MQ=255)
                                ttcttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGAtt                                                                                                                  <  2:2624738/150‑1 (MQ=255)
                                ttcttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGAtt                                                                                                                  <  1:1129567/150‑1 (MQ=255)
                                  cttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTc                                                                                                                <  1:1097317/150‑1 (MQ=255)
                                   ttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCg                                                                                                               <  2:2811930/150‑1 (MQ=255)
                                   ttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCg                                                                                                               <  2:1633594/150‑1 (MQ=255)
                                     ccTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGtt                                                                                                             <  2:3025015/150‑1 (MQ=255)
                                      cTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGttt                                                                                                            >  2:937028/1‑150 (MQ=255)
                                         ttttGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttct                                                                                                         <  2:2763246/150‑1 (MQ=255)
                                         ttttGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttct                                                                                                         <  2:3311505/150‑1 (MQ=255)
                                         ttttGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttct                                                                                                         <  1:3751593/150‑1 (MQ=255)
                                           ttgttgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttcttc                                                                                                       <  2:561209/150‑1 (MQ=255)
                                            tgttgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTATTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTATGAGTGTTTTTAATTGTGATTTCGTTTTATTCa                                                                                                      >  2:2668488/1‑150 (MQ=255)
                                            tgttgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCa                                                                                                      <  2:2318094/150‑1 (MQ=255)
                                             gttgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCAt                                                                                                     <  1:3698853/150‑1 (MQ=255)
                                              ttgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCAtt                                                                                                    <  2:3654457/150‑1 (MQ=255)
                                               tgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATtg                                                                                                   <  1:2396435/150‑1 (MQ=255)
                                                ggATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATtgt                                                                                                  >  2:137921/1‑150 (MQ=255)
                                                 gATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATtgtg                                                                                                 <  1:848267/150‑1 (MQ=255)
                                                    cTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGtgttt                                                                                              >  1:1543904/1‑150 (MQ=255)
                                                    cTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGtgttt                                                                                              >  1:3249186/1‑150 (MQ=255)
                                                       ggATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGtgtttgtt                                                                                           <  1:2486190/150‑1 (MQ=255)
                                                        gATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTttgttg                                                                                          >  2:1295660/1‑150 (MQ=255)
                                                         aTTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTtgttgt                                                                                         <  2:873305/150‑1 (MQ=255)
                                                         aTTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTtgttgt                                                                                         <  2:1006646/150‑1 (MQ=255)
                                                         aTTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTtgttgt                                                                                         <  2:384794/150‑1 (MQ=255)
                                                            aTCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTaa                                                                                      >  2:1612364/1‑150 (MQ=255)
                                                             tCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAAc                                                                                     >  1:73570/1‑150 (MQ=255)
                                                                aGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAActct                                                                                  >  1:1298120/1‑150 (MQ=255)
                                                                aGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAActct                                                                                  <  2:1807307/150‑1 (MQ=255)
                                                                aGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAActct                                                                                  <  1:2326825/150‑1 (MQ=255)
                                                                aGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAActct                                                                                  <  1:3599557/150‑1 (MQ=255)
                                                                aGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAActct                                                                                  >  2:405183/1‑150 (MQ=255)
                                                                 tgTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTg                                                                                 >  2:1159906/2‑150 (MQ=255)
                                                                      aGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTTAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTctct                                                                            >  2:2049713/1‑150 (MQ=255)
                                                                      aGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTctct                                                                            >  2:738279/1‑150 (MQ=255)
                                                                       gCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCtt                                                                           <  2:1433415/150‑1 (MQ=255)
                                                                        cGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCttt                                                                          >  1:2080619/1‑150 (MQ=255)
                                                                        cGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCttt                                                                          <  2:3469759/150‑1 (MQ=255)
                                                                         gTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTa                                                                         >  2:2988960/1‑150 (MQ=255)
                                                                           ttCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAAt                                                                       <  2:3565501/150‑1 (MQ=255)
                                                                           ttCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAAt                                                                       <  1:1617583/150‑1 (MQ=255)
                                                                           ttCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAAt                                                                       >  2:309564/1‑150 (MQ=255)
                                                                              tcCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAAtttt                                                                    <  1:1169041/150‑1 (MQ=255)
                                                                               ccATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAAttttt                                                                   <  2:1678393/150‑1 (MQ=255)
                                                                                 aTGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGt                                                                 <  2:659265/150‑1 (MQ=255)
                                                                                  tGCAGCCGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTc                                                                >  2:2560180/1‑150 (MQ=255)
                                                                                  tGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTc                                                                <  1:3797205/150‑1 (MQ=255)
                                                                                    cagcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCa                                                              <  1:2587885/150‑1 (MQ=255)
                                                                                     agcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAg                                                             >  2:1226906/1‑150 (MQ=255)
                                                                                     agcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAg                                                             >  1:1705912/1‑150 (MQ=255)
                                                                                     agcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAg                                                             >  1:2239817/1‑150 (MQ=255)
                                                                                     agcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAg                                                             >  1:2771369/1‑150 (MQ=255)
                                                                                     agcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAg                                                             >  1:2847862/1‑150 (MQ=255)
                                                                                     agcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAg                                                             <  1:933330/150‑1 (MQ=255)
                                                                                        agTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGata                                                          >  2:2304134/1‑150 (MQ=255)
                                                                                         gTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGatat                                                         <  2:432440/150‑1 (MQ=255)
                                                                                                    tAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGctct                                              >  2:788405/1‑150 (MQ=255)
                                                                                                     aGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCtctc                                             >  1:1478010/1‑150 (MQ=255)
                                                                                                     aGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCtctc                                             <  1:2988960/150‑1 (MQ=255)
                                                                                                          aaGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAAc                                        >  2:2475299/1‑150 (MQ=255)
                                                                                                              tatCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGttt                                    <  1:2652264/150‑1 (MQ=255)
                                                                                                                  gCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTAt                                <  1:2368683/150‑1 (MQ=255)
                                                                                                                    tataACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGa                              >  1:702497/1‑150 (MQ=255)
                                                                                                                    tataACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGa                              >  1:2786542/1‑150 (MQ=255)
                                                                                                                         cGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACtctc                         >  1:3161905/1‑150 (MQ=255)
                                                                                                                          gAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCg                        >  1:57865/1‑150 (MQ=255)
                                                                                                                          gAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCg                        >  2:3707237/1‑150 (MQ=255)
                                                                                                                            atatGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGtt                      <  1:1039829/150‑1 (MQ=255)
                                                                                                                               tgtgATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTc                   >  2:765044/1‑150 (MQ=255)
                                                                                                                                 tgATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCAc                 >  1:3363218/1‑150 (MQ=255)
                                                                                                                                  gATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCAcc                >  1:143876/1‑150 (MQ=255)
                                                                                                                                     atCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAtt             <  2:3113305/150‑1 (MQ=255)
                                                                                                                                        ggggTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAttttt          <  1:101786/150‑1 (MQ=255)
                                                                                                                                        ggggTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAttttt          <  2:1057553/150‑1 (MQ=255)
                                                                                                                                         gggTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTg         >  2:1271419/1‑150 (MQ=255)
                                                                                                                                           gTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGcc       >  2:708549/1‑150 (MQ=255)
                                                                                                                                             aCTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGcccc     >  1:2719591/1‑150 (MQ=255)
                                                                                                                                              cTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCa    <  1:196877/150‑1 (MQ=255)
                                                                                                                                               tttCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAg   >  1:3529437/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGGGCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGT  >  CP009273/1406992‑1407285

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: