Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,055,132 G→T 94.8% 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.8% 481.4 / 14.3 166E141* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (6/3);  new base T (78/79);  total (84/82)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

TCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGAT  >  CP009273/1054987‑1055278
                                                                                                                                                 |                                                                                                                                                  
tCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATgaagaa                                                                                                                                                  <  1:1915027/148‑1 (MQ=255)
  aCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                                <  1:2876477/148‑1 (MQ=255)
   cGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTg                                                                                                                                               >  2:743153/1‑148 (MQ=255)
   cGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTg                                                                                                                                               <  1:2250796/148‑1 (MQ=255)
   cGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTg                                                                                                                                               >  2:503899/1‑148 (MQ=255)
    gTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTgg                                                                                                                                              >  1:532543/1‑148 (MQ=255)
    gTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTgg                                                                                                                                              >  2:1711902/1‑148 (MQ=255)
     tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGa                                                                                                                                             <  2:1399737/148‑1 (MQ=255)
          aaCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAAcgcg                                                                                                                                        >  2:771605/1‑148 (MQ=255)
          aaCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAAcgcg                                                                                                                                        >  2:2856271/1‑148 (MQ=255)
           actcggcgCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGt                                                                                                                                       >  1:1571485/4‑148 (MQ=255)
           aCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGt                                                                                                                                       >  1:1287265/1‑148 (MQ=255)
           aCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGt                                                                                                                                       >  1:1287566/1‑148 (MQ=255)
              cggcgCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCa                                                                                                                                    >  2:765363/1‑148 (MQ=255)
                gcgcGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGa                                                                                                                                  >  2:2932166/1‑148 (MQ=255)
                gcgcGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGa                                                                                                                                  >  2:2931735/1‑148 (MQ=255)
                   cgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaaa                                                                                                                               >  2:723116/1‑148 (MQ=255)
                   cgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaaa                                                                                                                               <  1:681254/148‑1 (MQ=255)
                       ttCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGATCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGTAACGCGTGAAGAAAACTCa                                                                                                                           <  2:1563591/148‑1 (MQ=255)
                         cGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACg                                                                                                                         <  2:3249097/148‑1 (MQ=255)
                             aTCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGACAACTCACGGGcc                                                                                                                     >  2:2021501/1‑148 (MQ=255)
                                cGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAg                                                                                                                  <  2:2607039/148‑1 (MQ=255)
                                  aTGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTg                                                                                                                <  2:451713/148‑1 (MQ=255)
                                  aTGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTg                                                                                                                >  2:931446/1‑148 (MQ=255)
                                   tGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGc                                                                                                               >  1:1517898/1‑148 (MQ=255)
                                   tGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGc                                                                                                               >  1:3287254/1‑148 (MQ=255)
                                    ggTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGcc                                                                                                              >  1:1153480/1‑148 (MQ=255)
                                     gTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCt                                                                                                             >  1:1689729/1‑148 (MQ=255)
                                       aCGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTg                                                                                                           <  1:3270623/148‑1 (MQ=255)
                                       aCGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTg                                                                                                           >  1:1041202/1‑148 (MQ=255)
                                        cgcgTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGc                                                                                                          <  1:406156/148‑1 (MQ=255)
                                               aCTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGc                                                                                                   >  1:1608989/1‑148 (MQ=255)
                                                 tGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCa                                                                                                 <  1:1731015/148‑1 (MQ=255)
                                                 tGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCa                                                                                                 <  1:1731070/148‑1 (MQ=255)
                                                  ggCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAg                                                                                                <  1:2290808/148‑1 (MQ=255)
                                                   gCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGt                                                                                               >  1:1077802/1‑148 (MQ=255)
                                                    cttcGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTGTGTGCGCGTGAGCTGTGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTg                                                                                              <  1:2458353/145‑1 (MQ=255)
                                                    cTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTg                                                                                              <  1:2760873/148‑1 (MQ=255)
                                                      gCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGt                                                                                            <  1:384469/148‑1 (MQ=255)
                                                       cGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGtt                                                                                           <  2:955403/148‑1 (MQ=255)
                                                         tAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTgg                                                                                         <  2:432160/148‑1 (MQ=255)
                                                          aaGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTATTTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGAGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGc                                                                                        >  1:1608213/1‑148 (MQ=255)
                                                            gcgcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCaa                                                                                      >  2:1913025/1‑148 (MQ=255)
                                                                 aTCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGAc                                                                                 <  1:491936/148‑1 (MQ=255)
                                                                  tctcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACg                                                                                <  2:1363311/148‑1 (MQ=255)
                                                                   ctcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCACTCGACgg                                                                               >  2:2519222/1‑148 (MQ=255)
                                                                   ctcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgg                                                                               >  1:2200707/1‑148 (MQ=255)
                                                                     cAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                                                             <  2:225665/148‑1 (MQ=255)
                                                                     cAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                                                             >  1:318545/1‑148 (MQ=255)
                                                                      aGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCTTGAGCTGGGATGAAGCCCTCGACTTGTTATTTTAATAACTGGACCGCGTGAAGAAAACTCACGGGCCGATTGCCTTGCTGACCGCCAGTGGTTGGAAATCGACGGGGa                                                                            <  1:1974042/148‑1 (MQ=255)
                                                                      aGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGa                                                                            <  1:1974351/148‑1 (MQ=255)
                                                                       gCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGAt                                                                           <  2:3493062/148‑1 (MQ=255)
                                                                          aTACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGtt                                                                        >  1:2219585/1‑148 (MQ=255)
                                                                           tACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTc                                                                       >  2:2040283/1‑148 (MQ=255)
                                                                           tACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTc                                                                       <  1:2340271/148‑1 (MQ=255)
                                                                           tACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTc                                                                       <  1:986020/148‑1 (MQ=255)
                                                                            aCCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTcc                                                                      >  2:2655249/1‑148 (MQ=255)
                                                                            aCCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTcc                                                                      >  2:493957/1‑148 (MQ=255)
                                                                             ccTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCa                                                                     >  1:1136563/1‑148 (MQ=255)
                                                                               tCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATa                                                                   <  2:7482/148‑1 (MQ=255)
                                                                               tCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATa                                                                   >  1:276344/1‑148 (MQ=255)
                                                                                cccAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATaa                                                                  <  2:3320491/148‑1 (MQ=255)
                                                                                       cgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCg                                                           <  1:3519787/148‑1 (MQ=255)
                                                                                         gTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCggg                                                         >  2:3435455/1‑148 (MQ=255)
                                                                                         gTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCggg                                                         <  1:2915779/148‑1 (MQ=255)
                                                                                         gTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCggg                                                         >  1:1202497/1‑148 (MQ=255)
                                                                                         gTGATAACCGTTTTGGGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCggg                                                         <  1:3027727/148‑1 (MQ=255)
                                                                                             tAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATg                                                     <  2:521115/148‑1 (MQ=255)
                                                                                              aaCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGc                                                    <  1:901293/148‑1 (MQ=255)
                                                                                              aaCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGc                                                    >  1:2814275/1‑148 (MQ=255)
                                                                                              aaCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGc                                                    >  1:1033183/1‑148 (MQ=255)
                                                                                              aaCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGc                                                    >  1:2702322/1‑148 (MQ=255)
                                                                                               aCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCt                                                   >  2:602368/1‑148 (MQ=255)
                                                                                                 cGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTgg                                                 >  2:2162820/1‑148 (MQ=255)
                                                                                                  gTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGc                                                <  1:2136789/148‑1 (MQ=255)
                                                                                                  gTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTATATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGc                                                >  2:1027972/1‑148 (MQ=255)
                                                                                                     ttGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGaa                                             >  1:2047126/1‑148 (MQ=255)
                                                                                                       gtgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAg                                           <  2:1109027/148‑1 (MQ=255)
                                                                                                         gcgcGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCt                                         <  1:3209202/148‑1 (MQ=255)
                                                                                                           gcgTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTAt                                       <  1:2599616/148‑1 (MQ=255)
                                                                                                               gAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGcc                                   >  2:2812737/1‑148 (MQ=255)
                                                                                                               gAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGcc                                   <  1:2644531/148‑1 (MQ=255)
                                                                                                                aGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCt                                  <  1:381012/148‑1 (MQ=255)
                                                                                                                aGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCt                                  <  2:434997/148‑1 (MQ=255)
                                                                                                                aGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCt                                  <  2:434988/148‑1 (MQ=255)
                                                                                                                    gggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                              >  2:3357827/1‑148 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                             <  1:1609654/148‑1 (MQ=255)
                                                                                                                       aTGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATgg                           >  1:2375053/1‑148 (MQ=255)
                                                                                                                        tGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGt                          >  1:3296571/1‑148 (MQ=255)
                                                                                                                         gAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTa                         >  1:600775/1‑148 (MQ=255)
                                                                                                                         gAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTa                         <  2:1276613/148‑1 (MQ=255)
                                                                                                                           aGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAAt                       <  2:1033183/148‑1 (MQ=255)
                                                                                                                           aGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAAt                       <  2:2814275/148‑1 (MQ=255)
                                                                                                                           aGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAAt                       <  2:2702322/148‑1 (MQ=255)
                                                                                                                            gCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATa                      <  1:765363/148‑1 (MQ=255)
                                                                                                                                tCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGt                  <  1:2577001/148‑1 (MQ=255)
                                                                                                                                     aTGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa             >  1:375314/1‑148 (MQ=255)
                                                                                                                                       gTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACg           <  2:27841/148‑1 (MQ=255)
                                                                                                                                                agaaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGAt  <  1:771605/148‑1 (MQ=255)
                                                                                                                                                agaaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGAt  <  1:2856271/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGAT  >  CP009273/1054987‑1055278

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: