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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,441,8540CA16.5% 909.7 / 107.6 289R27S (CGT→AGT) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (146/95);  new base A (24/24);  total (170/119)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

AGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCC  >  CP009273/1441706‑1442003
                                                                                                                                                    |                                                                                                                                                     
aGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGgtcg                                                                                                                                                      >  1:1684918/1‑150 (MQ=255)
   gAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTcc                                                                                                                                                   >  2:4337895/1‑150 (MQ=255)
        aaGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGcc                                                                                                                                              >  2:4070695/1‑150 (MQ=255)
          gTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGcccc                                                                                                                                            <  2:1224764/150‑1 (MQ=255)
            gACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGc                                                                                                                                          >  1:837918/1‑150 (MQ=255)
            gACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGc                                                                                                                                          >  1:2655406/1‑150 (MQ=255)
             aCTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGca                                                                                                                                         >  1:3538604/1‑150 (MQ=255)
             aCTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAAGAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGca                                                                                                                                         >  1:1165949/1‑150 (MQ=255)
                tttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACaa                                                                                                                                      >  2:2225749/1‑150 (MQ=255)
                tttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACaa                                                                                                                                      >  1:3942092/1‑150 (MQ=255)
                 ttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACaaa                                                                                                                                     >  1:1001504/1‑150 (MQ=255)
                 ttCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACaaa                                                                                                                                     >  2:3694325/1‑150 (MQ=255)
                   cTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGt                                                                                                                                   >  2:2724260/1‑150 (MQ=255)
                   cTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGt                                                                                                                                   >  1:4099599/1‑150 (MQ=255)
                      gTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGaa                                                                                                                                >  1:2548229/1‑150 (MQ=255)
                        cGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGaaaa                                                                                                                              >  2:4000797/1‑150 (MQ=255)
                               tACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTg                                                                                                                       <  2:733607/150‑1 (MQ=255)
                                 cacTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGc                                                                                                                     >  1:302779/1‑150 (MQ=255)
                                       aaTCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGACCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTAGGCGATCtt                                                                                                               >  2:656031/1‑150 (MQ=255)
                                        aTCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATTGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCttt                                                                                                              >  1:1415098/1‑150 (MQ=255)
                                         tCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTg                                                                                                             <  2:2552135/150‑1 (MQ=255)
                                           acacTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGAt                                                                                                           >  1:2484559/1‑150 (MQ=255)
                                               tGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCgg                                                                                                       >  2:4449596/1‑150 (MQ=255)
                                               tGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCgg                                                                                                       <  1:3156786/150‑1 (MQ=255)
                                                  gTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCa                                                                                                    >  1:1735525/1‑150 (MQ=255)
                                                  gTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCa                                                                                                    >  1:1838225/1‑150 (MQ=255)
                                                  gTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCa                                                                                                    >  1:4164915/1‑150 (MQ=255)
                                                     aataatAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCg                                                                                                 >  2:277067/1‑150 (MQ=255)
                                                     aataatAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCg                                                                                                 >  2:622634/1‑150 (MQ=255)
                                                       taataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCggg                                                                                               >  2:3128817/1‑150 (MQ=255)
                                                        aataaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGc                                                                                              >  1:3349747/1‑150 (MQ=255)
                                                           aaGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAg                                                                                           <  2:883597/150‑1 (MQ=255)
                                                            aGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGa                                                                                          <  2:2641280/150‑1 (MQ=255)
                                                            aGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGa                                                                                          >  2:380550/1‑150 (MQ=255)
                                                            aGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGa                                                                                          >  2:2092374/1‑150 (MQ=255)
                                                            aGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGa                                                                                          <  1:4125104/150‑1 (MQ=255)
                                                            aGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGa                                                                                          >  1:1220028/1‑150 (MQ=255)
                                                             ggAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGaa                                                                                         <  1:2777934/150‑1 (MQ=255)
                                                                aaaGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAAtt                                                                                      >  2:1859456/1‑150 (MQ=255)
                                                                  aGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGc                                                                                    >  1:371349/1‑150 (MQ=255)
                                                                    tgatgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGt                                                                                  >  1:3502705/1‑150 (MQ=255)
                                                                      atgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCt                                                                                <  2:4588989/150‑1 (MQ=255)
                                                                      atgaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCt                                                                                <  2:925179/150‑1 (MQ=255)
                                                                        gaCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTAc                                                                              >  2:3300906/1‑150 (MQ=255)
                                                                         aCAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTAct                                                                             <  1:2246533/150‑1 (MQ=255)
                                                                          aagagCCGAATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGTCCCCGCACAAAGTGAAAAACTGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActg                                                                            <  1:4356868/149‑1 (MQ=255)
                                                                           agagCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgc                                                                           >  1:2990832/1‑150 (MQ=255)
                                                                           agagCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgc                                                                           >  1:3249225/1‑150 (MQ=255)
                                                                           agagCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgc                                                                           >  1:4253815/1‑150 (MQ=255)
                                                                            gagCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgct                                                                          >  2:2026548/1‑150 (MQ=255)
                                                                             agCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgctg                                                                         >  2:1082403/1‑150 (MQ=255)
                                                                             agCCACATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTActgctg                                                                         <  2:1455678/150‑1 (MQ=255)
                                                                               ccGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGAt                                                                       <  2:142118/150‑1 (MQ=255)
                                                                               ccGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGAt                                                                       <  1:496669/150‑1 (MQ=255)
                                                                                cGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATg                                                                      <  2:1678534/150‑1 (MQ=255)
                                                                                  cATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGcc                                                                    <  1:1975398/150‑1 (MQ=255)
                                                                                   aTGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCa                                                                   <  1:319172/150‑1 (MQ=255)
                                                                                   aTGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCa                                                                   >  1:2026948/1‑150 (MQ=255)
                                                                                     gTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAAc                                                                 <  2:1257813/150‑1 (MQ=255)
                                                                                      tAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACg                                                                >  1:455917/1‑150 (MQ=255)
                                                                                         cAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAg                                                             >  2:2871531/1‑150 (MQ=255)
                                                                                         cAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAg                                                             <  1:3917551/150‑1 (MQ=255)
                                                                                          aGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGc                                                            >  1:3686305/1‑150 (MQ=255)
                                                                                           gTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCg                                                           <  1:3666524/150‑1 (MQ=255)
                                                                                              ttAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGAt                                                        <  2:4312003/150‑1 (MQ=255)
                                                                                              ttAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGAt                                                        >  2:703259/1‑150 (MQ=255)
                                                                                               tAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATg                                                       >  2:4169902/1‑150 (MQ=255)
                                                                                                aaGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGt                                                      >  2:2597012/1‑150 (MQ=255)
                                                                                                   ccAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGa                                                   >  1:3458416/1‑150 (MQ=255)
                                                                                                      ggTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATaa                                                <  2:2126732/150‑1 (MQ=255)
                                                                                                        tCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACg                                              <  2:377066/150‑1 (MQ=255)
                                                                                                        tCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACg                                              <  1:1977248/150‑1 (MQ=255)
                                                                                                          cAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCa                                            >  2:49824/1‑150 (MQ=255)
                                                                                                           aaCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAg                                           <  1:3426520/150‑1 (MQ=255)
                                                                                                            aCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGt                                          <  1:1637423/150‑1 (MQ=255)
                                                                                                              aGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCa                                        >  2:3324284/1‑150 (MQ=255)
                                                                                                               gTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCAt                                       >  2:2475401/1‑150 (MQ=255)
                                                                                                               gTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCAt                                       >  2:2325140/1‑150 (MQ=255)
                                                                                                                 ttCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGt                                     >  1:3808540/1‑150 (MQ=255)
                                                                                                                 ttCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGt                                     <  2:2161960/150‑1 (MQ=255)
                                                                                                                    tcGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTg                                  >  2:3195329/1‑150 (MQ=255)
                                                                                                                      gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                                <  1:1459312/150‑1 (MQ=255)
                                                                                                                      gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                                >  2:2546488/1‑150 (MQ=255)
                                                                                                                      gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                                >  2:114225/1‑150 (MQ=255)
                                                                                                                      gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                                <  2:2794684/150‑1 (MQ=255)
                                                                                                                      gATCGTCAACACGGTCTTTATATTGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGcc                                >  2:1637109/1‑150 (MQ=255)
                                                                                                                       aTCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCt                               >  1:748867/1‑150 (MQ=255)
                                                                                                                       aTCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCt                               >  1:928257/1‑150 (MQ=255)
                                                                                                                           tCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                           >  1:95945/1‑150 (MQ=255)
                                                                                                                           tCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                           >  1:2518464/1‑150 (MQ=255)
                                                                                                                           tCAACACGGTCTGTATATTGATGGTCGTCCTGGCCCCGCACAAAGAGAAAAACGGTTGGCGATCTTTGAACCGGCCACCTGGAAAGAAATTGCGTCTACTGCTGATGCCAACGAAGAGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                           >  1:877178/1‑150 (MQ=255)
                                                                                                                             aaCACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCggg                         <  2:715855/150‑1 (MQ=255)
                                                                                                                                      tttATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGtct                <  1:4014389/150‑1 (MQ=255)
                                                                                                                                        tataTTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcg              <  1:4439203/150‑1 (MQ=255)
                                                                                                                                        tataTTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcg              >  2:1791044/1‑150 (MQ=255)
                                                                                                                                          tatTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcgcg            >  1:2320697/1‑150 (MQ=255)
                                                                                                                                           atTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGt           <  2:391080/150‑1 (MQ=255)
                                                                                                                                             tGATGGTAGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCg         <  2:1096862/150‑1 (MQ=255)
                                                                                                                                                  gtcgtcCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTggg    >  2:3231618/1‑150 (MQ=255)
                                                                                                                                                    cgtcCTGGCCACGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGcc  <  1:3122771/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCGAAAAAAGTGACTTTTCTTGTCGCTGCGTACACTGAAATCACACTGGGTAAATAATAAGGAAAAGTGATGACAGAGCCGCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCC  >  CP009273/1441706‑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: