Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,864,463 C→A 35.6% D210Y (GAT→TAT)  dgoT ← D‑galactonate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,864,4630CA35.6% 377.9 / 229.6 241D210Y (GAT→TAT) dgoTD‑galactonate transporter
Reads supporting (aligned to +/- strand):  ref base C (69/86);  new base A (41/44);  total (110/131)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 241 total aligned reads displayed.

AAAACCACAGTGTAGAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGAATCCAGA  >  CP009273/3864316‑3864609
                                                                                                                                                   |                                                                                                                                                  
     cacaGTGTAGAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATg                                                                                                                                             >  2:953294/1‑150 (MQ=255)
     cacaGTGTAGAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATg                                                                                                                                             >  1:2652657/1‑150 (MQ=255)
        aGTGTAGAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTaa                                                                                                                                          >  1:3794449/1‑150 (MQ=255)
        aGTGTAGAAGCCACCGCAAATTGCCCAAGAGAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTaa                                                                                                                                          >  1:2379332/1‑150 (MQ=255)
         gtgtAGAAGCCACCGCAAATTGCACAAGTTAAACGCAGATCAGTTTTCGATGGAACACCAGTTTCCAGTCTTTGGCTGTTATCGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATTACGAATTTAAt                                                                                                                                         <  2:1816934/150‑1 (MQ=255)
          tgtAGAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATc                                                                                                                                        >  2:4332028/1‑150 (MQ=255)
            tAGAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCa                                                                                                                                      <  1:1073744/150‑1 (MQ=255)
             aGAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAg                                                                                                                                     >  1:4258856/1‑150 (MQ=255)
              gAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGt                                                                                                                                    >  2:642121/1‑150 (MQ=255)
              gAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGt                                                                                                                                    <  2:2372329/150‑1 (MQ=255)
              gAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGt                                                                                                                                    >  2:477582/1‑150 (MQ=255)
               aaGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGtt                                                                                                                                   >  1:374133/1‑150 (MQ=255)
               aaGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGtt                                                                                                                                   >  2:31014/1‑150 (MQ=255)
               aaGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGtt                                                                                                                                   >  1:2041461/1‑150 (MQ=255)
                aGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTc                                                                                                                                  <  2:3384538/150‑1 (MQ=255)
                 gCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCa                                                                                                                                 <  2:4551873/150‑1 (MQ=255)
                        cAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGc                                                                                                                          >  2:2700697/1‑150 (MQ=255)
                         aaaTTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCt                                                                                                                         >  1:2544467/1‑150 (MQ=255)
                            ttGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGAt                                                                                                                      <  2:184491/150‑1 (MQ=255)
                               cccAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATAcc                                                                                                                   >  1:607708/1‑150 (MQ=255)
                                ccAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCt                                                                                                                  <  1:449877/150‑1 (MQ=255)
                                 cAAGATAAACGCCGATCCGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCtt                                                                                                                 >  2:3963425/1‑150 (MQ=255)
                                  aaGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCGTAACGAATGTAATCCAGTTCAGCTTTGCTGATACCttt                                                                                                                <  2:1361338/150‑1 (MQ=255)
                                  aaGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCttt                                                                                                                >  2:1943183/1‑150 (MQ=255)
                                  aaGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCttt                                                                                                                >  1:3057406/1‑150 (MQ=255)
                                   aGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTg                                                                                                               <  1:3287907/150‑1 (MQ=255)
                                   aGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTg                                                                                                               >  2:1525341/1‑150 (MQ=255)
                                   aGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTg                                                                                                               >  2:4168283/1‑150 (MQ=255)
                                   aGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTg                                                                                                               >  1:756099/1‑150 (MQ=255)
                                    gATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTgg                                                                                                              <  2:1344295/150‑1 (MQ=255)
                                      tAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTc                                                                                                            <  1:1093233/150‑1 (MQ=255)
                                       aaaCGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCa                                                                                                           >  1:1665434/1‑150 (MQ=255)
                                           gCCGATCAGTTTACGATGGAACACCATTTTCCAGTCTTTGGATGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACTAATGTAATCCATTTCAGCTTTGCTGATACCTTTGGTCAGgcg                                                                                                       <  1:1723234/150‑1 (MQ=255)
                                           gCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGgcg                                                                                                       <  1:4218141/150‑1 (MQ=255)
                                            ccGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGgcgc                                                                                                      >  2:2791202/1‑150 (MQ=255)
                                            ccGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAAAGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGgcgc                                                                                                      <  1:2624689/150‑1 (MQ=255)
                                                tCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCgcggcg                                                                                                  <  2:1237681/150‑1 (MQ=255)
                                                 cAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGcggcgg                                                                                                 >  2:1184823/1‑150 (MQ=255)
                                                     ttACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGa                                                                                             >  1:3858354/1‑150 (MQ=255)
                                                      tACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGAt                                                                                            <  2:4512507/150‑1 (MQ=255)
                                                      tACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGAt                                                                                            <  2:3479047/150‑1 (MQ=255)
                                                       aCGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATa                                                                                           <  1:4417232/150‑1 (MQ=255)
                                                        cGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATaa                                                                                          <  1:901389/150‑1 (MQ=255)
                                                        cGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATaa                                                                                          <  1:3837184/150‑1 (MQ=255)
                                                           tGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAAcc                                                                                       <  2:3243523/150‑1 (MQ=255)
                                                            ggAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCt                                                                                      <  1:4087686/150‑1 (MQ=255)
                                                             gAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCtt                                                                                     <  1:1634000/150‑1 (MQ=255)
                                                              aaCACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTa                                                                                    >  2:2008713/1‑150 (MQ=255)
                                                                cacCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTaaa                                                                                  <  1:3564426/150‑1 (MQ=255)
                                                                 acCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAAc                                                                                 <  2:4642620/150‑1 (MQ=255)
                                                                 acCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAAc                                                                                 <  1:1984603/150‑1 (MQ=255)
                                                                   cAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCa                                                                               <  2:1699163/150‑1 (MQ=255)
                                                                   cAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCAACATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAAACa                                                                               <  1:102887/150‑1 (MQ=255)
                                                                        tCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCa                                                                          <  1:3953051/150‑1 (MQ=255)
                                                                         ccAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTTAACCAAATCAg                                                                         <  2:3463607/150‑1 (MQ=255)
                                                                         ccAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAg                                                                         <  2:1140351/150‑1 (MQ=255)
                                                                         ccAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAg                                                                         >  2:3912870/1‑150 (MQ=255)
                                                                            gTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGa                                                                      <  1:2729394/150‑1 (MQ=255)
                                                                            gTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGa                                                                      <  1:3154449/150‑1 (MQ=255)
                                                                              cTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGAcc                                                                    <  1:1324895/150‑1 (MQ=255)
                                                                                 tGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGa                                                                 <  2:3812121/150‑1 (MQ=255)
                                                                                  ggCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGAt                                                                <  1:2237129/150‑1 (MQ=255)
                                                                                   gCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATa                                                               <  1:4543218/150‑1 (MQ=255)
                                                                                    cTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATaa                                                              <  2:4136916/150‑1 (MQ=255)
                                                                                     tGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAAt                                                             <  2:2447647/150‑1 (MQ=255)
                                                                                      gTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATg                                                            <  2:20163/150‑1 (MQ=255)
                                                                                       ttAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGc                                                           >  2:431163/1‑150 (MQ=255)
                                                                                       ttAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGc                                                           <  1:2078296/150‑1 (MQ=255)
                                                                                        tAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGcc                                                          >  2:3572082/1‑150 (MQ=255)
                                                                                          aCGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTTATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGa                                                        <  1:2798596/150‑1 (MQ=255)
                                                                                            ggCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATa                                                      <  2:3573100/150‑1 (MQ=255)
                                                                                               tGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATacca                                                   >  1:3531084/1‑150 (MQ=255)
                                                                                                 aCGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATaccacc                                                 <  1:2394258/150‑1 (MQ=255)
                                                                                                  cgcgCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATAccacca                                                <  2:2719678/150‑1 (MQ=255)
                                                                                                       ctctTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGAc                                           >  2:3895170/1‑150 (MQ=255)
                                                                                                        tctttcttCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACa                                          <  1:3895117/150‑1 (MQ=255)
                                                                                                        tctttcttCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACa                                          <  2:3472065/150‑1 (MQ=255)
                                                                                                         ctttcttCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACaa                                         >  1:4230826/1‑150 (MQ=255)
                                                                                                         ctttcttCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACaa                                         >  2:286709/1‑150 (MQ=255)
                                                                                                          tttcttCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGAGGCGGCTGATATAGCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAAt                                        <  1:2208595/150‑1 (MQ=255)
                                                                                                          tttcttCACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAAt                                        <  2:2854158/150‑1 (MQ=255)
                                                                                                           ttcttcACCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATg                                       <  2:3322844/150‑1 (MQ=255)
                                                                                                            tcttcACCGGCGCATCACCAGCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGa                                      <  2:588539/150‑1 (MQ=255)
                                                                                                                 aCCGGCGCATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACAcc                                 >  1:3076300/1‑150 (MQ=255)
                                                                                                                       gcATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGc                           >  1:1441650/1‑150 (MQ=255)
                                                                                                                       gcATCACCATCCACCAGACCGCCGCCATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGc                           <  2:3057406/150‑1 (MQ=255)
                                                                                                                                tCCACCAGACCGCCGCCATCACTAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCTACCAGATAATGCCGATACCACCAGTGACAATGAACACACAGTGCCAGCTCAAc                  <  2:3730338/150‑1 (MQ=255)
                                                                                                                                tCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAAc                  <  2:1840902/150‑1 (MQ=255)
                                                                                                                                      aGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTcc            <  1:3963361/150‑1 (MQ=255)
                                                                                                                                       gACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCt           >  2:2542319/1‑150 (MQ=255)
                                                                                                                                       gACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCt           <  2:4059057/150‑1 (MQ=255)
                                                                                                                                          cgccgcCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGaa        >  2:390676/1‑150 (MQ=255)
                                                                                                                                          cgccgcCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGaa        >  2:427212/1‑150 (MQ=255)
                                                                                                                                             cgccATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGAATcc     <  1:4168278/150‑1 (MQ=255)
                                                                                                                                              gccATCTCTATTGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCTGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCCGTGCCAGCTCAACATCTCCTGAATCCa    <  1:473400/150‑1 (MQ=255)
                                                                                                                                               ccATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGAATCCAg   >  1:392201/1‑150 (MQ=255)
                                                                                                                                               ccATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGAATCCAg   >  1:3738985/1‑150 (MQ=255)
                                                                                                                                                cATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGAATCCAGa  >  2:221437/1‑150 (MQ=255)
                                                                                                                                                cATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGAATCCAGa  >  1:3768154/1‑150 (MQ=255)
                                                                                                                                                cATAACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGAATCCAGa  <  1:3364010/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAACCACAGTGTAGAAGCCACCGCAAATTGCCCAAGATAAACGCCGATCAGTTTACGATGGAACACCAGTTTCCAGTCTTTGGCTGTTAACGGCTGACGCGCCTCTTTCTTCACCGGCGCATCACCATCCACCAGACCGCCGCCATCACGAATGTAATCCAGTTCAGCTTTGCTGATACCTTTGGTCAGGCGCGGCGGCTGATAAACCTTAAACCAAATCAGCGACCAGATAATGCCGATACCACCAGTGACAATGAACACCCAGTGCCAGCTCAACATCTCCTGAATCCAGA  >  CP009273/3864316‑3864609

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: