Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,924,898 T→A 8.6% D36V (GAT→GTT)  yebG ← DNA damage‑inducible protein regulated by LexA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,924,8980TA 8.6% 1044.6 / 63.8 309D36V (GAT→GTT) yebGDNA damage‑inducible protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (120/161);  new base A (19/7);  total (140/169)
Fisher's exact test for biased strand distribution p-value = 3.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 316 total aligned reads displayed.

CAGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATATTTC  >  CP009273/1924749‑1925047
                                                                                                                                                     |                                                                                                                                                     
cAGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCATCAAGGTGTCAAGCAGAt                                                                                                                                                       >  2:5057798/1‑150 (MQ=255)
cAGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAt                                                                                                                                                       <  1:4057832/150‑1 (MQ=255)
 aGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATc                                                                                                                                                      <  2:2530508/150‑1 (MQ=255)
  gCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATcc                                                                                                                                                     >  2:1581997/1‑150 (MQ=255)
   cACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAAccg                                                                                                                                                    >  1:2547226/1‑150 (MQ=255)
     ccGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccgcc                                                                                                                                                  <  1:5477943/150‑1 (MQ=255)
     ccGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccgcc                                                                                                                                                  <  1:74641/150‑1 (MQ=255)
     ccGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccgcc                                                                                                                                                  <  1:83472/150‑1 (MQ=255)
      cgacgaCCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccgccg                                                                                                                                                 >  1:1201114/1‑150 (MQ=255)
       gacgacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGt                                                                                                                                                >  1:959101/1‑150 (MQ=255)
         cgacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTAt                                                                                                                                              >  1:3582787/1‑150 (MQ=255)
          gagcTGTGGGGTCGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATc                                                                                                                                             <  2:3503157/147‑1 (MQ=255)
           acCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCa                                                                                                                                            <  2:3817316/150‑1 (MQ=255)
              tgtgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGc                                                                                                                                         <  1:2173221/150‑1 (MQ=255)
              tgtgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGc                                                                                                                                         >  1:1829111/1‑150 (MQ=255)
               gtgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCa                                                                                                                                        >  2:5477280/1‑150 (MQ=255)
                 ggggACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCAtt                                                                                                                                      >  2:3275032/1‑150 (MQ=255)
                  gggACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCAttt                                                                                                                                     <  1:1284770/150‑1 (MQ=255)
                   ggACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCAtttt                                                                                                                                    <  1:2500999/150‑1 (MQ=255)
                         tAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGTGCTTCACGTTGCTCGTCTTCAATTTGAACTGGAGAATTTGTCAGACAGGTGTCAAGAAGATCCGCCGAAACAAGCATTTTGTAATa                                                                                                                              >  1:2944004/1‑150 (MQ=255)
                         tAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATa                                                                                                                              <  2:1266528/150‑1 (MQ=255)
                            tttGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGc                                                                                                                           >  1:3827426/1‑150 (MQ=255)
                             ttGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCa                                                                                                                          >  2:4519959/1‑150 (MQ=255)
                                 cGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCgg                                                                                                                      >  2:2065901/1‑150 (MQ=255)
                                  ggTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGc                                                                                                                     >  2:2597327/1‑150 (MQ=255)
                                   gTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCt                                                                                                                    <  1:46739/150‑1 (MQ=255)
                                        cAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCtt                                                                                                               <  2:4759123/150‑1 (MQ=255)
                                         aGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttt                                                                                                              >  1:1044631/1‑150 (MQ=255)
                                         aGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttt                                                                                                              <  2:4335991/150‑1 (MQ=255)
                                           aaTGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttttt                                                                                                            <  1:4139568/150‑1 (MQ=255)
                                           aaTGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttttt                                                                                                            <  2:425681/150‑1 (MQ=255)
                                           aaTGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttttt                                                                                                            <  1:55074/150‑1 (MQ=255)
                                             tGGTTCTCAGCACATCATTTTTTTCTGCAAGCCATATCGAAATGGCTTCACGTTGCTAGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGc                                                                                                          <  2:2942472/150‑1 (MQ=255)
                                             tGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGc                                                                                                          <  1:1655618/150‑1 (MQ=255)
                                                 gCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGt                                                                                                      <  2:1821644/150‑1 (MQ=255)
                                                 gCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGt                                                                                                      >  2:3596984/1‑150 (MQ=255)
                                                     aGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAAc                                                                                                  >  2:5118645/1‑150 (MQ=255)
                                                     aGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGACCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAAc                                                                                                  >  1:2723977/1‑150 (MQ=255)
                                                      gCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACg                                                                                                 <  1:2707674/150‑1 (MQ=255)
                                                      gCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACg                                                                                                 >  2:5401603/1‑150 (MQ=255)
                                                             tttttGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAtttt                                                                                          <  1:568704/150‑1 (MQ=255)
                                                                 tGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTctct                                                                                      >  2:4889940/1‑150 (MQ=255)
                                                                   ttCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTc                                                                                    <  2:1139635/150‑1 (MQ=255)
                                                                        ccagccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCt                                                                               <  2:1056289/150‑1 (MQ=255)
                                                                         cagccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTc                                                                              <  1:147613/150‑1 (MQ=255)
                                                                         cagccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTc                                                                              <  1:2131116/150‑1 (MQ=255)
                                                                          agccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                             >  1:2821284/1‑150 (MQ=255)
                                                                          agccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                             <  1:2089796/150‑1 (MQ=255)
                                                                           gccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCGTCGCCCTCAc                                                                            <  2:2322093/150‑1 (MQ=255)
                                                                                  cGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGAc                                                                     <  2:2872345/150‑1 (MQ=255)
                                                                                   gAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACt                                                                    <  2:2676191/150‑1 (MQ=255)
                                                                                   gAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACt                                                                    >  2:2726814/1‑150 (MQ=255)
                                                                                   gAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACt                                                                    <  1:485751/150‑1 (MQ=255)
                                                                                    aaaGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTa                                                                   >  2:2945065/1‑150 (MQ=255)
                                                                                      aGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACg                                                                 <  2:1404/150‑1 (MQ=255)
                                                                                       gggCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGt                                                                <  1:5404198/150‑1 (MQ=255)
                                                                                           ttCACGTTGCTCGTCGTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTAttt                                                            <  2:5501863/150‑1 (MQ=255)
                                                                                             cACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTaa                                                          <  2:2145111/150‑1 (MQ=255)
                                                                                               cGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACt                                                        >  1:3260690/1‑150 (MQ=255)
                                                                                                gTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACtt                                                       <  2:4948000/150‑1 (MQ=255)
                                                                                                 ttGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTc                                                      >  2:1609191/1‑150 (MQ=255)
                                                                                                  tGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCa                                                     <  2:3023330/150‑1 (MQ=255)
                                                                                                      cGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCg                                                 >  1:867938/1‑150 (MQ=255)
                                                                                                      cGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCg                                                 >  2:5330796/1‑150 (MQ=255)
                                                                                                        tCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACAGcc                                               >  2:3745569/1‑150 (MQ=255)
                                                                                                        tCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGcc                                               >  2:2448710/1‑150 (MQ=255)
                                                                                                        tCTTCCATTTGCACTGGAGAATATGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTACCTGTTTGCTGGTAAAAGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAAACGcc                                               >  2:901976/1‑150 (MQ=255)
                                                                                                            ccATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCAttt                                           >  1:3711285/1‑150 (MQ=255)
                                                                                                            ccATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCAttt                                           >  2:755530/1‑150 (MQ=255)
                                                                                                             cATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTg                                          <  1:305188/150‑1 (MQ=255)
                                                                                                              aTTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGc                                         >  1:5294265/1‑150 (MQ=255)
                                                                                                              aTTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGc                                         >  1:303645/1‑150 (MQ=255)
                                                                                                                  gCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGcc                                     >  2:1952808/1‑150 (MQ=255)
                                                                                                                  gCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGcc                                     >  1:2966267/1‑150 (MQ=255)
                                                                                                                   cACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCct                                    >  1:853664/1‑150 (MQ=255)
                                                                                                                    aCTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCctc                                   <  2:2111789/150‑1 (MQ=255)
                                                                                                                       ggAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCa                                >  1:5549848/1‑150 (MQ=255)
                                                                                                                        gagaATTTGTCAGCATTGTGTATTGCTGATACGACTTATAAAGTTTTTTGTCATAAGAATCGGATTTCTTTTTGATGGTAAACGTCATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCat                               <  1:4387365/150‑1 (MQ=255)
                                                                                                                          gaATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataa                             >  2:3696767/1‑150 (MQ=255)
                                                                                                                          gaATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATAAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataa                             <  2:3970706/150‑1 (MQ=255)
                                                                                                                           aaTTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataat                            <  1:2726814/150‑1 (MQ=255)
                                                                                                                           aaTTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataat                            <  2:5197913/150‑1 (MQ=255)
                                                                                                                             tttGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCAtaataa                          <  1:1326619/150‑1 (MQ=255)
                                                                                                                               tGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACt                        >  1:4599736/1‑150 (MQ=255)
                                                                                                                                    gCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAt                   <  2:5376778/150‑1 (MQ=255)
                                                                                                                                     ccAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATACGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAtt                  <  2:799708/150‑1 (MQ=255)
                                                                                                                                     ccAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAtt                  <  1:1701179/150‑1 (MQ=255)
                                                                                                                                      cAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAttt                 >  2:743477/1‑150 (MQ=255)
                                                                                                                                       aGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAtttt                <  1:4987313/150‑1 (MQ=255)
                                                                                                                                        ggTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATGTta               >  2:3374594/1‑150 (MQ=255)
                                                                                                                                         gtgtCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTtat              <  1:994734/150‑1 (MQ=255)
                                                                                                                                           gtCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATAc            <  2:2799452/150‑1 (MQ=255)
                                                                                                                                               aGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGta        >  1:3397998/1‑150 (MQ=255)
                                                                                                                                               aGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGta        >  1:5201204/1‑150 (MQ=255)
                                                                                                                                                gCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtat       <  1:5091652/150‑1 (MQ=255)
                                                                                                                                                 cAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtata      <  2:4607878/150‑1 (MQ=255)
                                                                                                                                                 cAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtata      <  1:1255592/150‑1 (MQ=255)
                                                                                                                                                 cAGATACGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtata      <  1:5047597/150‑1 (MQ=255)
                                                                                                                                                  aGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTatat     <  1:1731096/150‑1 (MQ=255)
                                                                                                                                                   gATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATAtt    <  1:2412158/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATATTTC  >  CP009273/1924749‑1925047

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: