Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,441,854 C→A 8.3% R27S (CGT→AGT)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,441,8540CA 8.3% 1207.7 / 53.0 320R27S (CGT→AGT) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (144/149);  new base A (11/16);  total (155/165)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

CAGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCC  >  CP009273/1441705‑1442003
                                                                                                                                                     |                                                                                                                                                     
caGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGgtc                                                                                                                                                       <  1:1547805/150‑1 (MQ=255)
 aGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGgtcg                                                                                                                                                      <  2:1739555/150‑1 (MQ=255)
  gCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGgtagt                                                                                                                                                     <  2:2008603/150‑1 (MQ=255)
   cGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGgtcgtc                                                                                                                                                    >  1:4107094/1‑150 (MQ=255)
    gAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTcc                                                                                                                                                   >  1:1040053/1‑150 (MQ=255)
    gAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTcc                                                                                                                                                   >  1:861946/1‑150 (MQ=255)
       aaaaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTgg                                                                                                                                                >  2:2736845/1‑150 (MQ=255)
         aaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGcc                                                                                                                                              >  1:618165/1‑150 (MQ=255)
          aGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGccc                                                                                                                                             <  2:549360/150‑1 (MQ=255)
             gACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGc                                                                                                                                          >  1:263610/1‑150 (MQ=255)
              aCTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGca                                                                                                                                         >  1:682491/1‑150 (MQ=255)
               cTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGcac                                                                                                                                        <  1:2280561/150‑1 (MQ=255)
                 tttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACaa                                                                                                                                      <  1:3753342/150‑1 (MQ=255)
                 tttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACaa                                                                                                                                      >  2:5043890/1‑150 (MQ=255)
                  ttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACaaa                                                                                                                                     >  1:1688743/1‑150 (MQ=255)
                      tGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGa                                                                                                                                 <  1:1361018/150‑1 (MQ=255)
                         cGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGaaaa                                                                                                                              <  1:4743687/150‑1 (MQ=255)
                            tGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACg                                                                                                                           >  1:463194/1‑150 (MQ=255)
                                 acacTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTgg                                                                                                                      <  1:761864/150‑1 (MQ=255)
                                    cTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGa                                                                                                                   <  1:1442375/150‑1 (MQ=255)
                                     tGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGAt                                                                                                                  >  2:847419/1‑150 (MQ=255)
                                        aaTCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCtt                                                                                                               <  1:3985385/150‑1 (MQ=255)
                                          tCACACTGGGTAAATAATAATTAATAGTGATGACAGAGCCGAATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTAGATCGTCAACACGGTATTTATATTGATGGTCGTCCTGGCCCCGCACAAACTGAAAAACGGTTGGAGATCTTTg                                                                                                             <  2:1156418/150‑1 (MQ=255)
                                            acacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCTATCTTTGAt                                                                                                           >  2:4790189/1‑150 (MQ=255)
                                            acacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGAt                                                                                                           >  1:2758036/1‑150 (MQ=255)
                                            acacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGAt                                                                                                           >  1:2120375/1‑150 (MQ=255)
                                             cacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATc                                                                                                          <  2:2899544/150‑1 (MQ=255)
                                             cacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATc                                                                                                          >  1:2157664/1‑150 (MQ=255)
                                                tGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACACAGTGAAAAACGGTTGGCGATCTTTGATCCgg                                                                                                       >  2:3138507/1‑150 (MQ=255)
                                                tGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCgg                                                                                                       >  1:4186255/1‑150 (MQ=255)
                                                     aaaTAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCAcc                                                                                                  <  2:2087481/150‑1 (MQ=255)
                                                       ataataAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCgg                                                                                                <  2:442169/150‑1 (MQ=255)
                                                          ataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCa                                                                                             >  1:3912987/1‑150 (MQ=255)
                                                             aGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGa                                                                                          >  2:28152/1‑150 (MQ=255)
                                                              ggAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGaa                                                                                         <  2:463194/150‑1 (MQ=255)
                                                              ggAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGaa                                                                                         <  1:205491/150‑1 (MQ=255)
                                                                 aaaGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAAtt                                                                                      <  1:2344614/150‑1 (MQ=255)
                                                                  aaGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTg                                                                                     >  1:1209841/1‑150 (MQ=255)
                                                                        tgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTa                                                                               >  2:430139/1‑150 (MQ=255)
                                                                        tgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTa                                                                               >  1:1431278/1‑150 (MQ=255)
                                                                           cAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActg                                                                            >  2:3402953/1‑150 (MQ=255)
                                                                           cAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActg                                                                            >  2:2984863/1‑150 (MQ=255)
                                                                             gagCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgct                                                                          >  2:45914/1‑150 (MQ=255)
                                                                              agCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgctg                                                                         <  2:1945246/150‑1 (MQ=255)
                                                                               gCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGa                                                                        <  2:841504/150‑1 (MQ=255)
                                                                               gCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGa                                                                        >  2:2422854/1‑150 (MQ=255)
                                                                               gCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGa                                                                        <  2:3908153/150‑1 (MQ=255)
                                                                                ccGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACAGCTGAt                                                                       >  2:317124/1‑150 (MQ=255)
                                                                                 cGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATg                                                                      >  2:2054687/1‑150 (MQ=255)
                                                                                   cATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGcc                                                                    <  1:2792645/150‑1 (MQ=255)
                                                                                     tGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCaa                                                                  <  2:5087409/150‑1 (MQ=255)
                                                                                      gTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAAc                                                                 <  2:2292586/150‑1 (MQ=255)
                                                                                        aGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGa                                                               <  1:590613/150‑1 (MQ=255)
                                                                                        aGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGa                                                               >  1:921067/1‑150 (MQ=255)
                                                                                         gCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGaa                                                              >  2:3697201/1‑150 (MQ=255)
                                                                                         gCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGaa                                                              >  1:4994399/1‑150 (MQ=255)
                                                                                         gCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGaa                                                              <  1:4269943/150‑1 (MQ=255)
                                                                                           aGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGc                                                            <  1:745965/150‑1 (MQ=255)
                                                                                           aGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGc                                                            <  1:2918744/150‑1 (MQ=255)
                                                                                             tATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCgg                                                          <  1:5052583/150‑1 (MQ=255)
                                                                                               ttAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGAt                                                        <  2:3278723/150‑1 (MQ=255)
                                                                                                tAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCACCGAAGCGGATg                                                       >  2:4081345/1‑150 (MQ=255)
                                                                                                tAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGagg                                                       >  1:3260052/1‑148 (MQ=255)
                                                                                                  aGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTa                                                     >  2:3832226/1‑150 (MQ=255)
                                                                                                  aGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTa                                                     <  1:304222/150‑1 (MQ=255)
                                                                                                   gCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAg                                                    <  2:374193/150‑1 (MQ=255)
                                                                                                          ccAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGc                                             >  2:3286375/1‑150 (MQ=255)
                                                                                                           cAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCa                                            >  2:3529159/1‑150 (MQ=255)
                                                                                                            aaCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAg                                           <  1:3131401/150‑1 (MQ=255)
                                                                                                              cAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTc                                         >  1:2855426/1‑150 (MQ=255)
                                                                                                                 tttCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATg                                      >  2:1801986/1‑150 (MQ=255)
                                                                                                                 tttCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATg                                      >  2:1037269/1‑150 (MQ=255)
                                                                                                                    ctcGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCt                                   <  2:4133120/150‑1 (MQ=255)
                                                                                                                    ctcGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCt                                   <  1:4452368/150‑1 (MQ=255)
                                                                                                                     tcGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTg                                  >  1:4038598/1‑150 (MQ=255)
                                                                                                                     tcGATCGTCAACACGGTCTTTATATTGATGGTAGTCATGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTg                                  >  1:2564840/1‑150 (MQ=255)
                                                                                                                        aTCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCt                               >  2:3057924/1‑150 (MQ=255)
                                                                                                                        aTCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCt                               <  1:2432465/150‑1 (MQ=255)
                                                                                                                        aTCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCt                               >  2:4591622/1‑150 (MQ=255)
                                                                                                                        aTCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCt                               <  2:1561981/150‑1 (MQ=255)
                                                                                                                          cGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTgg                             <  1:1562554/150‑1 (MQ=255)
                                                                                                                           gTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGc                            >  2:1897408/1‑150 (MQ=255)
                                                                                                                            tCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                           <  1:1967418/150‑1 (MQ=255)
                                                                                                                             cAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCgg                          >  2:2825161/1‑150 (MQ=255)
                                                                                                                             cAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCgg                          <  1:4923775/150‑1 (MQ=255)
                                                                                                                              aaCACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCggg                         >  2:4186701/1‑150 (MQ=255)
                                                                                                                               acacGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAACGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGc                        <  2:224515/150‑1 (MQ=255)
                                                                                                                                cacGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGcc                       <  2:2604874/150‑1 (MQ=255)
                                                                                                                                 acGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGCTGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCt                      >  2:3104509/1‑150 (MQ=255)
                                                                                                                                 acGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCt                      >  1:2581558/1‑150 (MQ=255)
                                                                                                                                  cGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCtt                     <  2:1118201/150‑1 (MQ=255)
                                                                                                                                  cGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCtt                     >  2:4681468/1‑150 (MQ=255)
                                                                                                                                       tttATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGtct                >  1:3612916/1‑150 (MQ=255)
                                                                                                                                         tataTTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcg              <  1:1801986/150‑1 (MQ=255)
                                                                                                                                           tatTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcgcg            >  2:3978254/1‑150 (MQ=255)
                                                                                                                                             ttGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTc          >  1:1391931/1‑150 (MQ=255)
                                                                                                                                             ttGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTc          >  2:436976/1‑150 (MQ=255)
                                                                                                                                               gATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGc        <  1:3382149/150‑1 (MQ=255)
                                                                                                                                               gATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGc        >  1:172323/1‑150 (MQ=255)
                                                                                                                                                aTGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCt       >  1:4553266/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCC  >  CP009273/1441705‑1442003

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: