Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,055,132 G→T 94.1% E141* (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,055,1320GT94.1% 462.3 / 17.8 163E141* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (5/5);  new base T (80/73);  total (85/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

TCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAG  >  CP009273/1054987‑1055276
                                                                                                                                                 |                                                                                                                                                
tCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATgaataa                                                                                                                                                >  2:642579/1‑148 (MQ=255)
 cACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAAc                                                                                                                                               >  1:3292735/1‑148 (MQ=255)
 cACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAAc                                                                                                                                               >  2:403636/1‑148 (MQ=255)
  aCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                              <  2:2309507/148‑1 (MQ=255)
  aCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                              <  2:2309516/148‑1 (MQ=255)
    gTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTgg                                                                                                                                            >  1:2461360/1‑148 (MQ=255)
     tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGTCTGGCTGCGTAAGCGCCATCTCAGCGATACGTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGa                                                                                                                                           <  2:622103/148‑1 (MQ=255)
     tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGa                                                                                                                                           <  2:622002/148‑1 (MQ=255)
     tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGa                                                                                                                                           >  2:648274/1‑148 (MQ=255)
       cacaACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAc                                                                                                                                         >  1:2867779/1‑148 (MQ=255)
           aCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGt                                                                                                                                     >  1:1100904/1‑148 (MQ=255)
           aCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGt                                                                                                                                     <  2:2190752/148‑1 (MQ=255)
            cgcgGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTg                                                                                                                                    >  1:1247048/1‑148 (MQ=255)
                gcgcGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGa                                                                                                                                >  2:1636318/1‑148 (MQ=255)
                  gcgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaa                                                                                                                              <  2:2471216/148‑1 (MQ=255)
                   cgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaaa                                                                                                                             <  1:1624749/148‑1 (MQ=255)
                   cgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTGTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaaa                                                                                                                             >  2:3062866/1‑148 (MQ=255)
                     tATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACt                                                                                                                           >  1:590999/1‑148 (MQ=255)
                      aTTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTc                                                                                                                          >  1:2276162/1‑148 (MQ=255)
                      aTTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTc                                                                                                                          >  1:3304356/1‑148 (MQ=255)
                      aTTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTc                                                                                                                          >  1:84291/1‑148 (MQ=255)
                        tCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCAc                                                                                                                        >  2:2904584/1‑148 (MQ=255)
                        tCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCAc                                                                                                                        >  2:2904553/1‑148 (MQ=255)
                         cGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACg                                                                                                                       <  1:2301270/148‑1 (MQ=255)
                               ccGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGa                                                                                                                 >  1:2119265/1‑148 (MQ=255)
                               ccGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGa                                                                                                                 >  1:494185/1‑148 (MQ=255)
                               ccGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATATCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGa                                                                                                                 <  1:1534694/148‑1 (MQ=255)
                                  aTGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTg                                                                                                              <  2:2113856/148‑1 (MQ=255)
                                  aTGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTg                                                                                                              >  2:1035900/1‑148 (MQ=255)
                                    ggTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGcc                                                                                                            <  1:366590/148‑1 (MQ=255)
                                       aCGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTg                                                                                                         <  2:1930755/148‑1 (MQ=255)
                                       aCGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTg                                                                                                         >  2:534652/1‑148 (MQ=255)
                                       aCGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTg                                                                                                         <  1:599452/148‑1 (MQ=255)
                                        cgcgTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCCGAAAACTCACGGGCCGAGTGCCTTGc                                                                                                        >  2:297562/1‑148 (MQ=255)
                                              gACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCg                                                                                                  >  2:3107040/1‑148 (MQ=255)
                                               aCTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGc                                                                                                 >  1:2860790/1‑148 (MQ=255)
                                               aCTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGc                                                                                                 <  1:2935161/148‑1 (MQ=255)
                                               aCTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGc                                                                                                 <  1:2935145/148‑1 (MQ=255)
                                               aCTGGCTGCGTAAGCGCCATCTCAGCGATACCACCAAGAGCGGTGATAACCGATTTGTGCGCGTGAGCTGGGAGGAATCCCTCGACATGTACTATGAATAACTGGAACGCGAGCAGAAAACAAACGTGACGAGTGCCTGGCTGACCGc                                                                                                 >  2:933644/1‑148 (MQ=255)
                                                    cTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTg                                                                                            >  1:87703/1‑148 (MQ=255)
                                                     tGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTgg                                                                                           <  2:2748252/148‑1 (MQ=255)
                                                         tAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTgg                                                                                       <  1:2312719/148‑1 (MQ=255)
                                                         tAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTgg                                                                                       >  1:837963/1‑148 (MQ=255)
                                                           aGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                     >  2:671852/1‑148 (MQ=255)
                                                           aGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                     >  1:1263656/1‑148 (MQ=255)
                                                           aGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                     >  1:1263619/1‑148 (MQ=255)
                                                           aGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                     >  2:671856/1‑148 (MQ=255)
                                                           aGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                     >  2:2334405/1‑148 (MQ=255)
                                                           aGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                     >  2:2569141/1‑148 (MQ=255)
                                                            gcgcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCaa                                                                                    >  1:158607/1‑148 (MQ=255)
                                                               ccATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCg                                                                                 <  1:270458/148‑1 (MQ=255)
                                                               ccATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCg                                                                                 <  1:1412655/148‑1 (MQ=255)
                                                                      aGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGa                                                                          <  1:3253181/148‑1 (MQ=255)
                                                                            aCCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTcc                                                                    >  1:2191056/1‑148 (MQ=255)
                                                                              cTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCAt                                                                  <  2:3069632/148‑1 (MQ=255)
                                                                               tCCCCGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATa                                                                 <  2:546746/148‑1 (MQ=255)
                                                                               tCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATa                                                                 >  2:1497085/1‑148 (MQ=255)
                                                                                  cAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACg                                                              >  2:1126786/1‑148 (MQ=255)
                                                                                     cgcgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCtt                                                           >  2:1478211/1‑148 (MQ=255)
                                                                                        ggTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCgg                                                        <  1:3175537/148‑1 (MQ=255)
                                                                                          tGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCgggg                                                      >  1:1472735/1‑148 (MQ=255)
                                                                                           gATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGa                                                     >  1:2867537/1‑148 (MQ=255)
                                                                                            aTAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGAt                                                    >  1:967374/1‑148 (MQ=255)
                                                                                             tAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATg                                                   <  1:3131928/148‑1 (MQ=255)
                                                                                                ccGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTg                                                >  1:1506606/1‑148 (MQ=255)
                                                                                                ccGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTg                                                <  1:2279650/148‑1 (MQ=255)
                                                                                                  gTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGc                                              <  1:384588/148‑1 (MQ=255)
                                                                                                  gTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGc                                              <  1:2052186/148‑1 (MQ=255)
                                                                                                    tttGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGa                                            >  1:530164/1‑148 (MQ=255)
                                                                                                         gcgcGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCt                                       <  2:2414806/148‑1 (MQ=255)
                                                                                                           gcgTGAGCTGGGATGACGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCTGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTAt                                     <  1:967928/148‑1 (MQ=255)
                                                                                                             gTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTg                                   <  2:2706071/148‑1 (MQ=255)
                                                                                                             gTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTg                                   <  1:648274/148‑1 (MQ=255)
                                                                                                              tGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGc                                  <  1:1800221/148‑1 (MQ=255)
                                                                                                              tGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGc                                  <  1:403827/148‑1 (MQ=255)
                                                                                                               gAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGcc                                 <  2:1472735/148‑1 (MQ=255)
                                                                                                               gAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGcc                                 >  2:529666/1‑148 (MQ=255)
                                                                                                                  cTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTg                              <  1:762197/148‑1 (MQ=255)
                                                                                                                   tGGGATGAAGCCCTCGACATGTTCTAGGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGc                             <  2:1438540/148‑1 (MQ=255)
                                                                                                                    gggagccaGCCCGCGACCTGTTCTACGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                            <  1:3200637/141‑1 (MQ=255)
                                                                                                                    gggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                            <  2:494185/148‑1 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                           <  1:3107040/148‑1 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                           >  2:1118205/1‑148 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                           <  2:2243740/148‑1 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                           <  2:2243741/148‑1 (MQ=255)
                                                                                                                           aGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAAt                     >  2:2665154/1‑148 (MQ=255)
                                                                                                                           aGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAAt                     >  2:2665073/1‑148 (MQ=255)
                                                                                                                                   aCATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTgg             <  2:2276162/148‑1 (MQ=255)
                                                                                                                                   aCATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTgg             <  2:3258283/148‑1 (MQ=255)
                                                                                                                                   aCATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTgg             <  2:590999/148‑1 (MQ=255)
                                                                                                                                   aCATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTgg             <  2:591011/148‑1 (MQ=255)
                                                                                                                                     aTGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa           <  2:1770346/148‑1 (MQ=255)
                                                                                                                                     aTGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa           <  2:1770182/148‑1 (MQ=255)
                                                                                                                                     aTGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa           <  1:1033841/148‑1 (MQ=255)
                                                                                                                                      tGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTAc          <  1:534652/148‑1 (MQ=255)
                                                                                                                                      tGTTCTATGAATAACTGGAACCCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTAc          >  2:207993/1‑148 (MQ=255)
                                                                                                                                       gTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACg         <  1:642579/148‑1 (MQ=255)
                                                                                                                                            aTGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCgg    <  1:2042251/148‑1 (MQ=255)
                                                                                                                                              gaataaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGgag  >  2:867297/1‑148 (MQ=255)
                                                                                                                                              gaagaaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGGGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGgag  >  2:431406/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAG  >  CP009273/1054987‑1055276

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: