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,4530CT98.5% 861.1 / 2.3 262A227T (GCT→ACT) wcaEputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (123/135);  total (126/136)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-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 262 total aligned reads displayed.

TTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTAT  >  CP009273/2122310‑2122596
                                                                                                                                               |                                                                                                                                               
   aaTAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAg                                                                                                                                          <  1:1974209/148‑1 (MQ=255)
         gTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGc                                                                                                                                    <  1:3431612/148‑1 (MQ=255)
           tGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGcac                                                                                                                                  >  1:212575/1‑148 (MQ=255)
           tGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGcac                                                                                                                                  >  1:2410450/1‑148 (MQ=255)
             ccgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACAt                                                                                                                                <  2:2845184/148‑1 (MQ=255)
             ccgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACAt                                                                                                                                <  1:3475278/148‑1 (MQ=255)
             ccgccCCGGAACCCTTTCGGCACCGAGAAACAGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACAt                                                                                                                                >  2:1548927/1‑148 (MQ=255)
                ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTa                                                                                                                             >  2:1925255/1‑148 (MQ=255)
                ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTa                                                                                                                             <  2:1771864/148‑1 (MQ=255)
                  ccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAat                                                                                                                           <  1:485850/148‑1 (MQ=255)
                  ccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAat                                                                                                                           <  1:2161086/148‑1 (MQ=255)
                     gAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATAtt                                                                                                                        >  2:1482778/1‑148 (MQ=255)
                     gAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATAtt                                                                                                                        >  2:1001778/1‑148 (MQ=255)
                     gAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATAtt                                                                                                                        >  2:2757414/1‑148 (MQ=255)
                      aaCCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATAttt                                                                                                                       >  2:1213629/1‑148 (MQ=255)
                      aaCCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATAttt                                                                                                                       >  1:928713/1‑148 (MQ=255)
                        cccTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGt                                                                                                                     <  2:1748419/148‑1 (MQ=255)
                         ccTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTc                                                                                                                    <  1:1148198/148‑1 (MQ=255)
                          cTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCg                                                                                                                   <  1:3508056/148‑1 (MQ=255)
                             tCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTg                                                                                                                <  1:2093188/148‑1 (MQ=255)
                               ggCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGa                                                                                                              <  2:1065643/148‑1 (MQ=255)
                                 cACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACt                                                                                                            <  1:430131/148‑1 (MQ=255)
                                   ccGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACttt                                                                                                          <  2:849180/148‑1 (MQ=255)
                                     gagaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACttttt                                                                                                        <  1:536612/148‑1 (MQ=255)
                                      agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACtttttt                                                                                                       <  2:320643/148‑1 (MQ=255)
                                      agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACtttttt                                                                                                       >  2:26939/1‑148 (MQ=255)
                                       gaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTc                                                                                                      >  2:3405940/1‑148 (MQ=255)
                                       gaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTc                                                                                                      >  1:198458/1‑148 (MQ=255)
                                       gaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTCGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTAAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTc                                                                                                      >  1:197971/1‑148 (MQ=255)
                                        aaaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcg                                                                                                     <  2:2096372/148‑1 (MQ=255)
                                        aaaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcg                                                                                                     >  1:1247281/1‑148 (MQ=255)
                                         aaCCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTcgc                                                                                                    >  2:1661444/1‑148 (MQ=255)
                                           ccGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGt                                                                                                  >  1:871424/1‑148 (MQ=255)
                                                tAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGca                                                                                             <  2:2196653/148‑1 (MQ=255)
                                                 aaaTCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGcac                                                                                            <  1:1746398/148‑1 (MQ=255)
                                                     cTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCATTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACaa                                                                                        <  1:737369/148‑1 (MQ=255)
                                                      ttGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAAt                                                                                       >  1:1619887/1‑148 (MQ=255)
                                                      ttGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAAt                                                                                       <  2:2500746/148‑1 (MQ=255)
                                                      ttGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAAt                                                                                       <  2:1598238/148‑1 (MQ=255)
                                                       tGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAAtt                                                                                      >  1:2020874/1‑148 (MQ=255)
                                                         cATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTcc                                                                                    >  2:257241/1‑148 (MQ=255)
                                                          aTTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCa                                                                                   <  1:2398942/148‑1 (MQ=255)
                                                           ttGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCat                                                                                  >  1:1859648/1‑148 (MQ=255)
                                                           ttGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCat                                                                                  >  2:1646068/1‑148 (MQ=255)
                                                           ttGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACATTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCat                                                                                  >  1:1859653/1‑148 (MQ=255)
                                                            tGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCata                                                                                 >  1:373207/1‑148 (MQ=255)
                                                            tGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCata                                                                                 <  2:1322801/148‑1 (MQ=255)
                                                               ttttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCCtatta                                                                              >  2:3320937/1‑148 (MQ=255)
                                                                tttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACATAATTCCAtattat                                                                             >  2:1671718/1‑148 (MQ=255)
                                                                 ttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATattatt                                                                            <  2:2904397/148‑1 (MQ=255)
                                                                  tCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtg                                                                           <  1:2326417/148‑1 (MQ=255)
                                                                    cTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACCATTCCATATTATtggt                                                                         >  2:2273551/1‑148 (MQ=255)
                                                                    cTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggt                                                                         <  2:1094013/148‑1 (MQ=255)
                                                                     ttATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggtg                                                                        <  2:3052201/148‑1 (MQ=255)
                                                                         atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAg                                                                    >  2:2877906/1‑148 (MQ=255)
                                                                         atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAg                                                                    <  1:2648427/148‑1 (MQ=255)
                                                                          ttCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGa                                                                   <  1:2642617/148‑1 (MQ=255)
                                                                          ttCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGa                                                                   >  2:1579426/1‑148 (MQ=255)
                                                                           tCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGAt                                                                  <  1:126192/148‑1 (MQ=255)
                                                                            cAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATa                                                                 >  1:568397/1‑148 (MQ=255)
                                                                            cAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATa                                                                 <  1:1925255/148‑1 (MQ=255)
                                                                            cAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATa                                                                 <  2:3368467/148‑1 (MQ=255)
                                                                               aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAccc                                                              >  1:1344600/1‑148 (MQ=255)
                                                                               aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAccc                                                              <  1:157137/148‑1 (MQ=255)
                                                                               aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAccc                                                              <  1:789199/148‑1 (MQ=255)
                                                                                  ttGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACcccac                                                           >  2:2189063/1‑148 (MQ=255)
                                                                                     ttATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACcccaccca                                                        <  1:3189687/148‑1 (MQ=255)
                                                                                      tatataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCAt                                                       <  1:1991525/148‑1 (MQ=255)
                                                                                        tataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATgg                                                     >  2:1150352/1‑148 (MQ=255)
                                                                                          taAGGCTTTCGTCTTTGAGGTAGTACGTTGCCGTAAATCCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATATGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaa                                                   >  1:1825488/1‑148 (MQ=255)
                                                                                           aaGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaa                                                  >  2:1634026/1‑148 (MQ=255)
                                                                                           aaGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaa                                                  >  2:1393908/1‑148 (MQ=255)
                                                                                             ggCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAAt                                                >  1:1519772/1‑148 (MQ=255)
                                                                                              gCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAAtt                                               <  1:3136962/148‑1 (MQ=255)
                                                                                              gCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAAtt                                               >  2:2620009/1‑148 (MQ=255)
                                                                                                   cGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGac                                          >  2:69665/1‑148 (MQ=255)
                                                                                                     tCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGacac                                        >  1:1025283/1‑148 (MQ=255)
                                                                                                     tCTTTGAGGTAGTACGGTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGAACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGacac                                        <  2:1581943/148‑1 (MQ=255)
                                                                                                      cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                       >  1:2824091/1‑148 (MQ=255)
                                                                                                            ggTAGTACGTTTGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                 <  1:2148937/148‑1 (MQ=255)
                                                                                                            ggTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                 <  2:621532/148‑1 (MQ=255)
                                                                                                            ggTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCa                                 <  2:488563/148‑1 (MQ=255)
                                                                                                             gtagtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAt                                <  2:3515799/148‑1 (MQ=255)
                                                                                                               agtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTg                              <  1:2875038/148‑1 (MQ=255)
                                                                                                                gtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGa                             <  2:441981/148‑1 (MQ=255)
                                                                                                                 taCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAg                            >  1:701442/1‑148 (MQ=255)
                                                                                                                 taCGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAg                            >  1:2097777/1‑148 (MQ=255)
                                                                                                                  aCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGt                           <  2:1770476/148‑1 (MQ=255)
                                                                                                                      tGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGttttt                       <  1:3263271/148‑1 (MQ=255)
                                                                                                                      tGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGttttt                       <  1:2052461/148‑1 (MQ=255)
                                                                                                                        gCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTa                     >  1:727063/1‑148 (MQ=255)
                                                                                                                           tAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAAt                  >  2:914277/1‑148 (MQ=255)
                                                                                                                              aTGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCa               <  2:2095181/148‑1 (MQ=255)
                                                                                                                                gCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCGATTGAGTTTTTTAAATGCATa             <  2:2402228/148‑1 (MQ=255)
                                                                                                                                gCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATa             <  1:2195214/148‑1 (MQ=255)
                                                                                                                                      aTAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAGTGAGTTTTTTAAATGCATAACCGGc       >  1:2717587/1‑148 (MQ=255)
                                                                                                                                        aaTTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCtt     >  2:639129/1‑148 (MQ=255)
                                                                                                                                          ttCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTta   >  2:823567/1‑148 (MQ=255)
                                                                                                                                          ttCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtat  >  2:3098857/1‑148 (MQ=255)
                                                                                                                                           tCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtat  >  2:146/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: