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.5% 781.3 / 6.9 240A227T (GCT→ACT) wcaEputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (115/119);  total (119/121)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 240 total aligned reads displayed.

TTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATA  >  CP009273/2122311‑2122597
                                                                                                                                              |                                                                                                                                                
ttAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTGTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCaaa                                                                                                                                             >  2:12707/1‑148 (MQ=255)
 tAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCaaaa                                                                                                                                            <  1:23242/148‑1 (MQ=255)
  aaTAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAg                                                                                                                                           <  1:333411/148‑1 (MQ=255)
  aaTAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAg                                                                                                                                           <  1:305797/148‑1 (MQ=255)
   aTAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGc                                                                                                                                          >  2:3103947/1‑148 (MQ=255)
    tAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGcc                                                                                                                                         >  2:350604/1‑148 (MQ=255)
       cGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAgg                                                                                                                                      >  1:1107038/1‑148 (MQ=255)
       cGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAgg                                                                                                                                      >  1:2525384/1‑148 (MQ=255)
             cgccCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATg                                                                                                                                <  1:1023865/148‑1 (MQ=255)
               ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTa                                                                                                                              >  1:1898072/1‑148 (MQ=255)
               ccccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTa                                                                                                                              >  1:3450456/1‑148 (MQ=255)
                cccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTaa                                                                                                                             <  1:927944/148‑1 (MQ=255)
                 ccGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAat                                                                                                                            >  1:1724486/1‑148 (MQ=255)
                     aaCCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATAttt                                                                                                                        <  1:2171749/148‑1 (MQ=255)
                      aCCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTCATATTTg                                                                                                                       >  2:3418393/1‑148 (MQ=255)
                      aCCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTg                                                                                                                       >  2:81755/1‑148 (MQ=255)
                        ccTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTc                                                                                                                     >  1:807328/1‑148 (MQ=255)
                         cTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCg                                                                                                                    <  2:2156319/148‑1 (MQ=255)
                                cACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACt                                                                                                             <  2:809826/148‑1 (MQ=255)
                                    gagaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACttttt                                                                                                         <  2:323389/148‑1 (MQ=255)
                                     agaAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACtttttt                                                                                                        <  2:3057627/148‑1 (MQ=255)
                                            gCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCa                                                                                                 >  1:1273030/1‑148 (MQ=255)
                                             cTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAg                                                                                                >  1:877069/1‑148 (MQ=255)
                                              ttAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGc                                                                                               >  1:2501376/1‑148 (MQ=255)
                                              ttAAATCGTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGc                                                                                               <  1:945442/148‑1 (MQ=255)
                                               tAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGca                                                                                              >  1:2219444/1‑148 (MQ=255)
                                                aaaTCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGcac                                                                                             >  1:1974166/1‑148 (MQ=255)
                                                  aTCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGcacac                                                                                           <  1:1631158/148‑1 (MQ=255)
                                                   tCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGcacaca                                                                                          <  2:1484426/148‑1 (MQ=255)
                                                    cTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACaa                                                                                         <  2:2511531/148‑1 (MQ=255)
                                                    cTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACaa                                                                                         >  1:1716887/1‑148 (MQ=255)
                                                    cTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACaa                                                                                         <  1:2298740/148‑1 (MQ=255)
                                                    cTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACaa                                                                                         <  1:2106863/148‑1 (MQ=255)
                                                        cATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTcc                                                                                     <  2:1793696/148‑1 (MQ=255)
                                                        cATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTcc                                                                                     <  1:3025422/148‑1 (MQ=255)
                                                        cATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTcc                                                                                     <  1:3266748/148‑1 (MQ=255)
                                                           tGGTTTTCCTTATATTCAGACTTTTTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAAGTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCata                                                                                  <  2:3187841/148‑1 (MQ=255)
                                                           tGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCata                                                                                  <  2:178366/148‑1 (MQ=255)
                                                           tGGTTTTCCTTATATTCAGAATTTGTTATATAAGGCTTTCGTCTTTGAGGTATGACGTTGGCGTAACTGCCAGGATAATTCAGTCAAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCata                                                                                  <  1:2611871/148‑1 (MQ=255)
                                                              ttttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCAtatta                                                                               <  2:1447226/148‑1 (MQ=255)
                                                               tttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCAtattat                                                                              >  1:2063444/1‑148 (MQ=255)
                                                                ttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATattatt                                                                             <  2:656890/148‑1 (MQ=255)
                                                                ttCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATattatt                                                                             >  1:2439071/1‑148 (MQ=255)
                                                                 tCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtg                                                                            <  1:1462896/148‑1 (MQ=255)
                                                                  ccTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtgg                                                                           <  1:1922948/148‑1 (MQ=255)
                                                                    ttATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggtg                                                                         <  1:2883178/148‑1 (MQ=255)
                                                                     tataTTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggtgg                                                                        >  2:868956/1‑148 (MQ=255)
                                                                     tataTTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATtggtgg                                                                        >  2:2655677/1‑148 (MQ=255)
                                                                        atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAg                                                                     >  1:2672810/1‑148 (MQ=255)
                                                                        atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAg                                                                     >  1:524015/1‑148 (MQ=255)
                                                                        atTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAg                                                                     <  1:122308/148‑1 (MQ=255)
                                                                         ttCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGa                                                                    <  2:1995926/148‑1 (MQ=255)
                                                                          tCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGAt                                                                   >  2:536207/1‑148 (MQ=255)
                                                                             gACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAcc                                                                <  2:1107038/148‑1 (MQ=255)
                                                                             gACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAcc                                                                <  2:2885607/148‑1 (MQ=255)
                                                                              aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTTGTGGTAGATAccc                                                               >  2:2395163/1‑148 (MQ=255)
                                                                              aCTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAccc                                                               >  2:2698480/1‑148 (MQ=255)
                                                                               cTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATAcccc                                                              >  2:1166299/1‑148 (MQ=255)
                                                                                    ttATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACcccaccca                                                         >  2:923125/1‑148 (MQ=255)
                                                                                      atataAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATg                                                       <  2:1632162/148‑1 (MQ=255)
                                                                                         taAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaa                                                    >  1:2241536/1‑148 (MQ=255)
                                                                                           aGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGaaaa                                                  >  1:380816/1‑148 (MQ=255)
                                                                                            ggCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAAt                                                 >  2:3548262/1‑148 (MQ=255)
                                                                                             gCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAAtt                                                >  2:1586013/1‑148 (MQ=255)
                                                                                                  cGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGac                                           >  1:2925073/1‑148 (MQ=255)
                                                                                                    tCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGacac                                         <  2:1664659/148‑1 (MQ=255)
                                                                                                    tCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGacac                                         >  1:449842/1‑148 (MQ=255)
                                                                                                    tCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACAAAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAATTTCAGactc                                         <  2:1460867/148‑3 (MQ=255)
                                                                                                     cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                        <  2:625565/148‑1 (MQ=255)
                                                                                                     cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                        <  1:399947/148‑1 (MQ=255)
                                                                                                     cTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACAcc                                        >  2:3408372/1‑148 (MQ=255)
                                                                                                      tttGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCa                                       <  2:2267917/148‑1 (MQ=255)
                                                                                                          agatagaaCGTTGGCGTAAATGACAGGATAATTCAGTCCAAAAGACAGGAAAATGTAAAATTTGTCTTAGTACTTTTTTCGATTCACCACACAAATCCATATTATTGGTGGTAGATACAACACACAAGGAAAATACATAAACAAGGcc                                   >  2:1504952/4‑148 (MQ=255)
                                                                                                            gtagtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAt                                 >  2:910580/1‑148 (MQ=255)
                                                                                                            gtagtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAt                                 >  2:72680/1‑148 (MQ=255)
                                                                                                            gtagtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAt                                 >  2:1008689/1‑148 (MQ=255)
                                                                                                            gtagtaCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAt                                 >  2:1743587/1‑148 (MQ=255)
                                                                                                                 aCGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGt                            <  1:3103947/148‑1 (MQ=255)
                                                                                                                  cGTTGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGtt                           >  1:765160/1‑148 (MQ=255)
                                                                                                                    ttGGCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGtttt                         <  2:401907/148‑1 (MQ=255)
                                                                                                                      ggCGTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGtttttt                       >  1:2835006/1‑148 (MQ=255)
                                                                                                                         gTAAATGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTaaa                    <  2:2126256/148‑1 (MQ=255)
                                                                                                                             aTGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCa                >  2:606081/1‑148 (MQ=255)
                                                                                                                              tGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCAt               <  2:3533936/148‑1 (MQ=255)
                                                                                                                              tGCCAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCAt               <  2:1999846/148‑1 (MQ=255)
                                                                                                                               gCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATa              <  1:350604/148‑1 (MQ=255)
                                                                                                                                 cAGGATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATaaa            >  2:2951805/1‑147 (MQ=255)
                                                                                                                                    gATAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCgg         <  2:203560/148‑1 (MQ=255)
                                                                                                                                      tAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCt       >  1:586654/1‑148 (MQ=255)
                                                                                                                                      tAATTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCt       <  1:1820132/148‑1 (MQ=255)
                                                                                                                                        aTTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCttt     <  1:3316422/148‑1 (MQ=255)
                                                                                                                                        aTTCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCttt     <  2:449842/148‑1 (MQ=255)
                                                                                                                                         ttCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTta    <  1:2788980/148‑1 (MQ=255)
                                                                                                                                         ttCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTta    >  2:1205387/1‑148 (MQ=255)
                                                                                                                                          tCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtat   <  1:3199675/148‑1 (MQ=255)
                                                                                                                                          tCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtat   <  2:2834961/148‑1 (MQ=255)
                                                                                                                                          tCAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCAGTGAGTTTTTTAAATGCATAACCGGCTTtat   >  1:1237769/1‑148 (MQ=255)
                                                                                                                                           cAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtata  <  2:1723859/148‑1 (MQ=255)
                                                                                                                                           cAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtata  <  2:1724486/148‑1 (MQ=255)
                                                                                                                                           cAGTCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTtata  <  2:3205802/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAACCGCTTAAATCTTGCATTGGTTTTCCTTATATTCAGACTTTGTTATATAAGGCTTTCGTCTTTGAGGTAGTACGTTGGCGTAAATGCCAGGATAATTCAGCCCAAAAGCCAGGCACATGTAATATTTGTCGTTGGACTTTTTTCGCGTCAGCACACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATA  >  CP009273/2122311‑2122597

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: