Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,311,385 T→A 7.3% Q32L (CAG→CTG)  trpA ← tryptophan synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,311,3850TA 7.3% 988.0 / 49.8 282Q32L (CAG→CTG) trpAtryptophan synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (119/142);  new base A (3/18);  total (122/160)
Fisher's exact test for biased strand distribution p-value = 5.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 288 total aligned reads displayed.

AACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGC  >  CP009273/1311236‑1311534
                                                                                                                                                     |                                                                                                                                                     
aaCATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACt                                                                                                                                                       <  1:730669/150‑1 (MQ=255)
aaCATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACt                                                                                                                                                       >  1:3035979/1‑150 (MQ=255)
  cATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGc                                                                                                                                                     >  2:3760784/1‑150 (MQ=255)
   aTTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCt                                                                                                                                                    <  1:4395096/150‑1 (MQ=255)
   aTTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCt                                                                                                                                                    >  1:2969905/1‑150 (MQ=255)
     tgtgCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCa                                                                                                                                                  <  1:2869095/150‑1 (MQ=255)
     tgtgCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCa                                                                                                                                                  >  2:2106032/1‑150 (MQ=255)
         ccGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGc                                                                                                                                              <  1:993637/150‑1 (MQ=255)
           ggAGTCACACCTGCCGCAAAGGCGCGCAGTGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGccc                                                                                                                                            <  1:3114467/150‑1 (MQ=255)
              gTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGa                                                                                                                                         >  2:4023670/1‑150 (MQ=255)
                 acacCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCa                                                                                                                                      <  2:2594533/150‑1 (MQ=255)
                  cacCTGCCGCATAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGTGTCGGAGTAGTGGATACCTAACTCCAGCGCGTCAGCACCGGATTCAATTAGCGTATCGATATTTTTCAATGACTGCTCAATGCCCGGATCAc                                                                                                                                     <  2:3049307/150‑1 (MQ=255)
                   acCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCATTGACTGCTCAATGCCCGGATCAcc                                                                                                                                    >  2:4521075/1‑150 (MQ=255)
                    ccTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCg                                                                                                                                   >  2:1798656/1‑150 (MQ=255)
                     cTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGa                                                                                                                                  >  1:902649/1‑150 (MQ=255)
                     cTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGa                                                                                                                                  <  2:3529962/150‑1 (MQ=255)
                        ccGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCg                                                                                                                               <  2:743993/150‑1 (MQ=255)
                          gCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTg                                                                                                                             >  2:4755852/1‑150 (MQ=255)
                          gCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTg                                                                                                                             <  1:4173044/150‑1 (MQ=255)
                             aaGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACg                                                                                                                          <  1:5220852/150‑1 (MQ=255)
                              aGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGa                                                                                                                         >  2:161393/1‑150 (MQ=255)
                                  gcgcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAgg                                                                                                                     <  2:894684/150‑1 (MQ=255)
                                    gcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGaa                                                                                                                   <  1:194605/150‑1 (MQ=255)
                                      agagTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACg                                                                                                                 >  1:3744912/1‑150 (MQ=255)
                                       gagTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGa                                                                                                                >  1:5198494/1‑150 (MQ=255)
                                         gTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAAt                                                                                                              >  2:3906678/1‑150 (MQ=255)
                                          tGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATg                                                                                                             >  2:2998262/1‑150 (MQ=255)
                                           ggCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATgc                                                                                                            >  2:440141/1‑150 (MQ=255)
                                                tttGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCctt                                                                                                       >  1:1653595/1‑150 (MQ=255)
                                                  tGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCcttct                                                                                                     <  2:2253687/150‑1 (MQ=255)
                                                   gAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCcttctt                                                                                                    <  2:5012770/150‑1 (MQ=255)
                                                     aTCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTg                                                                                                  <  1:1815570/150‑1 (MQ=255)
                                                      tcgtcgGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTgc                                                                                                 <  1:1494863/150‑1 (MQ=255)
                                                           gggCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTcc                                                                                            >  1:16918/1‑150 (MQ=255)
                                                            ggCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCt                                                                                           <  1:4890886/150‑1 (MQ=255)
                                                              ccATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTc                                                                                         >  1:3489725/1‑150 (MQ=255)
                                                                aTCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCaa                                                                                       <  2:4931433/150‑1 (MQ=255)
                                                                 tCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAAc                                                                                      <  1:2101729/150‑1 (MQ=255)
                                                                  ccGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACt                                                                                     <  1:4808942/150‑1 (MQ=255)
                                                                    gCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGATCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTgg                                                                                   >  2:518895/1‑150 (MQ=255)
                                                                      cAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGc                                                                                 >  1:4913967/1‑150 (MQ=255)
                                                                         tGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCaaa                                                                              >  2:4188034/1‑150 (MQ=255)
                                                                         tGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCaaa                                                                              >  2:1341174/1‑150 (MQ=255)
                                                                         tGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCaaa                                                                              >  1:3704476/1‑150 (MQ=255)
                                                                          gggTCTGAGTAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCTCCTTCTTTGCGCTCCTTCAACTGGGCAAAc                                                                             <  2:3636136/150‑1 (MQ=255)
                                                                           ggTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACa                                                                            <  2:1955231/150‑1 (MQ=255)
                                                                           ggTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACa                                                                            <  1:4587329/150‑1 (MQ=255)
                                                                            gTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACag                                                                           >  2:2367300/1‑150 (MQ=255)
                                                                               ggAGACGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAgaga                                                                        >  2:671136/1‑150 (MQ=255)
                                                                                gagaAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGAt                                                                       <  1:242228/150‑1 (MQ=255)
                                                                                  gaAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTc                                                                     <  1:2676541/150‑1 (MQ=255)
                                                                                  gaAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTc                                                                     >  1:333101/1‑150 (MQ=255)
                                                                                  gaAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTc                                                                     <  2:1663592/150‑1 (MQ=255)
                                                                                   aaGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCg                                                                    <  1:4273019/150‑1 (MQ=255)
                                                                                    aGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGt                                                                   <  1:2677300/150‑1 (MQ=255)
                                                                                     ggggATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTa                                                                  >  2:2273695/1‑150 (MQ=255)
                                                                                      gggATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAg                                                                 <  1:325605/150‑1 (MQ=255)
                                                                                      gggATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAg                                                                 <  2:2107389/150‑1 (MQ=255)
                                                                                       ggATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGc                                                                <  1:214666/150‑1 (MQ=255)
                                                                                         aTACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGGTTCGTAGCGt                                                              >  2:1157440/1‑150 (MQ=255)
                                                                                         aTACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGt                                                              <  2:3035642/150‑1 (MQ=255)
                                                                                         aTACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGt                                                              >  2:3167889/1‑150 (MQ=255)
                                                                                          tACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGtt                                                             <  2:4472980/150‑1 (MQ=255)
                                                                                               aaCTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATc                                                        >  1:3974061/1‑150 (MQ=255)
                                                                                               aaCTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATc                                                        >  1:1361799/1‑150 (MQ=255)
                                                                                                  tCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGa                                                     >  1:4780721/1‑150 (MQ=255)
                                                                                                    cAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGAtt                                                   <  2:509267/150‑1 (MQ=255)
                                                                                                     aGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGAttt                                                  <  2:3400977/150‑1 (MQ=255)
                                                                                                           tCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTc                                            <  1:3985782/150‑1 (MQ=255)
                                                                                                           tCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTc                                            <  2:4723923/150‑1 (MQ=255)
                                                                                                            cAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCg                                           <  1:1191634/150‑1 (MQ=255)
                                                                                                              gCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTg                                         <  2:2882014/150‑1 (MQ=255)
                                                                                                              gCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTg                                         <  1:2280931/150‑1 (MQ=255)
                                                                                                               cACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGc                                        <  2:3663663/150‑1 (MQ=255)
                                                                                                                   ggCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTc                                    >  1:1925626/1‑150 (MQ=255)
                                                                                                                    gCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCa                                   <  2:4119963/150‑1 (MQ=255)
                                                                                                                     cTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCaa                                  >  2:2528038/1‑150 (MQ=255)
                                                                                                                        cAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAat                               <  2:5154187/150‑1 (MQ=255)
                                                                                                                        cAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAat                               <  2:4356979/150‑1 (MQ=255)
                                                                                                                          aTTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAatat                             >  1:619712/1‑150 (MQ=255)
                                                                                                                            tAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCg                           >  2:4875275/1‑150 (MQ=255)
                                                                                                                             aGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGt                          >  1:2323071/1‑150 (MQ=255)
                                                                                                                                gTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGaa                       >  2:1524555/1‑150 (MQ=255)
                                                                                                                                gTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTGTCACAATCTCGGGaa                       >  2:4847949/1‑150 (MQ=255)
                                                                                                                                 tATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAAc                      <  2:3834840/150‑1 (MQ=255)
                                                                                                                                 tATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAAc                      <  1:4099541/150‑1 (MQ=255)
                                                                                                                                  aTCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACg                     <  2:3256135/150‑1 (MQ=255)
                                                                                                                                  aTCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACg                     >  1:313562/1‑150 (MQ=255)
                                                                                                                                    cGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGt                   >  2:1911878/1‑150 (MQ=255)
                                                                                                                                     gATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTg                  <  2:406139/150‑1 (MQ=255)
                                                                                                                                      aTAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGa                 >  1:390238/1‑150 (MQ=255)
                                                                                                                                       tAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGaa                <  2:3617053/150‑1 (MQ=255)
                                                                                                                                       tAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTCCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGaa                >  1:1248408/1‑150 (MQ=255)
                                                                                                                                           tttCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATg            >  2:3263316/1‑150 (MQ=255)
                                                                                                                                               aaTGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCtt        <  2:1453695/150‑1 (MQ=255)
                                                                                                                                                aTGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCttt       >  1:1097006/1‑150 (MQ=255)
                                                                                                                                                aTGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCttt       <  2:2010402/150‑1 (MQ=255)
                                                                                                                                                   aCTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATc    <  2:1105496/150‑1 (MQ=255)
                                                                                                                                                     tGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATcgc  <  1:3225320/150‑1 (MQ=255)
                                                                                                                                                     tGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATcgc  <  1:598330/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGC  >  CP009273/1311236‑1311534

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: