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,1320GT97.5% 608.3 / 6.2 199E141* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (88/106);  total (91/108)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 201 total aligned reads displayed.

TCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGG  >  CP009273/1054987‑1055274
                                                                                                                                                 |                                                                                                                                              
tCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATgaataa                                                                                                                                              >  1:3455297/1‑148 (MQ=255)
 cACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAAc                                                                                                                                             <  1:3165278/148‑1 (MQ=255)
 cACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAAc                                                                                                                                             >  2:2533341/1‑148 (MQ=255)
 cACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAAc                                                                                                                                             <  1:821433/148‑1 (MQ=255)
  aCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                            >  2:362618/1‑148 (MQ=255)
     tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGa                                                                                                                                         >  1:1666198/1‑148 (MQ=255)
      acacAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGaa                                                                                                                                        >  1:2720995/1‑148 (MQ=255)
       cacaACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAc                                                                                                                                       >  1:2253262/1‑148 (MQ=255)
           aCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGt                                                                                                                                   >  2:656236/1‑148 (MQ=255)
            cgcgGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTg                                                                                                                                  <  2:1789443/148‑1 (MQ=255)
                gcgcGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGa                                                                                                                              >  1:1606863/1‑148 (MQ=255)
                   cgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaaa                                                                                                                           <  2:1553699/148‑1 (MQ=255)
                   cgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTAGACATGTTCTATGAATAACTGGAACGCGTGCAGaaaa                                                                                                                           >  2:3526050/1‑148 (MQ=255)
                      aTTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTc                                                                                                                        <  1:3097051/148‑1 (MQ=255)
                           ttATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACggg                                                                                                                   <  2:2016104/148‑1 (MQ=255)
                            tATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCGATGAATAACTGGAACGCGTGCAGAAAACTCACGGGc                                                                                                                  >  1:2397038/1‑148 (MQ=255)
                              tCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCg                                                                                                                >  2:3442826/1‑148 (MQ=255)
                               ccGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGa                                                                                                               >  2:3017862/1‑148 (MQ=255)
                                    ggTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTGCTCTGAATACCTGGTACGCGTGCAGAACACTCCCGGGCCGGGTGcg                                                                                                          >  2:1845411/1‑147 (MQ=255)
                                      tACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCtt                                                                                                        <  1:1847827/148‑1 (MQ=255)
                                           gTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGa                                                                                                   <  1:1530522/148‑1 (MQ=255)
                                            tGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGAc                                                                                                  <  2:1012619/148‑1 (MQ=255)
                                            tGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGAc                                                                                                  >  2:2266882/1‑148 (MQ=255)
                                              gACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCg                                                                                                >  2:2836565/1‑148 (MQ=255)
                                               aCTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGc                                                                                               >  1:3480765/1‑148 (MQ=255)
                                                    cTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTg                                                                                          >  2:341156/1‑148 (MQ=255)
                                                     tGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTgg                                                                                         >  2:639623/1‑148 (MQ=255)
                                                       cGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGtt                                                                                       >  2:2051747/1‑148 (MQ=255)
                                                        gTAAGCGCCATCTCAGCGATACCTCCCCGCGCGGTGATAACCGGTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGACCGCGTGCAGAAAACTCTCGGGCCGAGTGCGTTGCTGACCGCCCGGGGTTg                                                                                      >  2:3359403/1‑148 (MQ=255)
                                                         tAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTgg                                                                                     <  1:2907976/148‑1 (MQ=255)
                                                          aaGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGc                                                                                    >  2:3342993/1‑148 (MQ=255)
                                                          aaGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGc                                                                                    >  2:497802/1‑148 (MQ=255)
                                                          aaGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGc                                                                                    >  2:2599574/1‑148 (MQ=255)
                                                            gcgcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCaa                                                                                  >  2:788628/1‑148 (MQ=255)
                                                             cgcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGGGGTTGGCAAt                                                                                 >  2:2411912/1‑148 (MQ=255)
                                                               ccATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCg                                                                               <  2:3247845/148‑1 (MQ=255)
                                                               ccATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCg                                                                               <  2:477981/148‑1 (MQ=255)
                                                               ccATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCg                                                                               <  2:1666198/148‑1 (MQ=255)
                                                               ccATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCg                                                                               >  2:1580099/1‑148 (MQ=255)
                                                                cATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCACTCGa                                                                              >  2:3477604/1‑148 (MQ=255)
                                                                 aTCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGAc                                                                             >  2:538193/1‑148 (MQ=255)
                                                                  tctcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACg                                                                            <  1:2030921/148‑1 (MQ=255)
                                                                  tctcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACg                                                                            <  1:1746876/148‑1 (MQ=255)
                                                                    tcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTGGCTGACCGCCAGTGGTTGGCAATCGACggg                                                                          <  1:3090701/148‑1 (MQ=255)
                                                                      aGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGa                                                                        <  2:1114231/148‑1 (MQ=255)
                                                                       gggATACCTCCCCGCGCGGTGATAACCTTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACTTGTTCTATGAATAACGGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGAt                                                                       <  1:1283563/146‑1 (MQ=255)
                                                                         gATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGt                                                                     <  1:2846089/148‑1 (MQ=255)
                                                                            aCCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTcc                                                                  <  1:3108051/148‑1 (MQ=255)
                                                                              cgccccGCGCGGTGATAACCGTTTTGGGCGCGTGAGCTGGGATGAAGCCCTCGACATTTTATATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCAt                                                                <  1:1504430/146‑1 (MQ=255)
                                                                               tCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATa                                                               <  1:1059580/148‑1 (MQ=255)
                                                                                cccAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATaa                                                              >  1:2096329/1‑148 (MQ=255)
                                                                                  cAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACg                                                            <  2:1764867/148‑1 (MQ=255)
                                                                                   aGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGc                                                           <  2:2820620/148‑1 (MQ=255)
                                                                                   aGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGc                                                           <  2:2820643/148‑1 (MQ=255)
                                                                                   aGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGc                                                           <  1:2416124/148‑1 (MQ=255)
                                                                                   aGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGc                                                           <  2:3539102/148‑1 (MQ=255)
                                                                                    gcgcGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGCCGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCt                                                          <  2:1134248/148‑1 (MQ=255)
                                                                                         gTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCggg                                                     >  1:1313933/1‑148 (MQ=255)
                                                                                          tGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCgggg                                                    <  2:2534889/148‑1 (MQ=255)
                                                                                          tGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCgggg                                                    >  1:978448/1‑148 (MQ=255)
                                                                                           gATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGa                                                   >  1:3183049/1‑148 (MQ=255)
                                                                                             tAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATTTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATg                                                 <  1:332754/148‑1 (MQ=255)
                                                                                             tAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATg                                                 <  2:1461421/148‑1 (MQ=255)
                                                                                               aCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCt                                               >  1:831253/1‑148 (MQ=255)
                                                                                               aCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCt                                               <  1:2542275/148‑1 (MQ=255)
                                                                                                ccGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTg                                              <  1:798224/148‑1 (MQ=255)
                                                                                                ccGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTg                                              <  1:295423/148‑1 (MQ=255)
                                                                                                ccGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTg                                              >  1:1004540/1‑148 (MQ=255)
                                                                                                 cGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTgg                                             >  2:3014885/1‑148 (MQ=255)
                                                                                                    tttGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGa                                          >  1:2684545/1‑148 (MQ=255)
                                                                                                    tttGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGa                                          >  2:903019/1‑148 (MQ=255)
                                                                                                       gtgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAg                                       <  2:490768/148‑1 (MQ=255)
                                                                                                          cgcgTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTa                                    <  1:594928/148‑1 (MQ=255)
                                                                                                           gcgTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTAt                                   <  2:10707/148‑1 (MQ=255)
                                                                                                            cgTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTAtt                                  <  1:1678514/148‑1 (MQ=255)
                                                                                                             gTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAACGCTATTg                                 <  2:2471481/148‑1 (MQ=255)
                                                                                                             gTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTg                                 <  2:2471574/148‑1 (MQ=255)
                                                                                                             gTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGAAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTg                                 <  1:3332683/148‑1 (MQ=255)
                                                                                                              tGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGc                                <  1:1483810/148‑1 (MQ=255)
                                                                                                              tGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGc                                <  1:1788309/148‑1 (MQ=255)
                                                                                                              tGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGc                                >  1:2827736/1‑148 (MQ=255)
                                                                                                                   tGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGc                           <  2:2157198/148‑1 (MQ=255)
                                                                                                                    gggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                          <  1:3017681/148‑1 (MQ=255)
                                                                                                                    gggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                          <  1:2836525/148‑1 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                         >  1:2531686/1‑148 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                         >  2:2413509/1‑148 (MQ=255)
                                                                                                                      gATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATg                        <  2:2694005/148‑1 (MQ=255)
                                                                                                                       aTGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATgg                       >  1:380859/1‑148 (MQ=255)
                                                                                                                          aaGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTaa                    <  1:2031825/148‑1 (MQ=255)
                                                                                                                          aaGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTaa                    <  1:1459865/148‑1 (MQ=255)
                                                                                                                           aGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAAt                   <  1:2282006/148‑1 (MQ=255)
                                                                                                                            gCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAAAGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATa                  >  2:3204024/1‑148 (MQ=255)
                                                                                                                               cTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCg               <  1:656219/148‑1 (MQ=255)
                                                                                                                                 tgACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGtt             <  1:3526050/147‑1 (MQ=255)
                                                                                                                                    cATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGt          >  2:589643/1‑148 (MQ=255)
                                                                                                                                    cATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGt          <  2:254536/148‑1 (MQ=255)
                                                                                                                                     aTGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTa         >  1:929528/1‑148 (MQ=255)
                                                                                                                                       gTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACg       >  2:497001/1‑148 (MQ=255)
                                                                                                                                         tCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACggg     >  2:3093744/1‑148 (MQ=255)
                                                                                                                                            aTGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTCATAGCGTTGGTACGGGCgg  >  1:3022627/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGG  >  CP009273/1054987‑1055274

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: