Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,122,453 C→T 94.8% A227T (GCT→ACT)  wcaE ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,122,4530CT94.8% 738.8 / 24.2 229A227T (GCT→ACT) wcaEputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base T (113/104);  total (119/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 229 total aligned reads displayed.

TTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTAT  >  CP009273/2122310‑2122596
                                                                                                                                               |                                                                                                                                               
   aaTAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAg                                                                                                                                          <  2:2815417/148‑1 (MQ=255)
   aaTAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAg                                                                                                                                          >  1:151945/1‑148 (MQ=255)
      aGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCa                                                                                                                                              <  1:2438977/141‑1 (MQ=255)
       gCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAg                                                                                                                                      <  2:1074675/148‑1 (MQ=255)
        cGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAgg                                                                                                                                     >  2:1960188/1‑148 (MQ=255)
         gTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGc                                                                                                                                    <  2:2143734/148‑1 (MQ=255)
         gTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGc                                                                                                                                    <  2:681039/148‑1 (MQ=255)
          ttGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGca                                                                                                                                   <  2:2976241/148‑1 (MQ=255)
          ttGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGca                                                                                                                                   <  1:210714/148‑1 (MQ=255)
           tGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGcac                                                                                                                                  <  2:2394208/148‑1 (MQ=255)
             ccgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACAt                                                                                                                                <  1:2056120/148‑1 (MQ=255)
               gccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGt                                                                                                                              <  2:3000770/148‑1 (MQ=255)
                ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTa                                                                                                                             <  2:3163774/148‑1 (MQ=255)
                ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTa                                                                                                                             <  2:3163704/148‑1 (MQ=255)
                 cccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTaa                                                                                                                            <  1:2770871/148‑1 (MQ=255)
                   cGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAata                                                                                                                          <  1:602528/148‑1 (MQ=255)
                     gAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATAtt                                                                                                                        >  1:1235293/1‑148 (MQ=255)
                      aaCCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATAttt                                                                                                                       >  2:1303099/1‑148 (MQ=255)
                        cccTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGt                                                                                                                     <  2:2924723/148‑1 (MQ=255)
                        cccTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGt                                                                                                                     >  2:566220/1‑148 (MQ=255)
                         ccTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTc                                                                                                                    <  1:1391249/148‑1 (MQ=255)
                          cTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCg                                                                                                                   <  2:2197275/148‑1 (MQ=255)
                          cTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCg                                                                                                                   >  1:1157028/1‑148 (MQ=255)
                           tttCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGt                                                                                                                  <  1:2902726/148‑1 (MQ=255)
                            ttCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGtt                                                                                                                 >  2:1320114/1‑148 (MQ=255)
                              cGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTgg                                                                                                               <  2:2734106/148‑1 (MQ=255)
                               ggCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGa                                                                                                              <  1:1956687/148‑1 (MQ=255)
                                gCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGAc                                                                                                             <  1:1118968/148‑1 (MQ=255)
                                gCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGAc                                                                                                             <  2:2444994/148‑1 (MQ=255)
                                 cACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACt                                                                                                            <  1:581164/148‑1 (MQ=255)
                                  aCCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACtt                                                                                                           <  2:870687/148‑1 (MQ=255)
                                   ccGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACttt                                                                                                          >  2:2300275/1‑148 (MQ=255)
                                    cGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACtttt                                                                                                         <  1:421023/148‑1 (MQ=255)
                                    cGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACtttt                                                                                                         <  1:18717/148‑1 (MQ=255)
                                       gaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTc                                                                                                      >  1:1042104/1‑148 (MQ=255)
                                       gaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTc                                                                                                      >  2:2780058/1‑148 (MQ=255)
                                       gaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTc                                                                                                      >  2:1082484/1‑148 (MQ=255)
                                         aaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcgc                                                                                                    <  2:704027/148‑1 (MQ=255)
                                          aCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcgcg                                                                                                   >  2:2005950/1‑148 (MQ=255)
                                                     cTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACaa                                                                                        <  1:2181418/148‑1 (MQ=255)
                                                          aTTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCa                                                                                   <  1:1305150/148‑1 (MQ=255)
                                                              gTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCAtatt                                                                               <  2:3195318/148‑1 (MQ=255)
                                                              gTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCAtatt                                                                               <  2:2154493/148‑1 (MQ=255)
                                                               ttttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCAtatta                                                                              >  1:984697/1‑148 (MQ=255)
                                                                tttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCAtattat                                                                             >  1:652168/1‑148 (MQ=255)
                                                                     ttATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggtg                                                                        >  2:2172839/1‑148 (MQ=255)
                                                                      tataTTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggtgg                                                                       >  1:2365372/1‑148 (MQ=255)
                                                                       atatTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTggtggt                                                                      <  2:1376744/148‑1 (MQ=255)
                                                                         atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAg                                                                    >  2:3123084/1‑148 (MQ=255)
                                                                         atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAg                                                                    <  1:1785491/148‑1 (MQ=255)
                                                                          ttCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGa                                                                   <  1:2151468/148‑1 (MQ=255)
                                                                           tCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGAt                                                                  >  2:1800926/1‑148 (MQ=255)
                                                                           tCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGAt                                                                  >  2:1460318/1‑148 (MQ=255)
                                                                             aGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAc                                                                <  2:3302561/148‑1 (MQ=255)
                                                                               aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAccc                                                              >  1:1599930/1‑148 (MQ=255)
                                                                               aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAccc                                                              >  1:352600/1‑148 (MQ=255)
                                                                                  ttGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACcccac                                                           >  1:826257/1‑148 (MQ=255)
                                                                                  ttGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGATGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTAATGGTGGTAGATACcccac                                                           >  1:2139103/1‑148 (MQ=255)
                                                                                    gTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACcccaccc                                                         <  2:2524322/148‑1 (MQ=255)
                                                                                      tatataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCAt                                                       >  2:2669322/1‑148 (MQ=255)
                                                                                      tatataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCAt                                                       >  1:797217/1‑148 (MQ=255)
                                                                                         ataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGa                                                    >  2:773519/1‑148 (MQ=255)
                                                                                          taAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaa                                                   >  1:531056/1‑148 (MQ=255)
                                                                                           aaGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaa                                                  <  1:1499825/148‑1 (MQ=255)
                                                                                           aaGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaa                                                  >  1:2875533/1‑148 (MQ=255)
                                                                                           aaGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaa                                                  <  1:50985/148‑1 (MQ=255)
                                                                                            aGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaaa                                                 >  2:400731/1‑148 (MQ=255)
                                                                                                 ttCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAg                                            <  2:1457272/148‑1 (MQ=255)
                                                                                                      cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                       >  2:1753795/1‑148 (MQ=255)
                                                                                                      cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                       >  2:395809/1‑148 (MQ=255)
                                                                                                      cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                       <  2:490352/148‑1 (MQ=255)
                                                                                                      cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                       <  2:522266/148‑1 (MQ=255)
                                                                                                          gAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGc                                   <  1:3123084/148‑1 (MQ=255)
                                                                                                            ggTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                 <  1:22235/148‑1 (MQ=255)
                                                                                                            ggTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                 <  1:308420/148‑1 (MQ=255)
                                                                                                            ggTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                 >  2:357705/1‑148 (MQ=255)
                                                                                                            ggTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                 <  2:1629513/148‑1 (MQ=255)
                                                                                                            ggTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                 <  1:1663030/148‑1 (MQ=255)
                                                                                                             gtagtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAt                                <  2:1125625/148‑1 (MQ=255)
                                                                                                              tagtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCCGGCCAtt                               >  2:1133366/1‑148 (MQ=255)
                                                                                                               agtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTg                              <  2:1002546/148‑1 (MQ=255)
                                                                                                                 taCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAg                            >  2:3275714/1‑148 (MQ=255)
                                                                                                                     ttGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGtttt                        <  1:3129710/148‑1 (MQ=255)
                                                                                                                          gTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTaaa                   <  2:1599930/148‑1 (MQ=255)
                                                                                                                           tAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAAt                  <  2:1668255/148‑1 (MQ=255)
                                                                                                                              aTGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCa               <  1:484177/148‑1 (MQ=255)
                                                                                                                               tGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCAt              >  2:1157557/1‑148 (MQ=255)
                                                                                                                                gCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATa             >  1:3047002/1‑148 (MQ=255)
                                                                                                                                gCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATa             >  1:925617/1‑148 (MQ=255)
                                                                                                                                  cAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCCGGCCATTGAGTTTTTTAAATGCATAAc           >  1:910296/1‑148 (MQ=255)
                                                                                                                                  cAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAAc           <  1:2852631/148‑1 (MQ=255)
                                                                                                                                     gATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCgg        >  1:3271341/1‑148 (MQ=255)
                                                                                                                                      aTAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGc       <  2:1235293/148‑1 (MQ=255)
                                                                                                                                       tAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCt      <  2:2875533/148‑1 (MQ=255)
                                                                                                                                       tAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCt      <  2:2784847/148‑1 (MQ=255)
                                                                                                                                         aTTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCttt    <  2:470949/148‑1 (MQ=255)
                                                                                                                                         aTTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCttt    >  1:266826/1‑148 (MQ=255)
                                                                                                                                         aTTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCttt    >  1:2733534/1‑148 (MQ=255)
                                                                                                                                         aTTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCttt    >  1:2734187/1‑148 (MQ=255)
                                                                                                                                           tCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtat  >  2:1847270/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTAT  >  CP009273/2122310‑2122596

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: