Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,016 T→C 10.5% intergenic (+32/‑28) secA → / → mutT preprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,0160TC10.5% 440.0 / 34.2 153intergenic (+32/‑28)secA/mutTpreprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (65/72);  new base C (7/9);  total (72/81)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 154 total aligned reads displayed.

TGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATT  >  NC_000913/110870‑111163
                                                                                                                                                  |                                                                                                                                                   
tGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCt                                                                                                                                                   <  4:1317952/149‑1 (MQ=255)
  cAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                                                                                 <  4:863963/149‑1 (MQ=255)
   aGCCGCAGCTGCACTGGCGGCTCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGACCGTGCGGCTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgcg                                                                                                                                                >  4:1511460/1‑149 (MQ=255)
      cGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGAAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGGGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                                                                             <  2:589391/149‑1 (MQ=255)
        cAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttt                                                                                                                                           <  2:902791/149‑1 (MQ=255)
          gCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttttt                                                                                                                                         <  4:22066/149‑1 (MQ=255)
           cTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCtttttt                                                                                                                                        >  4:1412413/1‑149 (MQ=255)
               aCTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                                                                    <  2:1106054/149‑1 (MQ=255)
               aCTGGCGGCGCAAACCGAAGAGCGCAAAGGAGGCCGTAACGATACTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGTCTGCAATAAAAGCTAACGGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                                                                    <  1:1365679/149‑1 (MQ=255)
                   gcggcgCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGttt                                                                                                                                <  1:750832/149‑1 (MQ=255)
                     ggcgCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTaa                                                                                                                              <  1:929387/149‑1 (MQ=255)
                        gcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                                                                           >  2:652035/1‑149 (MQ=255)
                        gcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGAc                                                                                                                           >  1:90337/1‑149 (MQ=255)
                          aaaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACaa                                                                                                                         <  3:1504125/149‑1 (MQ=255)
                           aaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAAt                                                                                                                        >  1:497265/1‑149 (MQ=255)
                           aaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCt                                                                                                                                             <  3:229121/128‑1 (MQ=255)
                           aaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCACTGCCATGGCTGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTACAGGTTTAAGACAAt                                                                                                                        <  3:1556667/149‑1 (MQ=255)
                            aCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATg                                                                                                                       >  4:1266672/1‑149 (MQ=255)
                                   agCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGc                                                                                                                <  4:607276/149‑1 (MQ=255)
                                      gcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTGGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGc                                                                                                             <  3:906869/149‑1 (MQ=255)
                                       cAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCa                                                                                                            >  4:902864/1‑149 (MQ=255)
                                         aaGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCaaa                                                                                                          <  3:667341/149‑1 (MQ=255)
                                             aGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGc                                                                                                      <  2:465161/149‑1 (MQ=255)
                                               gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                                                                    >  2:1014963/1‑149 (MQ=255)
                                               gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCgg                                                                                                    <  3:496581/149‑1 (MQ=255)
                                                  gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtag                                                                                                 >  3:546485/1‑149 (MQ=255)
                                                  gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtag                                                                                                 <  1:1459621/149‑1 (MQ=255)
                                                  gTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtag                                                                                                 <  1:1110176/149‑1 (MQ=255)
                                                   tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGTGCAGGATTCTGCGCCTTTTTTATAGGATTAAGACAATGAAAAAGCTGCAAATTGCggtagg                                                                                                >  4:473073/1‑149 (MQ=255)
                                                    aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                               >  4:568934/1‑149 (MQ=255)
                                                    aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                               <  3:852225/149‑1 (MQ=255)
                                                    aaCGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                               >  3:408234/1‑149 (MQ=255)
                                                       gATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtatt                                                                                            >  3:813789/1‑149 (MQ=255)
                                                           cTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                                                                        >  1:648171/1‑149 (MQ=255)
                                                               cccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCaa                                                                                    <  3:456540/149‑1 (MQ=255)
                                                               cccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCaa                                                                                    <  4:546485/149‑1 (MQ=255)
                                                                ccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAAc                                                                                   >  4:255873/1‑149 (MQ=255)
                                                                  gTGCGGTTCTGGTAAAAAAGACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTGACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACga                                                                                 >  4:1085305/1‑149 (MQ=255)
                                                                     cGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGaa                                                                              <  1:236173/149‑1 (MQ=255)
                                                                       gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAAGTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACga                                                                                 >  1:1228962/1‑144 (MQ=255)
                                                                       gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACga                                                                                 <  2:1228962/144‑1 (MQ=255)
                                                                       gTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACa                                                                            >  3:327908/1‑149 (MQ=255)
                                                                        ttCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtat                                                                                                                                      >  1:181670/1‑90 (MQ=255)
                                                                         tCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAAt                                                                          <  3:255873/149‑1 (MQ=255)
                                                                         tCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAAt                                                                          <  4:330577/149‑1 (MQ=255)
                                                                         tCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAAt                                                                          <  4:303725/149‑1 (MQ=255)
                                                                         tCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAAt                                                                          <  2:737576/149‑1 (MQ=255)
                                                                             gTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaaa                                                                      <  3:1412413/149‑1 (MQ=255)
                                                                             gTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaaa                                                                      <  2:648171/149‑1 (MQ=255)
                                                                              tAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAt                                                                     <  1:684764/149‑1 (MQ=255)
                                                                              tAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAt                                                                     <  4:182443/149‑1 (MQ=255)
                                                                              tAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAt                                                                     >  2:1455390/1‑149 (MQ=255)
                                                                                aaaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCt                                                                   <  1:1092603/149‑1 (MQ=255)
                                                                                  aaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                                 >  4:1493120/1‑149 (MQ=255)
                                                                                   aaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                                                                <  2:283599/149‑1 (MQ=255)
                                                                                    aTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtat                                                               <  1:24326/149‑1 (MQ=255)
                                                                                    aTACAAGAAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtat                                                               <  4:922530/149‑1 (MQ=255)
                                                                                        aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcg                                                           >  4:1180203/1‑149 (MQ=255)
                                                                                        aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcg                                                           >  1:668264/1‑149 (MQ=255)
                                                                                        aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtt                                                                                                                                            >  4:676778/1‑68 (MQ=38)
                                                                                        aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtt                                                                                                                                            <  3:676778/68‑1 (MQ=38)
                                                                                         aGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcgc                                                          <  1:1118381/149‑1 (MQ=255)
                                                                                          gCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcgcg                                                         >  1:1368185/1‑149 (MQ=255)
                                                                                                  aTGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgcag                                                 <  4:1007073/149‑1 (MQ=255)
                                                                                                    ggCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgcagca                                               <  1:374335/149‑1 (MQ=255)
                                                                                                       cgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATg                                                                         >  1:521582/1‑120 (MQ=255)
                                                                                                        gccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATg                                           >  3:607137/1‑149 (MQ=255)
                                                                                                               aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                                                                >  3:1371239/1‑121 (MQ=255)
                                                                                                               aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACAt                                    <  4:1535077/149‑1 (MQ=255)
                                                                                                               aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACAt                                    <  3:726709/149‑1 (MQ=255)
                                                                                                               aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACAt                                    >  2:1383985/1‑149 (MQ=255)
                                                                                                               aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACAt                                    >  2:160288/1‑149 (MQ=255)
                                                                                                                tAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAAc                                                                                   >  2:745042/1‑101 (MQ=255)
                                                                                                                 aaaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATgg                                  <  2:531257/149‑1 (MQ=255)
                                                                                                                  aaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGc                                 >  1:748462/1‑149 (MQ=255)
                                                                                                                      cTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATTCGCACATGGCGAAt                             <  1:1196774/149‑1 (MQ=255)
                                                                                                                           tGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGca                                       <  3:1504354/134‑1 (MQ=255)
                                                                                                                            gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTg                       <  4:942224/149‑1 (MQ=255)
                                                                                                                            gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTg                       >  1:383855/1‑149 (MQ=255)
                                                                                                                            gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTg                       <  3:182276/149‑1 (MQ=255)
                                                                                                                            gTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTg                       >  2:102011/1‑149 (MQ=255)
                                                                                                                             ttGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                                                                >  1:189328/1‑107 (MQ=255)
                                                                                                                             ttGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                                                                <  2:189328/107‑1 (MQ=255)
                                                                                                                              tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                     <  1:62199/149‑1 (MQ=255)
                                                                                                                                aaGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGt                   >  1:1391916/1‑149 (MQ=255)
                                                                                                                                   aaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTc                >  3:1183086/2‑149 (MQ=255)
                                                                                                                                    aaaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTcc               >  1:126328/1‑149 (MQ=255)
                                                                                                                                    aaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTcc               >  1:1247493/1‑149 (MQ=255)
                                                                                                                                     aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                     >  4:1100572/1‑142 (MQ=255)
                                                                                                                                     aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                     <  3:1100572/142‑1 (MQ=255)
                                                                                                                                     aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTccc              >  2:1398515/1‑149 (MQ=255)
                                                                                                                                      aaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGTGAATAAACTGGAGTTTCCcg             >  3:471550/1‑149 (MQ=255)
                                                                                                                                      aaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcg             <  4:813789/149‑1 (MQ=255)
                                                                                                                                       aGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATaa                                                             >  4:458184/1‑100 (MQ=255)
                                                                                                                                         gcgcAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggcg          >  1:1056601/1‑149 (MQ=255)
                                                                                                                                           gcAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGt        >  4:1332334/1‑149 (MQ=255)
                                                                                                                                           gcAGGATTCTGCGCCTTTTTTATAGGGTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAAGCGTCGCGCACCAGATGCGCACAGGACGAATAAACTGGAGTTTCTCGGCGGt        >  2:298152/1‑149 (MQ=255)
                                                                                                                                             aGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTaa      >  1:701166/1‑149 (MQ=255)
                                                                                                                                               gATTCTGCGCCTTCTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTaaaa    >  1:1477778/1‑149 (MQ=255)
                                                                                                                                                 ttCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAAtt  >  2:420458/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAATT  >  NC_000913/110870‑111163

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