Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,826,454 C→A 5.3% D482Y (GAT→TAT)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,826,4540CA 5.3% 861.9 / 14.8 224D482Y (GAT→TAT) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (90/122);  new base A (7/5);  total (97/127)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 228 total aligned reads displayed.

CAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCATTTTCTGCGG  >  CP009273/3826307‑3826601
                                                                                                                                                   |                                                                                                                                                   
cacacGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATc                                                                                                                                                     <  2:3485048/145‑1 (MQ=255)
 aaacGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCg                                                                                                                                                    >  2:1085121/1‑148 (MQ=255)
 aaacGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCg                                                                                                                                                    <  1:2063519/148‑1 (MQ=255)
   acGCGCACCAGCGTTTGTACACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGACAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACACTCGcc                                                                                                                                                  <  2:2801104/148‑1 (MQ=255)
   acGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGcc                                                                                                                                                  >  1:2354132/1‑148 (MQ=255)
   acGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGcc                                                                                                                                                  <  2:866554/148‑1 (MQ=255)
    cgcgCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCTCCATTGATTCGTAGTTAGCGTAACAATCGCCa                                                                                                                                                 <  2:3145333/148‑1 (MQ=255)
       gcACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCAccc                                                                                                                                              <  1:284439/148‑1 (MQ=255)
        cACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCAcccc                                                                                                                                             <  2:2290581/148‑1 (MQ=255)
            aGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAgtg                                                                                                                                         <  2:1200054/148‑1 (MQ=255)
             gCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAgtgt                                                                                                                                        <  1:3439923/148‑1 (MQ=255)
             gCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAgtgt                                                                                                                                        <  1:3363761/148‑1 (MQ=255)
              cGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTa                                                                                                                                       <  1:3127642/148‑1 (MQ=255)
               gTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTAc                                                                                                                                      <  1:1770562/148‑1 (MQ=255)
               gTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTAc                                                                                                                                      <  2:3533521/148‑1 (MQ=255)
                 ttGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCg                                                                                                                                    >  1:1354657/1‑148 (MQ=255)
                    tAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGa                                                                                                                                 >  2:2256656/1‑148 (MQ=255)
                     aaaCGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGaa                                                                                                                                <  1:1491826/148‑1 (MQ=255)
                     aaaCGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGaa                                                                                                                                >  1:817418/1‑148 (MQ=255)
                       aCGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAAtt                                                                                                                              >  2:3010046/1‑148 (MQ=255)
                              ccGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTgcgc                                                                                                                       >  2:2858850/1‑148 (MQ=255)
                                ggAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCACTAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCAc                                                                                                                     <  2:2144412/148‑1 (MQ=255)
                                 gAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCCCCCCAGTGTACCGGGAATTTCTGCGCAcc                                                                                                                    <  1:425751/148‑1 (MQ=255)
                                  aGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCg                                                                                                                   <  1:2319937/148‑1 (MQ=255)
                                       tATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGa                                                                                                              >  1:874036/1‑148 (MQ=255)
                                       tATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGa                                                                                                              >  1:3015674/1‑148 (MQ=255)
                                       tATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATAGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGa                                                                                                              <  1:2254800/148‑1 (MQ=255)
                                       tATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATAGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGa                                                                                                              >  2:1673083/1‑148 (MQ=255)
                                         ttttCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAgg                                                                                                            <  2:3554691/148‑1 (MQ=255)
                                           ttCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGcc                                                                                                          <  1:2199626/148‑1 (MQ=255)
                                            tCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCg                                                                                                         <  2:3276675/148‑1 (MQ=255)
                                            tCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCg                                                                                                         <  2:1959483/148‑1 (MQ=255)
                                              aaaGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAg                                                                                                       <  1:1609043/148‑1 (MQ=255)
                                                aGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCg                                                                                                     >  1:2288315/1‑148 (MQ=255)
                                                aGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCg                                                                                                     <  2:2631818/148‑1 (MQ=255)
                                                  ccgccgATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCggg                                                                                                   >  2:1317523/1‑148 (MQ=255)
                                                     ccgATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAACCCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCaa                                                                                                <  2:687944/148‑1 (MQ=255)
                                                        atatCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACa                                                                                             <  2:813639/148‑1 (MQ=255)
                                                        atatCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACa                                                                                             <  2:813907/148‑1 (MQ=255)
                                                         tatCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACaa                                                                                            >  1:1647711/1‑148 (MQ=255)
                                                          atCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAg                                                                                           <  1:3513673/148‑1 (MQ=255)
                                                           tCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGa                                                                                          <  1:1229103/148‑1 (MQ=255)
                                                            cGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGAc                                                                                         <  2:946208/148‑1 (MQ=255)
                                                                gCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGct                                                                                     >  2:1531939/1‑148 (MQ=255)
                                                                    cAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGcttctt                                                                                 >  1:3304139/1‑148 (MQ=255)
                                                                      gAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTcc                                                                               >  1:2268324/1‑148 (MQ=255)
                                                                       aaGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCt                                                                              >  2:2159656/1‑148 (MQ=255)
                                                                         gCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCa                                                                            >  1:774255/1‑148 (MQ=255)
                                                                             aaaCCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCaa                                                                        <  1:1748104/148‑1 (MQ=255)
                                                                             aaaCCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATAGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCaa                                                                        >  2:447592/1‑148 (MQ=255)
                                                                             aaaCCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATAGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCaa                                                                        >  2:2237486/1‑148 (MQ=255)
                                                                               aCCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACg                                                                      <  1:3473112/148‑1 (MQ=255)
                                                                               aCCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATAGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACg                                                                      >  1:2938796/1‑148 (MQ=255)
                                                                                ccTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACgg                                                                     >  2:711227/1‑148 (MQ=255)
                                                                                ccTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACgg                                                                     <  1:242958/148‑1 (MQ=255)
                                                                                 cTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGgt                                                                    <  1:2377690/148‑1 (MQ=255)
                                                                                  tGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGgtg                                                                   <  2:954792/148‑1 (MQ=255)
                                                                                  tGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGgtg                                                                   <  2:954642/148‑1 (MQ=255)
                                                                                   gAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGgtgt                                                                  <  1:2367608/148‑1 (MQ=255)
                                                                                       ggCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCtt                                                              >  2:2054831/1‑148 (MQ=255)
                                                                                        gCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTg                                                             <  2:3204656/148‑1 (MQ=255)
                                                                                          cAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAg                                                           >  1:684628/1‑148 (MQ=255)
                                                                                           aaTAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGc                                                          >  1:3408859/1‑148 (MQ=255)
                                                                                           aaTAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGc                                                          >  1:324163/1‑148 (MQ=255)
                                                                                           aaTAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGc                                                          >  1:323743/1‑148 (MQ=255)
                                                                                           aaTAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGc                                                          >  1:3051825/1‑148 (MQ=255)
                                                                                             tAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCAc                                                        >  2:1935366/1‑148 (MQ=255)
                                                                                                acacacCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGtt                                                     <  2:397595/148‑1 (MQ=255)
                                                                                                 cAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTc                                                    >  2:1164507/1‑148 (MQ=255)
                                                                                                    accaccGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCcac                                                 >  2:1309694/1‑148 (MQ=255)
                                                                                                    accaccGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCcac                                                 >  2:2727747/1‑148 (MQ=255)
                                                                                                    accaccGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCcac                                                 <  1:1380123/148‑1 (MQ=255)
                                                                                                     ccaccGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCcaca                                                >  1:2620880/1‑148 (MQ=255)
                                                                                                      caccGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCacac                                               <  2:1910294/148‑1 (MQ=255)
                                                                                                        ccGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCa                                             <  2:656029/148‑1 (MQ=255)
                                                                                                          gcgcAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGt                                           <  2:1363747/148‑1 (MQ=255)
                                                                                                           cgcAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGtt                                          <  2:1591498/148‑1 (MQ=255)
                                                                                                               ggCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGtt                                      >  1:2247784/1‑148 (MQ=255)
                                                                                                                gCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTg                                     >  1:533515/1‑148 (MQ=255)
                                                                                                                gCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTg                                     >  1:3346440/1‑148 (MQ=255)
                                                                                                                gCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTg                                     <  2:2336071/148‑1 (MQ=255)
                                                                                                                   ttCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAg                                  <  2:3157707/148‑1 (MQ=255)
                                                                                                                   ttCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAg                                  >  1:400435/1‑148 (MQ=255)
                                                                                                                       gCCATTGATTCGTAGTTAGCGTAACAATAGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGt                              <  1:1539956/148‑1 (MQ=255)
                                                                                                                        ccATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTa                             <  1:39256/148‑1 (MQ=255)
                                                                                                                           ttGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGc                          <  1:278570/148‑1 (MQ=255)
                                                                                                                            tGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCa                         >  1:599215/1‑148 (MQ=255)
                                                                                                                                   tAGTTAGCGTAACACTCGCCACCCCAGTGTACCGGGAATGTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGt                  <  1:1818669/148‑1 (MQ=255)
                                                                                                                                   tAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGt                  <  1:1317523/148‑1 (MQ=255)
                                                                                                                                    aGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGtt                 <  2:306749/148‑1 (MQ=255)
                                                                                                                                     gTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTa                <  2:3304739/148‑1 (MQ=255)
                                                                                                                                     gTTAGCGTAACAATAGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTa                <  1:2698834/148‑1 (MQ=255)
                                                                                                                                       tAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATg              <  1:1085121/148‑1 (MQ=255)
                                                                                                                                         gCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCa            <  2:1354707/148‑1 (MQ=255)
                                                                                                                                          cGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCAt           <  1:1953334/148‑1 (MQ=255)
                                                                                                                                             aaCAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCAtttt        >  1:3446973/1‑148 (MQ=255)
                                                                                                                                             aaCAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCAtttt        >  1:3337380/1‑148 (MQ=255)
                                                                                                                                                 aTCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTCGATGTACGCATAATGGTTATGCATTTTCTGc    >  2:2412366/1‑148 (MQ=255)
                                                                                                                                                   cGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCATTTTCTGCgg  >  1:915248/1‑148 (MQ=255)
                                                                                                                                                   cGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCATTTTCTGCgg  >  1:2020899/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAAACGCGCACCAGCGTTTGTAAACGTGCGCCGGAGCGGTATTTTCAAAGCCGCCGATATCGTGGCTCCAGAAGCCAAAACCTGAAAGGCCAATAGACAAACCACCGCGCAGGCTTTCCGCCATTGATTCGTAGTTAGCGTAACAATCGCCACCCCAGTGTACCGGGAATTTCTGCGCACCGACGGAGGCCGAGCGGGCAAACAAGACAGCTTCTTCCTCACCAACGGTGTCCTTGAGCACGTTCCACACCAGTTCGTTGTAGATGTACGCATAATGGTTATGCATTTTCTGCGG  >  CP009273/3826307‑3826601

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: ‑