Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,317,200 G→T 11.1% Q2K (CAA→AAA)  trpE ← component I of anthranilate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,317,2000GT11.1% 614.0 / 2.3 202Q2K (CAA→AAA) trpEcomponent I of anthranilate synthase
Reads supporting (aligned to +/- strand):  ref base G (49/121);  new base T (30/2);  total (79/123)
Fisher's exact test for biased strand distribution p-value = 2.55e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

GCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCATTGTCGATACC  >  CP009273/1317051‑1317349
                                                                                                                                                     |                                                                                                                                                      
gCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTg                                                                                                                                                        >  1:3503184/1‑150 (MQ=255)
  tttttAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCa                                                                                                                                                      >  1:2357710/1‑150 (MQ=255)
  tttttAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCa                                                                                                                                                      >  1:1878902/1‑150 (MQ=255)
   ttttAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTTCAt                                                                                                                                                     >  1:2153894/1‑150 (MQ=255)
     ttAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTg                                                                                                                                                   <  2:1606716/150‑1 (MQ=255)
     ttAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTg                                                                                                                                                   <  1:468962/150‑1 (MQ=255)
       aaaTCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGtt                                                                                                                                                 <  2:2351315/150‑1 (MQ=255)
       aaaTCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGtt                                                                                                                                                 <  2:4729741/150‑1 (MQ=255)
       aaaTCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGtt                                                                                                                                                 <  1:604559/150‑1 (MQ=255)
        aaTCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGACTTGTTa                                                                                                                                                >  2:1468125/1‑150 (MQ=255)
         atcatcTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTAt                                                                                                                                               <  1:3254287/150‑1 (MQ=255)
          tcatcTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTAtt                                                                                                                                              <  1:428187/150‑1 (MQ=255)
          tcatcTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTAtt                                                                                                                                              <  1:4381541/150‑1 (MQ=255)
           catcTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATtc                                                                                                                                             <  1:3604147/150‑1 (MQ=255)
            atcTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATtct                                                                                                                                            >  2:766225/1‑150 (MQ=255)
            atcTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATtct                                                                                                                                            <  2:3515264/150‑1 (MQ=255)
              cTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTctct                                                                                                                                          <  2:4531019/150‑1 (MQ=255)
              cTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGGTGTTTTGTTTTGATTTTTTTTgtct                                                                                                                                          >  1:461765/1‑150 (MQ=255)
               tttGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTTCATTGTTTTTCTCTa                                                                                                                                         >  2:2295270/1‑150 (MQ=255)
               tttGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTTCATTGTTATTCTCTa                                                                                                                                         >  2:4043621/1‑150 (MQ=255)
                 tGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAAt                                                                                                                                       <  1:868659/150‑1 (MQ=255)
                 tGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAAt                                                                                                                                       <  1:1924448/150‑1 (MQ=255)
                 tGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGGTTTTGTGGTTTCATTGTTATTCTCTAAt                                                                                                                                       >  2:1299004/1‑150 (MQ=255)
                  gCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAAtt                                                                                                                                      <  2:436427/150‑1 (MQ=255)
                   cTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTTTGTTTGGATTGTTATTCTCTAAttt                                                                                                                                     >  2:3307222/1‑150 (MQ=255)
                       cGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGtt                                                                                                                                 <  1:917317/150‑1 (MQ=255)
                       cGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGtt                                                                                                                                 <  1:2180334/150‑1 (MQ=255)
                       cGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGtt                                                                                                                                 <  2:2595129/150‑1 (MQ=255)
                           atCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCaaa                                                                                                                            <  1:2976638/150‑1 (MQ=255)
                              tGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCaaaaaa                                                                                                                         <  1:2490523/150‑1 (MQ=255)
                                cGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGGGTTTTTGTTTTGTTTGTTATTGTATAAATTTTTTAaaaaaaaa                                                                                                                       >  1:806442/1‑150 (MQ=255)
                                     tCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCg                                                                                                                  <  1:4094708/150‑1 (MQ=255)
                                      ccAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTTCAGTTTTATTCTCTAATTTTTTTAAAAAAAAATCCCGc                                                                                                                 >  1:314902/1‑150 (MQ=255)
                                      ccAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAAGCCCg                                                                                                                  <  1:4300047/150‑1 (MQ=255)
                                       cagcagcagCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGGATGGTTTTTTTGGAATTTTGTTCAAAAAAAATCCCGCt                                                                                                                >  2:3187624/1‑150 (MQ=255)
                                        agcagcagcGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTc                                                                                                               <  1:3969025/150‑1 (MQ=255)
                                          cagcagcGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTTGATGGTTATTCTCTAATTTTGTTCAAAAAAAATCCCCCGCAt                                                                                                             >  1:2665667/1‑150 (MQ=255)
                                          cagcagcGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCAt                                                                                                             <  2:346498/150‑1 (MQ=255)
                                           agcagcGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCCTTGTTATTCTCTAATTTTTTTCAAAAAAAATCCCGCTCAtt                                                                                                            >  1:3488652/1‑150 (MQ=255)
                                           agcagcGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGGTTTTGTGTTTGGTTTTTTTTTCTCGTATTTTGTTGACAAAATATCCCGCTCagt                                                                                                            >  2:3602369/1‑148 (MQ=255)
                                             cagcGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTTTGTTTTGATTGTTATTAGCTAATTTTGTTAAAAAAAAAACCCCCGCAATAg                                                                                                          >  1:2454260/1‑150 (MQ=255)
                                                 gTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCgg                                                                                                      <  1:1984230/150‑1 (MQ=255)
                                                  gtGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCggg                                                                                                     <  1:3402612/149‑1 (MQ=255)
                                                   tGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGc                                                                                                    <  1:3818129/150‑1 (MQ=255)
                                                    gCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCt                                                                                                   <  2:3965266/150‑1 (MQ=255)
                                                    gCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCt                                                                                                   <  2:1577033/150‑1 (MQ=255)
                                                     ccGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTTCATTGTTATTCTCTAATTTTTTTAAAAAAAAATCCCGCTCATTAGGCGGGCTg                                                                                                  >  2:2164284/1‑150 (MQ=255)
                                                      cGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTGTTGTGTTTTGATTGTTATTCTCTAAATTTGTTCAAAAAAAATCCCGCTCATTAGGCGGGCTgg                                                                                                 >  2:430151/1‑150 (MQ=255)
                                                        gACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGt                                                                                               <  2:2153894/150‑1 (MQ=255)
                                                        gACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGt                                                                                               <  1:2621215/150‑1 (MQ=255)
                                                          cGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTTTGTTTTCATTTTTATTATATAATTTTGTTGAAAAAAAATACCGCTCATTAGGCGGGCTTGGGAt                                                                                             >  2:2506805/1‑150 (MQ=255)
                                                             tCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTg                                                                                          <  1:1816521/150‑1 (MQ=255)
                                                             tCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTg                                                                                          <  2:11901/150‑1 (MQ=255)
                                                               cccACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGGTATTATCTAAATTTTTTCAAAAAAAAACCCCCTCCATAGGCGGGCTGGGGATCTGAt                                                                                        >  2:1250352/1‑150 (MQ=255)
                                                                ccACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGGATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCAATAGGCGGGCGGGGTATATGAtt                                                                                       >  2:150248/1‑150 (MQ=255)
                                                                 cacacaACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTg                                                                                      <  1:2764334/150‑1 (MQ=255)
                                                                  acacaACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAATCCCGCTCATTAGGCGGGCTGGGTATCTGATTGc                                                                                     >  2:4521738/1‑150 (MQ=255)
                                                                  acacaACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTATCTAATTTTGTTAAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATATGATTGc                                                                                     >  2:4521455/1‑150 (MQ=255)
                                                                    acaACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTAAAAAAAAATCCCCCTCATTAGGCGGGCTGGGTATCTGATTGCtt                                                                                   >  2:1537205/1‑150 (MQ=255)
                                                                     caACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTTTTTTTTGATTGTTTTTCTCTAATTTTTTTAAAAAAAAATCCCGCTCATAAGGCGGGCGGGGTTTCTGTATGCttt                                                                                  >  1:3796074/1‑150 (MQ=255)
                                                                     caACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTa                                                                                 <  2:2437399/150‑1 (MQ=255)
                                                                      aaCTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTa                                                                                 <  2:2739194/150‑1 (MQ=255)
                                                                                  aGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAAt                                                                     <  1:1609780/150‑1 (MQ=255)
                                                                                    cgcgGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATac                                                                   <  1:928315/150‑1 (MQ=255)
                                                                                      cgGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATacac                                                                 <  1:2614998/150‑1 (MQ=255)
                                                                                            gATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCg                                                           <  1:4447535/150‑1 (MQ=255)
                                                                                              ttGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGtt                                                         <  1:2627937/150‑1 (MQ=255)
                                                                                               tGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGttt                                                        <  2:3925098/150‑1 (MQ=255)
                                                                                                gTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTTCATTGTTATTCTCTAATTTTGTTAAAAAAAAATCCCGCGCATTATGCGGGCTTGGTATATGATTGCTTTAAGCATGGTGAATACACTGCCCGgtta                                                       >  2:1854235/1‑149 (MQ=255)
                                                                                                  cgcgATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAg                                                     <  1:4669503/150‑1 (MQ=255)
                                                                                                   gcgATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAgg                                                    <  1:3180449/150‑1 (MQ=255)
                                                                                                   gcgATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAgg                                                    <  1:3883537/150‑1 (MQ=255)
                                                                                                    cgATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGa                                                   <  1:459101/150‑1 (MQ=255)
                                                                                                      aTAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAATCCCCCGCATTATGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGTATACACAGCCCGTTTAAGGAAg                                                 >  1:2306891/1‑150 (MQ=255)
                                                                                                      aTAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAg                                                 <  2:817157/150‑1 (MQ=255)
                                                                                                        aaGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTg                                               <  2:842052/150‑1 (MQ=255)
                                                                                                        aaGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTg                                               <  1:1150853/150‑1 (MQ=255)
                                                                                                        aaGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTg                                               <  2:3213055/150‑1 (MQ=255)
                                                                                                         aGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTATAATTTTGTTCAAAAAAAATCCCGCTCATTATGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATAAACTGCCCGTTTCAGGAAGTgc                                              >  2:3901050/1‑150 (MQ=255)
                                                                                                                tCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCGTTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGccacca                                       <  2:2319696/150‑1 (MQ=255)
                                                                                                                tCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGccacca                                       <  2:2319596/150‑1 (MQ=255)
                                                                                                                tCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGccacca                                       <  2:2051820/150‑1 (MQ=255)
                                                                                                                   cAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAAcc                                    <  2:3373359/150‑1 (MQ=255)
                                                                                                                   cAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAAcc                                    <  2:2783230/150‑1 (MQ=255)
                                                                                                                        tAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCa                               <  2:2430633/150‑1 (MQ=255)
                                                                                                                           cAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTa                            <  2:911355/150‑1 (MQ=255)
                                                                                                                            aGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTAc                           <  2:963813/150‑1 (MQ=255)
                                                                                                                                 gagaGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAATAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGaaaa                      <  2:1643538/150‑1 (MQ=255)
                                                                                                                                 gagaGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGaaaa                      <  1:4501066/150‑1 (MQ=255)
                                                                                                                                    agTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTg                   <  2:3790154/150‑1 (MQ=255)
                                                                                                                                     gTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGc                  >  2:2209961/1‑150 (MQ=255)
                                                                                                                                     cacGGTTTTTGTGATTGCAATCATATTCTCTAATAATCATCAAAAAAAAGCCCGCTCAATACGCGGGCTGGGTATCAGAAAGCTTAACGAATGGTGAATACACTGCCCGTAACAGGAAGTGCGCCACCAACCTTACAGAACGAAAATAGc                  >  2:4534249/3‑150 (MQ=255)
                                                                                                                                        ggTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCttt               <  1:4066940/150‑1 (MQ=255)
                                                                                                                                           ttttGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCAt            <  2:3157815/150‑1 (MQ=255)
                                                                                                                                            tttGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCAtt           >  1:4117790/1‑150 (MQ=255)
                                                                                                                                                  tttGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCATTGTCGAt     <  1:2873639/150‑1 (MQ=255)
                                                                                                                                                   ttGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCATTGTCGATa    >  1:1007730/1‑150 (MQ=255)
                                                                                                                                                    tGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCATTGTCGATAc   <  1:2263898/150‑1 (MQ=255)
                                                                                                                                                    tGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCATTGTCGATAc   >  1:1474673/1‑150 (MQ=255)
                                                                                                                                                     gCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTACCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCATTGTCGATAcc  >  1:4211734/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCTGGGTATCTGATTGCTTTACGCATGGTGAATACACTGCCCGTTTCAGGAAGTGCGCCACCAACCTTTCAGTACGAAAATTGCTTTCATTGTCGATACC  >  CP009273/1317051‑1317349

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: