Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,061,284 G→T 7.4% E677* (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
*NC_0009131,061,2840GT 7.4% 285.4 / 7.3 81E677* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (30/45);  new base T (3/3);  total (33/48)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCC  >  NC_000913/1061146‑1061420
                                                                                                                                          |                                                                                                                                        
cGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGa                                                                                                                                         <  1:682706/140‑1 (MQ=255)
cGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGa                                                                                                                                         <  1:434140/140‑1 (MQ=255)
cGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGa                                                                                                                                         <  1:434074/140‑1 (MQ=255)
   caccaGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACg                                                                                                                                      <  2:8861/140‑1 (MQ=255)
       aGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTcc                                                                                                                                  <  2:175218/140‑1 (MQ=255)
       aGGCATTCCGCGAAGATCCGGATATCGAACCGCTGGGAACGCCGAGTGGCCTGATTGATATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTcc                                                                                                                                  <  2:159486/140‑1 (MQ=255)
             tCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTATGAGAAAATCGAACGCTCCCACgct                                                                                                                            >  1:1634146/1‑138 (MQ=255)
                         cGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCg                                                                                                                >  1:1206807/1‑140 (MQ=255)
                         cGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCg                                                                                                                <  1:586256/140‑1 (MQ=255)
                         cGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCg                                                                                                                >  1:1554522/1‑140 (MQ=255)
                                 gAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTa                                                                                                        <  1:1274975/140‑1 (MQ=255)
                                    ccGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATcc                                                                                                     >  1:956761/1‑140 (MQ=255)
                                      gCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGt                                                                                                   <  1:1259329/140‑1 (MQ=255)
                                        tGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTg                                                                                                 <  1:481725/140‑1 (MQ=255)
                                         gggCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGc                                                                                                <  2:379489/140‑1 (MQ=255)
                                           gCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCAt                                                                                              <  2:182232/140‑1 (MQ=255)
                                                    gTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCt                                                                                     <  2:1284248/140‑1 (MQ=255)
                                                        ccTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGc                                                                                 <  2:906642/140‑1 (MQ=255)
                                                        ccTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGc                                                                                 <  1:1012689/140‑1 (MQ=255)
                                                         cTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCa                                                                                >  2:1194663/1‑140 (MQ=255)
                                                          tGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCAt                                                                               >  1:616635/1‑140 (MQ=255)
                                                          tGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCAt                                                                               >  1:616724/1‑140 (MQ=255)
                                                           gATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATc                                                                              >  2:418324/1‑140 (MQ=255)
                                                            aTTGGGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATcc                                                                             <  2:1276932/140‑1 (MQ=255)
                                                            aTTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATcc                                                                             <  2:1276912/140‑1 (MQ=255)
                                                             ttGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCg                                                                            >  1:726810/1‑140 (MQ=255)
                                                                  aTCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCTAACGCTCCCACgg                                                                                                                             <  2:637382/86‑1 (MQ=255)
                                                                  aTCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCTAACGCTCCCACgg                                                                                                                             <  2:639659/86‑1 (MQ=255)
                                                                  aTCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCTAACGCTCCCACgg                                                                                                                             <  2:639807/86‑1 (MQ=255)
                                                                  aTCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCTAACGCTCCCACgg                                                                                                                             >  1:637382/1‑86 (MQ=255)
                                                                  aTCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCTAACGCTCCCACgg                                                                                                                             >  1:639659/1‑86 (MQ=255)
                                                                  aTCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCTAACGCTCCCACgg                                                                                                                             >  1:639807/1‑86 (MQ=255)
                                                                  aTCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATAGAACGCTACCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGAttt                                                                       >  1:1244336/1‑140 (MQ=255)
                                                                    cTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTcc                                                                     <  2:813645/140‑1 (MQ=255)
                                                                      aCTCGAAAACTATCGCCTATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCTTCTGCAATCTGTGCCTCCGGAATTCCGa                                                                   >  2:1683750/1‑140 (MQ=255)
                                                                      aCTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGa                                                                   >  1:545452/1‑140 (MQ=255)
                                                                      aCTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGa                                                                   >  1:575875/1‑140 (MQ=255)
                                                                      aCTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGa                                                                   >  1:1017028/1‑140 (MQ=255)
                                                                      aCTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGa                                                                   >  1:1015546/1‑140 (MQ=255)
                                                                       cTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGAc                                                                  >  1:76742/1‑140 (MQ=255)
                                                                       cTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGAc                                                                  >  1:590228/1‑140 (MQ=255)
                                                                       cTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGAc                                                                  >  1:501871/1‑140 (MQ=255)
                                                                       cTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCACCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGAc                                                                  >  1:336769/1‑140 (MQ=255)
                                                                        tCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACt                                                                 >  1:1463708/1‑140 (MQ=255)
                                                                          gAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTc                                                               >  2:70834/1‑140 (MQ=255)
                                                                            aaaCTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCAc                                                             >  2:949903/1‑140 (MQ=255)
                                                                            aaaCTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGATTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCAc                                                             <  2:427637/140‑1 (MQ=255)
                                                                             aaCTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACt                                                            <  1:1402245/140‑1 (MQ=255)
                                                                             aaCTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACt                                                            <  2:616724/140‑1 (MQ=255)
                                                                             aaCTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACt                                                            <  2:616635/140‑1 (MQ=255)
                                                                             aaCTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACt                                                            <  1:1470080/140‑1 (MQ=255)
                                                                              aCTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTc                                                           <  2:905623/140‑1 (MQ=255)
                                                                                 aTCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCa                                                        <  1:645353/140‑1 (MQ=255)
                                                                                  tCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAg                                                       <  1:634227/140‑1 (MQ=255)
                                                                                       gATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTAtg                                                  <  2:963846/140‑1 (MQ=255)
                                                                                                tACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGa                                         >  2:951799/1‑140 (MQ=255)
                                                                                                       aTTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTg                                  >  1:1218206/1‑140 (MQ=255)
                                                                                                       aTTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTg                                  >  1:1296621/1‑140 (MQ=255)
                                                                                                          gTCAGGGGCATCCGATGTGGTTTGAGAATATCGATCGCTCCCACGGTGGGCCTGGCTCGCAAATGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCTCTCGCAGTTATGTGAGTCGGTAACGCTGCGt                               >  2:316634/1‑140 (MQ=255)
                                                                                                          gTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGt                               <  1:667289/140‑1 (MQ=255)
                                                                                                          gTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGt                               <  1:1116596/140‑1 (MQ=255)
                                                                                                          gTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGt                               <  1:55489/140‑1 (MQ=255)
                                                                                                          gTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGt                               <  1:958764/140‑1 (MQ=255)
                                                                                                          gTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCAGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGt                               >  1:1607804/1‑140 (MQ=255)
                                                                                                             aGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAg                            <  1:1208861/140‑1 (MQ=255)
                                                                                                             aGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAg                            <  2:1306952/140‑1 (MQ=255)
                                                                                                               tggCATCCGATGTGGTTTGAGAAAATCGAACGCACCCCCGGTGGTCCTGGCTCGCAAATGTATCCGTTGCATCTGCAATCTTTGCAACCGGATTTCCGACTTCAATCGGAGTTATGTGAGTCGGAAACGCTGCGTCAGCa                          <  2:871062/139‑1 (MQ=255)
                                                                                                               gggCATGCGATGTGTTTTGAGATATTCGTACGCTCCCACGGTGGGCCTGTCTCTCTCAAGTATCCGTTGGATCTGCAATCTGTGCTTCAGGTTTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCa                          <  2:420210/140‑1 (MQ=255)
                                                                                                               gggCATCCGATGTGGTTTGAGATAATCGATCGCTCCCACGGTGGGCCTGGCTCGCAACACTATCCGTTGCATCTGCTATGTGTGCATCCGGATTTCCGACTTCAGTCGCAGTTATGTTAGTCGGAAACGCTGCGTAAGCa                          <  2:853965/140‑1 (MQ=255)
                                                                                                               gggCATACGATGTGGATTTAGAAAATCGATCGCTCCCAAGGTGGGCATGGCTCGAACAAGTATTAGTTGCATCTGCAATCTGTGCATCCGGATTTCCGTGTTCACTCGCAGTTATGTGAGTCGGAATCGCTGCGTCAGCa                          <  2:1668002/140‑1 (MQ=255)
                                                                                                                ggCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCaa                         <  2:286648/140‑1 (MQ=255)
                                                                                                                      cGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACg                   <  1:1490376/140‑1 (MQ=255)
                                                                                                                          gtgGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAg               <  2:1015546/140‑1 (MQ=255)
                                                                                                                          gtgGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAg               <  2:1017028/140‑1 (MQ=255)
                                                                                                                          gtgGATTGAGAAAATCGAACGCTCCCACGGTGGGCCTGTCTCGCAAAAGTATCCGTTTCATCAGCAATCTGTGCATCCGGATTTCCGAATTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAg               <  1:316634/140‑1 (MQ=255)
                                                                                                                                tGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTa         <  2:1047989/140‑1 (MQ=255)
                                                                                                                                  agaAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTaaa       >  2:629583/1‑140 (MQ=255)
                                                                                                                                   gaAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCTATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAag      <  2:1214030/140‑1 (MQ=255)
                                                                                                                                      aaTCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGc   <  1:198228/140‑1 (MQ=255)
                                                                                                                                       aTCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGcc  >  2:120405/1‑140 (MQ=255)
                                                                                                                                       aTCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGcc  >  2:1216488/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                        
CGCCACCAGGCATTCCGCGAAGATCCGGATCTCGAACCGCTGGGCACGCCGAGTGGCCTGATTGAGATCTACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCC  >  NC_000913/1061146‑1061420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: