Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,060,123 A→T 8.2% N290Y (AAC→TAC) ‡ 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,1230AT 8.2% 249.4 / 9.1 73N290Y (AAC→TAC) ‡torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (37/30);  new base T (3/3);  total (40/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 = 9.90e-01

CCGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCG  >  NC_000913/1059982‑1060269
                                                                                                                                              |                                                                                                                                                  
ccGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCa                                                                                                                                             >  2:611250/1‑149 (MQ=255)
 cGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAAccgcaa                                                                                                                                             <  2:473027/149‑5 (MQ=255)
 cGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCaa                                                                                                                                            >  1:164129/1‑149 (MQ=255)
   gATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGGTTAACCGCaatg                                                                                                                                          <  2:411581/149‑3 (MQ=255)
     tCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGGTTAACCGCaatgtg                                                                                                                                        <  4:328204/149‑5 (MQ=255)
     tCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGGTTAACCGCaatgtg                                                                                                                                        >  3:328204/1‑145 (MQ=255)
         gATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGt                                                                                                                                    >  4:231326/1‑149 (MQ=255)
         gATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGt                                                                                                                                    >  2:334469/1‑149 (MQ=255)
           tGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGc                                                                                                                                  >  4:418695/1‑149 (MQ=255)
                aTGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTg                                                                                                                             >  4:39202/1‑149 (MQ=255)
                aTGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTg                                                                                                                             >  4:232495/1‑149 (MQ=255)
                 tGACTATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGc                                                                                                                            >  4:97500/1‑149 (MQ=255)
                    atatTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAAc                                                                                                                         >  4:584099/1‑149 (MQ=255)
                     tatTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACt                                                                                                                        >  1:127296/1‑149 (MQ=255)
                       ttACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTgg                                                                                                                      <  3:246253/149‑1 (MQ=255)
                           gcgcAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCt                                                                                                                  >  3:613484/1‑149 (MQ=255)
                             gcAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTgg                                                                                                                <  3:243190/149‑1 (MQ=255)
                                 cTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGcaca                                                                                                            <  1:167212/149‑1 (MQ=255)
                                  tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAAccgcaat                                                                                                                                            >  4:100518/1‑112 (MQ=255)
                                  tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGCGGGTTAAccgcaat                                                                                                                                            <  3:100518/117‑6 (MQ=255)
                                  tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGGTTAACCGCaat                                                                                                                                           <  3:42689/117‑2 (MQ=255)
                                  tAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCAGGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGGTTAACCGCaat                                                                                                                                           >  4:42689/1‑116 (MQ=255)
                                    aaaGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATAc                                                                                                         >  1:198409/1‑149 (MQ=255)
                                     aaGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACg                                                                                                        <  2:164129/149‑1 (MQ=255)
                                      aGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGc                                                                                                       <  4:194608/149‑1 (MQ=255)
                                       gCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCt                                                                                                      <  1:413231/149‑1 (MQ=255)
                                       gCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCt                                                                                                      <  3:584099/149‑1 (MQ=255)
                                       gCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCt                                                                                                      <  2:198409/149‑1 (MQ=255)
                                         gAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGt                                                                                                    >  1:514195/1‑149 (MQ=255)
                                          aaaGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTa                                                                                                   >  2:83363/1‑149 (MQ=255)
                                                      ggTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCt                                                                                       <  3:312793/149‑1 (MQ=255)
                                                        tGAAATTGAGGTCATCAGCATCGATCCTGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGt                                                                                     <  3:97500/149‑1 (MQ=255)
                                                        tGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGaa                                                                                            >  3:119839/1‑142 (MQ=255)
                                                        tGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGaa                                                                                            <  4:119839/142‑1 (MQ=255)
                                                          aaaTTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTAc                                                                                   <  2:328379/149‑1 (MQ=255)
                                                                aGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACaaa                                                                             <  4:512690/149‑1 (MQ=255)
                                                                       caGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTcc                                                                      <  4:79716/149‑1 (MQ=255)
                                                                       caGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTcc                                                                      >  1:109768/1‑149 (MQ=255)
                                                                         gCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                    >  2:79611/1‑149 (MQ=255)
                                                                         gCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                    <  1:181591/149‑1 (MQ=255)
                                                                              gATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTaa                                                               <  3:544004/149‑1 (MQ=255)
                                                                               aTCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAac                                                              >  3:488622/1‑149 (MQ=255)
                                                                               aTCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAac                                                              >  4:152838/1‑149 (MQ=255)
                                                                                 ccGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAacta                                                            <  2:486867/149‑1 (MQ=255)
                                                                                          acaTCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACtgtgtg                                                     >  1:491947/1‑147 (MQ=255)
                                                                                          acaTCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACtgtgtg                                                     <  2:491947/147‑1 (MQ=255)
                                                                                             tCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGttt                                                <  3:60758/149‑1 (MQ=255)
                                                                                                aCCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                    >  1:146298/1‑126 (MQ=255)
                                                                                                aCCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCtt                                                                    <  2:146298/126‑1 (MQ=255)
                                                                                                aCCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGa                                             <  3:152838/149‑1 (MQ=255)
                                                                                                  ccATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGc                                           >  3:334147/1‑149 (MQ=255)
                                                                                                  ccATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGc                                           >  3:540964/1‑149 (MQ=255)
                                                                                                   cATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGGTTAACCGCaat                                                                                                                                           >  4:197044/1‑51 (MQ=37)
                                                                                                    aTGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAg                                         >  2:347939/1‑149 (MQ=255)
                                                                                                        gTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTcc                                     <  1:347939/149‑1 (MQ=255)
                                                                                                                   gcgAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctg                          >  4:413312/1‑149 (MQ=255)
                                                                                                                   gcgAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctg                          >  2:615648/1‑149 (MQ=255)
                                                                                                                     gAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctgct                        >  4:245046/1‑149 (MQ=255)
                                                                                                                      aGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATctgctg                       <  3:231326/149‑1 (MQ=255)
                                                                                                                         aTGTGAAGCACATTGC‑GGTTAACCCGCAGACTGACGTGCCGCTGAATCTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACTACAAAAACTTCATTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGt                    >  4:381149/1‑149 (MQ=255)
                                                                                                                         aTGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTTCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGt                    <  1:60215/149‑1 (MQ=255)
                                                                                                                            tgAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATGCGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGCTgag                 >  4:386272/1‑149 (MQ=255)
                                                                                                                              aaGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGaa               >  3:275209/1‑149 (MQ=255)
                                                                                                                                  acaTTGC‑GGGTAACCCGCAAACTAACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGAc           >  4:207084/1‑149 (MQ=255)
                                                                                                                                    aTTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTc                                                                       >  4:183646/1‑87 (MQ=255)
                                                                                                                                    aTTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTc                                                                       <  3:183646/87‑1 (MQ=255)
                                                                                                                                       gC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCa      <  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)
                                                                                                                                              |                                                                                                                                                  
CCGGATCACGATGTTTATGAATATTACGCGCAGCTAAAAGCGAAAGTCGCCGCCGGTGAAATTGAGGTCATCAGCATCGATCCGGTTGTCACATCCACCCATGAGTATCTGGGGCGCGAGCATGTGAAGCACATTGC‑GGTTAACCCGCAAACTGACGTGCCGCTGCAACTGGCGCTGGCACATACGCTGTACAGTGAAAACCTGTACGACAAAAACTTCCTTGCTAACTACTGTGTGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCG  >  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: ‑