Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,060,121 T→G 9.1% V289G (GTT→GGT)  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,060,1210TG 9.1% 241.8 / 12.3 72V289G (GTT→GGT) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (35/30);  new base G (4/3);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCG  >  NC_000913/1059982‑1060269
                                                                                                                                            |                                                                                                                                                    
ccGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCa                                                                                                                                             >  2:611250/1‑149 (MQ=255)
 cGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAAccgcaa                                                                                                                                             <  2:473027/149‑5 (MQ=255)
 cGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCaa                                                                                                                                            >  1:164129/1‑149 (MQ=255)
   gATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAACCGCaatg                                                                                                                                          <  2:411581/149‑3 (MQ=255)
     tCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAACCGCaatgtg                                                                                                                                        <  4:328204/149‑5 (MQ=255)
     tCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAACCGCaatgtg                                                                                                                                        >  3:328204/1‑145 (MQ=255)
         gATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGt                                                                                                                                    >  4:231326/1‑149 (MQ=255)
         gATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGt                                                                                                                                    >  2:334469/1‑149 (MQ=255)
           tGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGc                                                                                                                                  >  4:418695/1‑149 (MQ=255)
                aTGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTg                                                                                                                             >  4:39202/1‑149 (MQ=255)
                aTGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTg                                                                                                                             >  4:232495/1‑149 (MQ=255)
                 tGACTATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGc                                                                                                                            >  4:97500/1‑149 (MQ=255)
                    atatTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAAc                                                                                                                         >  4:584099/1‑149 (MQ=255)
                     tatTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACt                                                                                                                        >  1:127296/1‑149 (MQ=255)
                       ttACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTgg                                                                                                                      <  3:246253/149‑1 (MQ=255)
                           gcgcAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCt                                                                                                                  >  3:613484/1‑149 (MQ=255)
                             gcAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTgg                                                                                                                <  3:243190/149‑1 (MQ=255)
                                 cTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGcaca                                                                                                            <  1:167212/149‑1 (MQ=255)
                                  tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAAccgcaat                                                                                                                                            >  4:100518/1‑112 (MQ=255)
                                  tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAAccgcaat                                                                                                                                            <  3:100518/117‑6 (MQ=255)
                                  tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAACCGCaat                                                                                                                                           <  3:42689/117‑2 (MQ=255)
                                  tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCAGGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAACCGCaat                                                                                                                                           >  4:42689/1‑116 (MQ=255)
                                    aaaGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATAc                                                                                                         >  1:198409/1‑149 (MQ=255)
                                     aaGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACg                                                                                                        <  2:164129/149‑1 (MQ=255)
                                      aGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGc                                                                                                       <  4:194608/149‑1 (MQ=255)
                                       gCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCt                                                                                                      <  1:413231/149‑1 (MQ=255)
                                       gCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCt                                                                                                      <  3:584099/149‑1 (MQ=255)
                                       gCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCt                                                                                                      <  2:198409/149‑1 (MQ=255)
                                         gAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGt                                                                                                    >  1:514195/1‑149 (MQ=255)
                                          aaaGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTa                                                                                                   >  2:83363/1‑149 (MQ=255)
                                                      ggTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCt                                                                                       <  3:312793/149‑1 (MQ=255)
                                                        tGAAATTGAGGTCATCAGCATCGATCCTGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGt                                                                                     <  3:97500/149‑1 (MQ=255)
                                                        tGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGaa                                                                                            >  3:119839/1‑142 (MQ=255)
                                                        tGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGaa                                                                                            <  4:119839/142‑1 (MQ=255)
                                                          aaaTTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTAc                                                                                   <  2:328379/149‑1 (MQ=255)
                                                                aGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACaaa                                                                             <  4:512690/149‑1 (MQ=255)
                                                                       caGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTcc                                                                      <  4:79716/149‑1 (MQ=255)
                                                                       caGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTcc                                                                      >  1:109768/1‑149 (MQ=255)
                                                                         gCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                    >  2:79611/1‑149 (MQ=255)
                                                                         gCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                    <  1:181591/149‑1 (MQ=255)
                                                                              gATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTaa                                                               <  3:544004/149‑1 (MQ=255)
                                                                               aTCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAac                                                              >  3:488622/1‑149 (MQ=255)
                                                                               aTCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAac                                                              >  4:152838/1‑149 (MQ=255)
                                                                                 ccGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAacta                                                            <  2:486867/149‑1 (MQ=255)
                                                                                          acaTCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACtgtgtg                                                     >  1:491947/1‑147 (MQ=255)
                                                                                          acaTCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACtgtgtg                                                     <  2:491947/147‑1 (MQ=255)
                                                                                             tCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGttt                                                <  3:60758/149‑1 (MQ=255)
                                                                                                aCCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                    >  1:146298/1‑126 (MQ=255)
                                                                                                aCCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                    <  2:146298/126‑1 (MQ=255)
                                                                                                aCCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGa                                             <  3:152838/149‑1 (MQ=255)
                                                                                                  ccATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGc                                           >  3:334147/1‑149 (MQ=255)
                                                                                                  ccATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGc                                           >  3:540964/1‑149 (MQ=255)
                                                                                                   cATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAACCGCaat                                                                                                                                           >  4:197044/1‑51 (MQ=37)
                                                                                                    aTGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAg                                         >  2:347939/1‑149 (MQ=255)
                                                                                                        gTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTcc                                     <  1:347939/149‑1 (MQ=255)
                                                                                                                   gcgAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctg                          >  4:413312/1‑149 (MQ=255)
                                                                                                                   gcgAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctg                          >  2:615648/1‑149 (MQ=255)
                                                                                                                     gAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctgct                        >  4:245046/1‑149 (MQ=255)
                                                                                                                      aGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctgctg                       <  3:231326/149‑1 (MQ=255)
                                                                                                                         aTGTGAAGCACATTGCGGTTAACCCGCAGACTGACGTGCCGCTGAATCTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACTACAAAAACTTCATTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGt                    >  4:381149/1‑149 (MQ=255)
                                                                                                                         aTGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTTCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGt                    <  1:60215/149‑1 (MQ=255)
                                                                                                                            tgAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATGCGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGCTgag                 >  4:386272/1‑149 (MQ=255)
                                                                                                                              aaGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGaa               >  3:275209/1‑149 (MQ=255)
                                                                                                                                  acaTTGCGGGTAACCCGCAAACTAACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGAc           >  4:207084/1‑149 (MQ=255)
                                                                                                                                    aTTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTc                                                                       >  4:183646/1‑87 (MQ=255)
                                                                                                                                    aTTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTc                                                                       <  3:183646/87‑1 (MQ=255)
                                                                                                                                       gCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCa      <  4:613484/149‑1 (MQ=255)
                                                                                                                                          ggTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGTTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGc    >  2:17275/1‑149 (MQ=255)
                                                                                                                                          ggTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACgg         >  4:561666/1‑144 (MQ=255)
                                                                                                                                          ggTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACgg         <  3:561666/144‑1 (MQ=255)
                                                                                                                                          ggTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGc    <  1:83363/149‑1 (MQ=255)
                                                                                                                                           gTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGcc   >  4:67143/1‑149 (MQ=255)
                                                                                                                                            ttAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCg  >  4:515641/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CCGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCG  >  NC_000913/1059982‑1060269

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