Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,055,132 G→T 100% 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,1320GT98.1% 825.5 / 2.6 262E141* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base T (113/143);  total (115/147)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 269 total aligned reads displayed.

GATCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGAT  >  CP009273/1054985‑1055278
                                                                                                                                                   |                                                                                                                                                  
gatCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATgaataa                                                                                                                                                  <  1:3697774/150‑1 (MQ=255)
gatCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATgaataa                                                                                                                                                  >  2:2014183/1‑150 (MQ=255)
  tCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                                >  1:3534154/1‑150 (MQ=255)
      gTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGaa                                                                                                                                            >  2:1122997/1‑150 (MQ=255)
       tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAc                                                                                                                                           <  1:3352021/150‑1 (MQ=255)
       tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAc                                                                                                                                           <  2:4237879/150‑1 (MQ=255)
         cacaACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAcgc                                                                                                                                         <  2:1470777/150‑1 (MQ=255)
          acaACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAcgcg                                                                                                                                        >  1:1268500/1‑150 (MQ=255)
          acaACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAcgcg                                                                                                                                        <  2:4089785/150‑1 (MQ=255)
               gcggcgCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAg                                                                                                                                   <  2:3364625/150‑1 (MQ=255)
                  gcgcGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaa                                                                                                                                >  2:4186389/1‑150 (MQ=255)
                   cgcgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaaa                                                                                                                               <  2:3930795/150‑1 (MQ=255)
                   cgcgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaaa                                                                                                                               >  2:1544357/1‑150 (MQ=255)
                       tATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCa                                                                                                                           >  1:1773791/1‑150 (MQ=255)
                        aTTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCAc                                                                                                                          <  2:4370148/150‑1 (MQ=255)
                          tCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACgg                                                                                                                        <  1:1458562/150‑1 (MQ=255)
                          tCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACgg                                                                                                                        <  1:4389233/150‑1 (MQ=255)
                              tATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCg                                                                                                                    <  2:4940285/150‑1 (MQ=255)
                               aTCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGCTGAAGCCCTCGACATGTTCTATGAATCACTGGAACGCGTGCAGAACACTCTCGTGCCGa                                                                                                                   >  2:2639241/1‑150 (MQ=255)
                                  cGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTg                                                                                                                >  1:4415680/1‑150 (MQ=255)
                                   gATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGc                                                                                                               >  1:410237/1‑150 (MQ=255)
                                   gATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGc                                                                                                               <  1:1103686/150‑1 (MQ=255)
                                    aTGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGcc                                                                                                              <  1:4798128/150‑1 (MQ=255)
                                     tGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCt                                                                                                             <  2:1209998/150‑1 (MQ=255)
                                        tACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGc                                                                                                          >  1:1866521/1‑150 (MQ=255)
                                         aCGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCt                                                                                                         >  2:1073299/1‑150 (MQ=255)
                                           gcgTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGa                                                                                                       <  2:2331066/150‑1 (MQ=255)
                                            cgTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGAc                                                                                                      >  1:1130782/1‑150 (MQ=255)
                                               ggACTGGCTGCGTAAGCGCCATCTCAGAGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGc                                                                                                   <  1:4677592/150‑1 (MQ=255)
                                                  cTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAg                                                                                                <  2:964762/150‑1 (MQ=255)
                                                  cTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAg                                                                                                <  2:484643/150‑1 (MQ=255)
                                                   tGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGACCGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGt                                                                                               >  2:4627602/1‑150 (MQ=255)
                                                   tGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGt                                                                                               >  1:1401949/1‑150 (MQ=255)
                                                    ggCTGCGTCAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCTTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTg                                                                                              <  2:2661577/150‑1 (MQ=255)
                                                      cTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGt                                                                                            >  2:3790213/1‑150 (MQ=255)
                                                       tGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGtt                                                                                           >  2:308978/1‑150 (MQ=255)
                                                        gCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTg                                                                                          >  1:2324374/1‑150 (MQ=255)
                                                           tAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                       <  2:840329/150‑1 (MQ=255)
                                                            aaGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCaa                                                                                      >  1:3935280/1‑150 (MQ=255)
                                                               cgcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCg                                                                                   >  1:2609633/1‑150 (MQ=255)
                                                                gcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGa                                                                                  >  1:2647454/1‑150 (MQ=255)
                                                                gcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGa                                                                                  >  2:3983970/1‑150 (MQ=255)
                                                                gcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGa                                                                                  <  1:4036465/150‑1 (MQ=255)
                                                                  cATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACg                                                                                <  1:3474042/150‑1 (MQ=255)
                                                                     ctcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                                                             <  1:49348/150‑1 (MQ=255)
                                                                     ctcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                                                             >  2:516202/1‑150 (MQ=255)
                                                                          cGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGtt                                                                        >  1:1980992/1‑150 (MQ=255)
                                                                            aTACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTcc                                                                      <  2:1360324/150‑1 (MQ=255)
                                                                                  cccAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACg                                                                <  2:2236391/150‑1 (MQ=255)
                                                                                  cccAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACg                                                                <  2:1198125/150‑1 (MQ=255)
                                                                                     aGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCtt                                                             <  2:1527674/150‑1 (MQ=255)
                                                                                      gcgcGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTc                                                            <  2:4080691/150‑1 (MQ=255)
                                                                                       cgcgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTca                                                           <  1:1092081/150‑2 (MQ=255)
                                                                                        gcgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCgg                                                          >  1:1030371/1‑150 (MQ=255)
                                                                                               tAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCt                                                   >  2:794395/1‑150 (MQ=255)
                                                                                               tAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCt                                                   <  2:2812766/150‑1 (MQ=255)
                                                                                                 aCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTgg                                                 >  2:2026170/1‑150 (MQ=255)
                                                                                                   cGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCg                                               <  1:244173/150‑1 (MQ=255)
                                                                                                   cGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCg                                               <  2:2741431/150‑1 (MQ=255)
                                                                                                      tttGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGaaa                                            >  2:4643534/1‑150 (MQ=255)
                                                                                                        tgtgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGc                                          <  2:1437797/150‑1 (MQ=255)
                                                                                                        tgtgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGc                                          <  2:3338658/150‑1 (MQ=255)
                                                                                                         gtgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCt                                         <  1:3136917/150‑1 (MQ=255)
                                                                                                          tgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTa                                        <  1:2964682/150‑1 (MQ=255)
                                                                                                           gcgcGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTAt                                       <  1:1642327/150‑1 (MQ=255)
                                                                                                                 gAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGTACGCGTGCAGAAAACTCACGGGCCGATTGCATTGCTGACCGCCAGTTGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCtt                                 <  1:2699368/150‑1 (MQ=255)
                                                                                                                 gAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCtt                                 >  2:2663514/1‑150 (MQ=255)
                                                                                                                 gAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCtt                                 >  1:641892/1‑150 (MQ=255)
                                                                                                                 gAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCtt                                 <  2:3664805/150‑1 (MQ=255)
                                                                                                                      gggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATg                            <  2:4910993/150‑1 (MQ=255)
                                                                                                                      gggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATg                            <  2:290766/150‑1 (MQ=255)
                                                                                                                      gggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATg                            >  2:576437/1‑150 (MQ=255)
                                                                                                                       tgATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCATAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATgg                           <  2:4532852/149‑1 (MQ=255)
                                                                                                                       ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATgg                           <  2:1773777/150‑1 (MQ=255)
                                                                                                                       ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATgg                           <  2:5074190/150‑1 (MQ=255)
                                                                                                                        gATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGt                          >  1:4534929/1‑150 (MQ=255)
                                                                                                                        gATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGt                          >  1:3645546/1‑150 (MQ=255)
                                                                                                                        gATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGt                          <  1:2129139/150‑1 (MQ=255)
                                                                                                                          tGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTaa                        <  1:1714971/150‑1 (MQ=255)
                                                                                                                            aaGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCTAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATa                      <  2:2359855/150‑1 (MQ=255)
                                                                                                                            aaGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTCGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATa                      <  1:4957899/150‑1 (MQ=255)
                                                                                                                             aGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAg                     >  1:1315578/1‑150 (MQ=255)
                                                                                                                              gCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGc                    <  2:1806652/150‑1 (MQ=255)
                                                                                                                                ccTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGt                  >  2:4437945/1‑150 (MQ=255)
                                                                                                                                ccTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGt                  >  2:3865697/1‑150 (MQ=255)
                                                                                                                                 cTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGtt                 <  2:1625551/150‑1 (MQ=255)
                                                                                                                                 cTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGtt                 >  1:4065154/1‑150 (MQ=255)
                                                                                                                                  tCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTg                <  1:5103058/150‑1 (MQ=255)
                                                                                                                                   cGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGTATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTgg               >  1:1192150/1‑150 (MQ=255)
                                                                                                                                     aCATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa             >  2:4637896/1‑150 (MQ=255)
                                                                                                                                     aCATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa             <  1:2873969/150‑1 (MQ=255)
                                                                                                                                      cATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTAc            >  1:3391952/1‑150 (MQ=255)
                                                                                                                                      cATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTAc            <  2:4741835/150‑1 (MQ=255)
                                                                                                                                      cATGTTCTAGGAATAACGGGAACGCTTTCAGAAAACTCACGGGCCGTGTGCCTTGCTGACCGCCTGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTAc            <  1:2241384/150‑1 (MQ=255)
                                                                                                                                       aTGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACg           <  2:511142/150‑1 (MQ=255)
                                                                                                                                        tGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACgg          <  1:2026170/150‑1 (MQ=255)
                                                                                                                                         gTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCACTCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACggg         <  2:3450882/150‑1 (MQ=255)
                                                                                                                                         gTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACggg         <  2:1157365/150‑1 (MQ=255)
                                                                                                                                             tATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGga     <  1:5163073/150‑1 (MQ=255)
                                                                                                                                               tGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGgaga   >  2:4279740/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAGAT  >  CP009273/1054985‑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: