Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,057,226 G→T 11.8% E839* (GAA→TAA)  torA → trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,057,2260GT11.8% 651.7 / 42.8 194E839* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (80/91);  new base T (14/9);  total (94/100)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 194 total aligned reads displayed.

CCAACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAT  >  CP009273/1057080‑1057372
                                                                                                                                                  |                                                                                                                                                  
ccAACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGa                                                                                                                                                   <  2:2789099/148‑1 (MQ=255)
  aaCGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAat                                                                                                                                                 <  1:1467995/148‑1 (MQ=255)
      tgtTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGt                                                                                                                                             <  2:368931/148‑1 (MQ=255)
       gtTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGACCAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGtt                                                                                                                                            >  2:1934178/1‑148 (MQ=255)
       gtTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGtt                                                                                                                                            >  2:2450170/1‑148 (MQ=255)
        ttGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTc                                                                                                                                           <  2:1930175/148‑1 (MQ=255)
        ttGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTc                                                                                                                                           <  2:1301994/148‑1 (MQ=255)
        ttGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTc                                                                                                                                           <  1:2089746/148‑1 (MQ=255)
             cATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCcgcg                                                                                                                                      <  2:1576/148‑1 (MQ=255)
              aTCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGt                                                                                                                                     >  2:1143731/1‑148 (MQ=255)
              aTCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGt                                                                                                                                     <  1:2137334/148‑1 (MQ=255)
                cGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGACCAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCg                                                                                                                                   >  1:1202821/1‑148 (MQ=255)
                cGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCg                                                                                                                                   >  1:1737662/1‑148 (MQ=255)
                 gACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGc                                                                                                                                  <  2:2156035/148‑1 (MQ=255)
                  aCATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCa                                                                                                                                 <  2:443226/148‑1 (MQ=255)
                   cATCGGTACATCGCAGCTGGCGCAGGCGACCCGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAg                                                                                                                                >  1:391974/1‑148 (MQ=255)
                   cATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAg                                                                                                                                >  1:3522996/1‑148 (MQ=255)
                   cATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAg                                                                                                                                >  2:75500/1‑148 (MQ=255)
                   cATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAACTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAg                                                                                                                                <  1:553171/148‑1 (MQ=255)
                    aTCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATCTGTTCCCGCGTCGCAgg                                                                                                                               >  2:1361977/1‑148 (MQ=255)
                    aTCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                                                                                                                               >  2:1361688/1‑148 (MQ=255)
                    aTCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAgg                                                                                                                               >  1:438677/1‑148 (MQ=255)
                     tCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGt                                                                                                                              >  2:302242/1‑148 (MQ=255)
                      cGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTg                                                                                                                             >  1:693866/1‑148 (MQ=255)
                      cGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTg                                                                                                                             <  2:2198611/148‑1 (MQ=255)
                        gTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGaa                                                                                                                           <  1:506045/148‑1 (MQ=255)
                        gTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGaa                                                                                                                           <  1:2631778/148‑1 (MQ=255)
                        gTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGaa                                                                                                                           <  1:1928979/148‑1 (MQ=255)
                         tACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGaaa                                                                                                                          <  1:1664631/148‑1 (MQ=255)
                         tACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGaaa                                                                                                                          <  1:75153/148‑1 (MQ=255)
                              cGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGGGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATg                                                                                                                     >  1:3295512/1‑148 (MQ=255)
                               gCAGCTGGCGCCGGCGACCAGTGCGCCCACGACGCTGGTGGAAATTGAGAAGTACAACGGCACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGACTATGTTCCCGCGTCGCAGGTGAAATCATGa                                                                                                                    <  2:2967043/148‑1 (MQ=255)
                                 aGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGAcc                                                                                                                  <  2:2745650/148‑1 (MQ=255)
                                    tGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACg                                                                                                               <  1:2542255/148‑1 (MQ=255)
                                       cgcAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTg                                                                                                            <  2:2847848/148‑1 (MQ=255)
                                         cAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGAc                                                                                                          <  1:2330112/148‑1 (MQ=255)
                                            gCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGc                                                                                                       >  1:3318969/1‑148 (MQ=255)
                                               aCCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGcaca                                                                                                    <  2:1558071/148‑1 (MQ=255)
                                                   gTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAg                                                                                                >  2:1772090/1‑148 (MQ=255)
                                                    tGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGa                                                                                               <  1:1000020/148‑1 (MQ=255)
                                                    tGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGa                                                                                               <  2:1703206/148‑1 (MQ=255)
                                                    tGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGa                                                                                               <  2:3464122/148‑1 (MQ=255)
                                                     gcgcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGAt                                                                                              >  1:3263207/1‑148 (MQ=255)
                                                      cgcACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGAtt                                                                                             >  2:2642587/1‑148 (MQ=255)
                                                        cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGc                                                                                           <  2:1578576/148‑1 (MQ=255)
                                                        cacaCTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGc                                                                                           >  1:2765393/1‑148 (MQ=255)
                                                         acacTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGcc                                                                                          <  2:3331732/148‑1 (MQ=255)
                                                            cTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCtgt                                                                                       >  1:2496902/1‑148 (MQ=255)
                                                              aCGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTgtgt                                                                                     >  2:710520/1‑148 (MQ=255)
                                                                gCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGttt                                                                                   <  1:2717079/148‑1 (MQ=255)
                                                                 cTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTa                                                                                  <  1:752392/148‑1 (MQ=255)
                                                                      ggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCt                                                                             <  1:3544970/148‑1 (MQ=255)
                                                                        aaaTTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTgg                                                                           <  1:2397249/148‑1 (MQ=255)
                                                                         aaTTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGc                                                                          <  1:1068803/148‑1 (MQ=255)
                                                                         aaTTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCCCAGTGCGNATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGc                                                                          <  1:1807558/148‑1 (MQ=255)
                                                                           ttGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTa                                                                        >  1:356705/1‑148 (MQ=255)
                                                                            tGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAg                                                                       <  1:1287642/148‑1 (MQ=255)
                                                                            tGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAg                                                                       <  1:2070064/148‑1 (MQ=255)
                                                                                 aGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAg                                                                  >  1:2809743/1‑148 (MQ=255)
                                                                                 aGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAg                                                                  <  1:167372/148‑1 (MQ=255)
                                                                                 aGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAg                                                                  >  1:2808141/1‑148 (MQ=255)
                                                                                  gTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgt                                                                 >  2:1296851/1‑148 (MQ=255)
                                                                                     cAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgttgt                                                              <  1:1925491/148‑1 (MQ=255)
                                                                                      aaCGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAgttgtt                                                             <  2:597475/148‑1 (MQ=255)
                                                                                        cGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTtct                                                           >  1:3086754/1‑148 (MQ=255)
                                                                                          gAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTcc                                                         >  2:1554739/1‑148 (MQ=255)
                                                                                           aaCAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTccc                                                        >  2:1545041/1‑148 (MQ=255)
                                                                                            aCAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCg                                                       <  2:3137366/148‑1 (MQ=255)
                                                                                               gTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGa                                                    <  1:2633985/148‑1 (MQ=255)
                                                                                                tGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAg                                                   >  1:2669778/1‑148 (MQ=255)
                                                                                                    gCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTgg                                               >  1:1381808/1‑148 (MQ=255)
                                                                                                       ggTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACg                                            <  1:3434109/148‑1 (MQ=255)
                                                                                                       ggTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACg                                            <  2:2288720/148‑1 (MQ=255)
                                                                                                       ggTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACg                                            <  1:1827103/148‑1 (MQ=255)
                                                                                                          gACGGCGTTTAACTGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCTCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCATTTGTTCTCCCGTGAGCTGGACGATg                                         <  1:2258830/148‑1 (MQ=255)
                                                                                                          gACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATg                                         <  1:1181451/148‑1 (MQ=255)
                                                                                                           aCGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGa                                        >  2:2621154/1‑148 (MQ=255)
                                                                                                           aCGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGa                                        >  2:840247/1‑148 (MQ=255)
                                                                                                           aCGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGa                                        >  2:2224599/1‑148 (MQ=255)
                                                                                                           aCGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGa                                        >  2:1748250/1‑148 (MQ=255)
                                                                                                              gCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaaca                                     >  2:3368501/1‑148 (MQ=255)
                                                                                                               cGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaacaa                                    >  2:1520990/1‑148 (MQ=255)
                                                                                                               cGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaacaa                                    >  1:2677957/1‑148 (MQ=255)
                                                                                                                gTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGaacaac                                   <  1:2450170/148‑1 (MQ=255)
                                                                                                                 tttAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACt                                  <  2:1228503/148‑1 (MQ=255)
                                                                                                                         ccccGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCaaa                          >  2:1481665/1‑148 (MQ=255)
                                                                                                                          cccGTGGAGATGGTGGCGCAGTGCGAATCTGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAAt                         <  1:1205661/148‑1 (MQ=255)
                                                                                                                            cGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCg                       >  2:1232925/1‑148 (MQ=255)
                                                                                                                            cGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCg                       >  2:3459057/1‑148 (MQ=255)
                                                                                                                              tGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGca                     <  1:1046079/148‑2 (MQ=255)
                                                                                                                                gagaTGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAg                   >  1:2893338/1‑148 (MQ=255)
                                                                                                                                  gaTGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTg                 >  1:866931/1‑148 (MQ=255)
                                                                                                                                    tggtggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTgcg               <  2:111508/148‑1 (MQ=255)
                                                                                                                                        ggCGCAGTGCGAATATGGTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGAGTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGa           <  2:2378443/148‑1 (MQ=255)
                                                                                                                                          cgcAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATg         <  2:2401350/148‑1 (MQ=255)
                                                                                                                                           gcAGTGCTAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATgg        <  1:1143731/148‑1 (MQ=255)
                                                                                                                                           gcAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACACCACAACAGAGGCCCGGGGTGCCCGCCGGGCGAGCGCAGTTGTTCCCCCCACACCTGGACGATGAACAACACACGCAAATCGCCAGTGCGCAGATgg        <  1:1629659/148‑1 (MQ=255)
                                                                                                                                            cAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGc       >  2:3084948/1‑148 (MQ=255)
                                                                                                                                            cAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGc       <  2:106225/148‑1 (MQ=255)
                                                                                                                                                 cGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAt  >  1:3536380/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCAACGTGTTGACCATCGACATCGGTACATCGCAGCTGGCGCAGGCGACCAGTGCGCACACTACGCTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAAT  >  CP009273/1057080‑1057372

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: