Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,122,453 C→T 100% 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,4530CT97.1% 903.9 / 10.7 275A227T (GCT→ACT) wcaEputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (1/7);  new base T (144/123);  total (145/130)
Fisher's exact test for biased strand distribution p-value = 2.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 275 total aligned reads displayed.

TTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACA  >  CP009273/2122310‑2122599
                                                                                                                                               |                                                                                                                                                  
 ttAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAg                                                                                                                                             <  2:226651/150‑1 (MQ=255)
  tAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGc                                                                                                                                            >  1:621080/1‑150 (MQ=255)
    aTAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCa                                                                                                                                          <  2:4137404/150‑1 (MQ=255)
    aTAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCa                                                                                                                                          <  2:3133329/150‑1 (MQ=255)
     tAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAg                                                                                                                                         >  2:4934072/1‑150 (MQ=255)
        cGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGca                                                                                                                                      >  1:4603555/1‑150 (MQ=255)
        cGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGca                                                                                                                                      <  1:536566/150‑1 (MQ=255)
         gTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGcac                                                                                                                                     >  2:1015417/1‑150 (MQ=255)
         gTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGcac                                                                                                                                     >  2:1481444/1‑150 (MQ=255)
          ttGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGcaca                                                                                                                                    >  2:4367315/1‑150 (MQ=255)
          ttGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGcaca                                                                                                                                    <  1:4295264/150‑1 (MQ=255)
            gccgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATg                                                                                                                                  <  1:3159118/150‑1 (MQ=255)
              cgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTa                                                                                                                                <  2:1862315/150‑1 (MQ=255)
              cgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTa                                                                                                                                <  2:1760005/150‑1 (MQ=255)
              cgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTa                                                                                                                                <  2:3554113/150‑1 (MQ=255)
               gccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTaa                                                                                                                               >  2:288260/1‑150 (MQ=255)
                ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAat                                                                                                                              <  1:1517526/150‑1 (MQ=255)
                 cccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAata                                                                                                                             >  2:4729505/1‑150 (MQ=255)
                       aCCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTc                                                                                                                       <  1:1480401/150‑1 (MQ=255)
                         ccTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGt                                                                                                                     <  1:4942210/150‑1 (MQ=255)
                            ttCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATTTTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTgg                                                                                                                  <  1:725789/150‑1 (MQ=255)
                               ggCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTGTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACt                                                                                                               >  2:1299995/1‑150 (MQ=255)
                                gCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACtt                                                                                                              <  1:1590453/150‑1 (MQ=255)
                                     gagaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCATACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATCATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTc                                                                                                         >  1:2977667/1‑150 (MQ=255)
                                      agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcg                                                                                                        <  1:2621418/150‑1 (MQ=255)
                                      agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTGTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCCGGATAATTCAGTCCAAAAGCCAGGCCCATGTAATATTTGTCGATGGACTTTTTTcg                                                                                                        >  2:4344756/1‑150 (MQ=255)
                                       gaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcgc                                                                                                       >  1:3562799/1‑150 (MQ=255)
                                        aaaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcgcg                                                                                                      >  2:4472724/1‑150 (MQ=255)
                                        aaaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcgcg                                                                                                      <  1:4499818/150‑1 (MQ=255)
                                         aaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGt                                                                                                     <  2:1703044/150‑1 (MQ=255)
                                         aaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGt                                                                                                     >  2:3828832/1‑150 (MQ=255)
                                          aCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTc                                                                                                    <  1:2600576/150‑1 (MQ=255)
                                          aCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTc                                                                                                    <  1:79416/150‑1 (MQ=255)
                                          aCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTc                                                                                                    <  1:2552189/150‑1 (MQ=255)
                                           tcGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCa                                                                                                   >  2:5204/2‑150 (MQ=255)
                                           ccGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCa                                                                                                   <  1:1679221/150‑1 (MQ=255)
                                           ccGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCa                                                                                                   <  1:3569717/150‑1 (MQ=255)
                                            cGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAg                                                                                                  >  1:2181501/1‑150 (MQ=255)
                                                 aaaTCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGcacac                                                                                             >  2:1076195/1‑150 (MQ=255)
                                                  aaTCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGcacaca                                                                                            >  1:4213106/1‑150 (MQ=255)
                                                    tCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAAt                                                                                          >  2:7468/1‑150 (MQ=255)
                                                         cATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCat                                                                                     >  2:2816332/1‑150 (MQ=255)
                                                         cATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCat                                                                                     >  1:3828856/1‑150 (MQ=255)
                                                          aTTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCata                                                                                    >  2:4442229/1‑150 (MQ=255)
                                                              gTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCAtattat                                                                                >  1:3325142/1‑150 (MQ=255)
                                                              gTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCAtattat                                                                                <  2:2216844/150‑1 (MQ=255)
                                                                tttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtg                                                                              <  1:4191387/150‑1 (MQ=255)
                                                                   ccTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggtg                                                                           >  1:4067085/1‑150 (MQ=255)
                                                                    cTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggtgg                                                                          <  2:1246952/150‑1 (MQ=255)
                                                                     ttATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTggtggt                                                                         >  1:4861561/1‑150 (MQ=255)
                                                                     ttATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTggtggt                                                                         <  2:2990665/150‑1 (MQ=255)
                                                                         atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTTATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGAt                                                                     <  1:1859729/150‑1 (MQ=255)
                                                                         atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGAt                                                                     >  2:3727568/1‑150 (MQ=255)
                                                                          ttCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATa                                                                    <  1:3599366/150‑1 (MQ=255)
                                                                            cAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAcc                                                                  >  2:1912770/1‑150 (MQ=255)
                                                                            cAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAcc                                                                  >  2:4271464/1‑150 (MQ=255)
                                                                               aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACccca                                                               <  1:678152/150‑1 (MQ=255)
                                                                               aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACccca                                                               <  2:3752130/150‑1 (MQ=255)
                                                                                  ttGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACcccaccc                                                            <  1:3890481/150‑1 (MQ=255)
                                                                                   tGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACcccaccca                                                           <  2:3912637/150‑1 (MQ=255)
                                                                                     ttATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATg                                                         <  2:2328120/150‑1 (MQ=255)
                                                                                       atataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGa                                                       >  2:1652850/1‑150 (MQ=255)
                                                                                       atataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGa                                                       <  1:2051943/150‑1 (MQ=255)
                                                                                       anataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCCTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGa                                                       >  1:3468115/1‑150 (MQ=255)
                                                                                          taAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaaa                                                    <  2:1619280/150‑1 (MQ=255)
                                                                                          taAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaaa                                                    <  2:4056556/150‑1 (MQ=255)
                                                                                            aGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAAtt                                                  <  2:1982271/150‑1 (MQ=255)
                                                                                             ggCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTc                                                 >  1:834101/1‑150 (MQ=255)
                                                                                             ggCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTc                                                 <  1:1167297/150‑1 (MQ=255)
                                                                                               cTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAg                                               >  2:1131231/1‑150 (MQ=255)
                                                                                                tttCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGACAATTCAGa                                              >  2:4831113/1‑150 (MQ=255)
                                                                                                tttCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGa                                              <  1:1014887/150‑1 (MQ=255)
                                                                                                tttCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGa                                              >  1:1739462/1‑150 (MQ=255)
                                                                                                 ttCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGac                                             <  1:4271465/150‑1 (MQ=255)
                                                                                                 ttCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGac                                             >  1:1393222/1‑150 (MQ=255)
                                                                                                 ttCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAAATCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGac                                             <  1:4808462/150‑1 (MQ=255)
                                                                                                   cGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGacac                                           >  1:1377350/1‑150 (MQ=255)
                                                                                                    gTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                          >  2:377279/1‑150 (MQ=255)
                                                                                                    gTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                          >  2:4562371/1‑150 (MQ=255)
                                                                                                     tCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCa                                         <  1:4333489/150‑1 (MQ=255)
                                                                                                     tCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCa                                         <  1:3629002/150‑1 (MQ=255)
                                                                                                      cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAg                                        <  1:1717712/150‑1 (MQ=255)
                                                                                                      cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAg                                        >  1:3732642/1‑150 (MQ=255)
                                                                                                      cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAg                                        <  2:2661404/150‑1 (MQ=255)
                                                                                                          gAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                    >  2:32942/1‑150 (MQ=255)
                                                                                                            ggTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAtt                                  <  1:444097/150‑1 (MQ=255)
                                                                                                                gtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGt                              >  1:3926603/1‑150 (MQ=255)
                                                                                                                  aCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGttt                            >  2:303729/1‑150 (MQ=255)
                                                                                                                   cGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGtttt                           <  2:1595664/150‑1 (MQ=255)
                                                                                                                    gTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGttttt                          <  2:901967/150‑1 (MQ=255)
                                                                                                                     ttGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGtttttt                         <  1:288260/150‑1 (MQ=255)
                                                                                                                     ttGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGtttttt                         >  1:1028218/1‑150 (MQ=255)
                                                                                                                      tGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTa                        >  2:1599629/1‑150 (MQ=255)
                                                                                                                       ggCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTaa                       >  1:3420550/1‑150 (MQ=255)
                                                                                                                          gTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATg                    >  2:2731310/1‑150 (MQ=255)
                                                                                                                                gCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAAc              >  1:2392535/1‑150 (MQ=255)
                                                                                                                                  cAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCg            <  2:1014536/150‑1 (MQ=255)
                                                                                                                                   aGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCgg           <  1:363990/150‑1 (MQ=255)
                                                                                                                                      aTAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCtt        >  2:541176/1‑150 (MQ=255)
                                                                                                                                         aTTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtat     <  1:1974689/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACA  >  CP009273/2122310‑2122599

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: