Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,456,739 G→T 6.0% D153Y (GAT→TAT)  gspK → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,456,7390GT 6.0% 1120.5 / 27.4 300D153Y (GAT→TAT) gspKgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (153/128);  new base T (12/7);  total (165/135)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

GTTTTAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTT  >  CP009273/3456590‑3456888
                                                                                                                                                     |                                                                                                                                                     
gTTTTAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTg                                                                                                                                                       >  1:610966/1‑150 (MQ=255)
  tttAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGAt                                                                                                                                                     >  2:1732171/1‑150 (MQ=255)
   ttAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATg                                                                                                                                                    <  1:3296464/150‑1 (MQ=255)
       tctcAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGa                                                                                                                                                >  2:3347444/1‑150 (MQ=255)
       tctcAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGa                                                                                                                                                >  2:1073727/1‑150 (MQ=255)
        ctcAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGAt                                                                                                                                               >  2:26935/1‑150 (MQ=255)
        ctcAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGAt                                                                                                                                               >  2:27332/1‑150 (MQ=255)
           aaTGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTTATGTCGATGAc                                                                                                                                            >  2:2201188/1‑150 (MQ=255)
               cGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGAATAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACAAAATTGATGTCGATGACTcac                                                                                                                                        >  2:558018/1‑150 (MQ=255)
                  tAGCAAAAATCTCTGACTACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCATTCGCTGACTATATTGATGTCGATGACTcaccac                                                                                                                                     <  1:4465845/150‑1 (MQ=255)
                  tAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTcaccac                                                                                                                                     >  1:2783185/1‑150 (MQ=255)
                   aGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACg                                                                                                                                    >  2:1225877/1‑150 (MQ=255)
                     cAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCt                                                                                                                                  <  2:3043798/150‑1 (MQ=255)
                        aaaTCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTc                                                                                                                               <  1:1523198/150‑1 (MQ=255)
                          aTCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCAc                                                                                                                             <  2:2724446/150‑1 (MQ=255)
                           tctcTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACg                                                                                                                            <  1:1638210/150‑1 (MQ=255)
                           tctcTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACg                                                                                                                            >  1:2395275/1‑150 (MQ=255)
                             tctGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATCTTGATGTCGATGACTCACCACGCTTTCACGgc                                                                                                                          >  2:3799729/1‑150 (MQ=255)
                             tctGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGgc                                                                                                                          >  2:3799542/1‑150 (MQ=255)
                              ctGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGgcg                                                                                                                         >  1:345803/1‑150 (MQ=255)
                              ctGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGgcg                                                                                                                         >  1:2202922/1‑150 (MQ=255)
                                 acgacCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCag                                                                                                                      <  2:3672625/150‑1 (MQ=255)
                                  cgacCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCaga                                                                                                                     <  1:1561227/150‑1 (MQ=255)
                                   gacCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGGTTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCagaa                                                                                                                    >  2:3738318/1‑150 (MQ=255)
                                   gacCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCagaa                                                                                                                    >  2:3576932/1‑150 (MQ=255)
                                   gacCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCagaa                                                                                                                    >  1:2284274/1‑150 (MQ=255)
                                    acCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCagaag                                                                                                                   <  1:4293877/150‑1 (MQ=255)
                                       cGCTGGCTTCCCCTGAGTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATa                                                                                                                <  2:3447039/150‑1 (MQ=255)
                                        gCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAg                                                                                                               >  1:1463895/1‑150 (MQ=255)
                                         cTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGc                                                                                                              >  2:4080384/1‑150 (MQ=255)
                                              ttCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTa                                                                                                         <  2:4425309/150‑1 (MQ=255)
                                              ttCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTa                                                                                                         >  1:2229069/1‑150 (MQ=255)
                                                     gATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGc                                                                                                  >  2:1289709/1‑150 (MQ=255)
                                                      aTTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGcc                                                                                                 >  2:256192/1‑150 (MQ=255)
                                                      aTTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGcc                                                                                                 >  2:2850780/1‑150 (MQ=255)
                                                      aTTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGcc                                                                                                 <  2:868130/150‑1 (MQ=255)
                                                          ccccGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGAc                                                                                             >  2:4447874/1‑150 (MQ=255)
                                                          ccccGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGAc                                                                                             <  1:1883261/150‑1 (MQ=255)
                                                           cccGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGAcg                                                                                            <  1:2961344/150‑1 (MQ=255)
                                                           cccGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGAcg                                                                                            <  2:1850798/150‑1 (MQ=255)
                                                           cccGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGAcg                                                                                            >  2:3984787/1‑150 (MQ=255)
                                                                  aGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGcccc                                                                                     <  1:222960/150‑1 (MQ=255)
                                                                      aTTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCa                                                                                 >  1:4082370/1‑150 (MQ=255)
                                                                      aTTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCa                                                                                 <  1:4383901/150‑1 (MQ=255)
                                                                      aTTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCa                                                                                 >  2:2874710/1‑150 (MQ=255)
                                                                      aTTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCa                                                                                 >  1:2157139/1‑150 (MQ=255)
                                                                        tCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATa                                                                               >  2:2715268/1‑150 (MQ=255)
                                                                        tCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATa                                                                               <  2:1179874/150‑1 (MQ=255)
                                                                         cAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAg                                                                              <  2:461399/150‑1 (MQ=255)
                                                                            tGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTTATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGc                                                                           <  1:3904423/150‑1 (MQ=255)
                                                                                ctgctgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAAt                                                                       >  1:2397636/1‑150 (MQ=255)
                                                                                 tgctgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTTATGACGATGACTCACCACGCTTTCACGGCGAAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATc                                                                      >  2:1755777/1‑150 (MQ=255)
                                                                                  gctgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCa                                                                     <  1:2437935/150‑1 (MQ=255)
                                                                                     gATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTTATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGAt                                                                  <  1:3603985/150‑1 (MQ=255)
                                                                                     gATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGAt                                                                  >  1:1167425/1‑150 (MQ=255)
                                                                                     gATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGAt                                                                  <  1:3139679/150‑1 (MQ=255)
                                                                                      aTTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATg                                                                 >  1:2005657/1‑150 (MQ=255)
                                                                                       ttAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGc                                                                <  2:2955073/150‑1 (MQ=255)
                                                                                        tAAAGCCTGTATCGATCGCGGAAACACCGATGAGATCTTTCAATCAATCGCTGACTATATTTATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCCGATGCt                                                               <  1:2103274/150‑1 (MQ=255)
                                                                                          aCGCCGGTATCGATCGCGGAAACACCGATTAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCtat                                                             >  2:1907848/1‑148 (MQ=255)
                                                                                              cGGTATCGATCGCGGAAACACCGATGAGATCGTGCTATCAATCGCTGACTTTATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGTTGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCt                                                         <  2:3970085/150‑1 (MQ=255)
                                                                                              cGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCt                                                         >  1:4499114/1‑150 (MQ=255)
                                                                                              cGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCt                                                         <  2:3499383/150‑1 (MQ=255)
                                                                                               ggTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTg                                                        >  2:4032293/1‑150 (MQ=255)
                                                                                               ggTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTg                                                        <  2:429667/150‑1 (MQ=255)
                                                                                                 tATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGAc                                                      >  1:2915035/1‑150 (MQ=255)
                                                                                                   tcgatcgCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTg                                                    <  1:1910033/150‑1 (MQ=255)
                                                                                                         gcgGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTTATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAAt                                              <  2:3255270/150‑1 (MQ=255)
                                                                                                          cgGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAAtt                                             >  1:1080282/1‑150 (MQ=255)
                                                                                                          cgGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAAtt                                             >  1:259817/1‑150 (MQ=255)
                                                                                                           ggAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTa                                            >  1:4233230/1‑150 (MQ=255)
                                                                                                                   cGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGAt                                    >  2:4170986/1‑150 (MQ=255)
                                                                                                                   cGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGAt                                    <  1:2839448/150‑1 (MQ=255)
                                                                                                                    gATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATa                                   <  1:3961858/150‑1 (MQ=255)
                                                                                                                      tGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATaaa                                 <  1:465137/150‑1 (MQ=255)
                                                                                                                        agaTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAg                               <  1:851773/150‑1 (MQ=255)
                                                                                                                        agaTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAg                               <  1:3721325/150‑1 (MQ=255)
                                                                                                                         gaTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAgg                              <  1:1073633/150‑1 (MQ=255)
                                                                                                                         gaTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAgg                              <  1:3347444/150‑1 (MQ=255)
                                                                                                                            cGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCAt                           <  2:4097127/150‑1 (MQ=255)
                                                                                                                              tGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATaa                         >  1:2074850/1‑150 (MQ=255)
                                                                                                                               gCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAAc                        >  2:100209/1‑150 (MQ=255)
                                                                                                                               gCAATCAATCGCTGACTATATTGATGTCGATGAATCACCACGATTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATaaa                        >  2:4325500/1‑149 (MQ=255)
                                                                                                                                caatcaatCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAAcc                       >  1:3758115/1‑150 (MQ=255)
                                                                                                                                caatcaatCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAAcc                       <  1:698505/150‑1 (MQ=255)
                                                                                                                                  atcaatcGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGa                     <  2:2322972/150‑1 (MQ=255)
                                                                                                                                   tcaatcGCTGACTATATTTATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGaa                    >  2:249483/1‑150 (MQ=255)
                                                                                                                                   tcaatcGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGaa                    >  1:3199235/1‑150 (MQ=255)
                                                                                                                                   tcaatcGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGaa                    <  2:1340652/150‑1 (MQ=255)
                                                                                                                                        cGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACAt               <  2:1056704/150‑1 (MQ=255)
                                                                                                                                         gCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATc              <  2:2783185/150‑1 (MQ=255)
                                                                                                                                         gCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATc              <  2:3893106/150‑1 (MQ=255)
                                                                                                                                          cTGACTATATTTATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCt             >  1:3248394/1‑150 (MQ=255)
                                                                                                                                            gACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGACAACATCTAt           >  2:3468150/1‑150 (MQ=255)
                                                                                                                                             aCTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATc          >  1:1693770/1‑150 (MQ=255)
                                                                                                                                              cTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCa         >  1:184879/1‑150 (MQ=255)
                                                                                                                                                atatTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGc       >  2:1970336/1‑150 (MQ=255)
                                                                                                                                                atatTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGc       <  1:1289709/150‑1 (MQ=255)
                                                                                                                                                  atTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCgg     >  1:1776063/1‑150 (MQ=255)
                                                                                                                                                     gATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCtt  <  2:883875/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTTAATCTCAATGCGTTAGCAAAAATCTCTGACGACCCGCTGGCTTCCCCTGATTTCCCCGCCCAGGTATTCAGTGCGCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTT  >  CP009273/3456590‑3456888

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: