Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,013 G→A 8.4% intergenic (+29/‑31) 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,0130GA 8.4% 356.7 / 15.0 119intergenic (+29/‑31)secA/mutTpreprotein translocase subunit, ATPase/dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (57/51);  new base A (6/4);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 123 total aligned reads displayed.

CGACTCTGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAAT  >  NC_000913/110864‑111162
                                                                                                                                                     |                                                                                                                                                     
cgacTCTGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAgg                                                                                                                                                       >  1:174653/1‑150 (MQ=255)
    tctGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTc                                                                                                                                                   <  2:1769285/150‑1 (MQ=255)
       gCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                                                                                <  1:623199/150‑1 (MQ=255)
       gCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTgc                                                                                                                                                <  2:75379/150‑1 (MQ=255)
         aGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAAGAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGTTTCTgcgc                                                                                                                                              >  1:334648/1‑150 (MQ=255)
          gCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGcc                                                                                                                                             >  1:102648/1‑150 (MQ=255)
            cGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtt                                                                                                                                           <  2:1155686/150‑1 (MQ=255)
             gCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttt                                                                                                                                          <  2:949518/150‑1 (MQ=255)
             gCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCCGCGCCttt                                                                                                                                          <  1:1377986/150‑1 (MQ=255)
               aGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCttttt                                                                                                                                        >  1:1225985/1‑150 (MQ=255)
               aGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCtttta                                                                                                                                        <  1:237065/150‑2 (MQ=255)
                gCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGAATCCTGCGCCttttac                                                                                                                                       <  2:1135917/150‑3 (MQ=255)
                   ggaCTGGCGGCGCAAACCGTCTAGCGCAAAGTAGGACGTAACGATCCTTCCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTTAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtata                                                                                                                                    <  2:1307641/148‑1 (MQ=255)
                   gCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtata                                                                                                                                    <  1:837277/150‑1 (MQ=255)
                   gCACGGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTtata                                                                                                                                    <  1:1030746/150‑1 (MQ=255)
                    cACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                                                                   >  1:1397116/1‑150 (MQ=255)
                    cACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAg                                                                                                                                   <  2:271471/150‑1 (MQ=255)
                     aCTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAgg                                                                                                                                  <  1:1037040/150‑1 (MQ=255)
                     aCTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAgg                                                                                                                                  <  2:1748284/150‑1 (MQ=255)
                        ggcggcGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTAGAGGttt                                                                                                                               >  2:732273/1‑150 (MQ=255)
                        ggcggcGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTGCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGCTTCTGCGCCTTTTTTATAGGttt                                                                                                                               >  1:124003/1‑150 (MQ=255)
                              gcAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACa                                                                                                                         <  1:1963240/150‑1 (MQ=255)
                               cAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACaa                                                                                                                        >  1:1113474/1‑150 (MQ=255)
                                aaaCCGGAGAGAGCAAAGTAGGACGTAACGATACTTGCCCGTTCGGTTCTGCTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAAt                                                                                                                       <  1:252552/150‑1 (MQ=255)
                                 aaCCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATg                                                                                                                      <  2:1733740/150‑1 (MQ=255)
                                   ccGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTGAAGACAATGaa                                                                                                                    >  1:1251019/1‑150 (MQ=255)
                                     ggAGCGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCATTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGACTGTTTTATAGGTTTAAGACAATGaaaa                                                                                                                  >  2:1054206/1‑150 (MQ=255)
                                     ggAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGaaaa                                                                                                                  >  2:1954728/1‑150 (MQ=255)
                                         agCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCt                                                                                                              <  2:298908/150‑1 (MQ=255)
                                           cgcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGc                                                                                                            <  2:1876756/150‑1 (MQ=255)
                                            gcAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCa                                                                                                           <  1:412241/150‑1 (MQ=255)
                                                   aGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                                    <  1:1501704/150‑1 (MQ=255)
                                                   aGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCg                                                                                                    >  2:652404/1‑150 (MQ=255)
                                                     gACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggt                                                                                                  >  2:469714/1‑150 (MQ=255)
                                                       cGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtag                                                                                                >  1:1656825/1‑150 (MQ=255)
                                                        gTAACGATGCTTGCCCGTGCGGTTCTGGTAAAAAATTCACGCAGTGACATGGCCTCCTGGTTTAATAGCTAACTGTTGAAGTAAAAGGGGCATTTTTCTGCTCCTTTTTTATAGGTTTAAGACCATGTTAAAGCGGCTAATTGCggtagg                                                                                               <  2:137863/150‑1 (MQ=255)
                                                         tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                              >  1:1303494/1‑150 (MQ=255)
                                                         tAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCggtaggt                                                                                              >  2:930892/1‑150 (MQ=255)
                                                              aTCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGtattat                                                                                         <  2:954873/150‑1 (MQ=255)
                                                                ccTTGCCCGTGCGGTTCTGGTAAAAAATACAATCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTAGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGTTattatta                                                                                       >  1:1517529/1‑149 (MQ=255)
                                                                ccTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTc                                                                                       <  2:1653814/150‑1 (MQ=255)
                                                                 cTTGCCCGTGCTGTTCTGGTAAAAAATACAAGAAGTGCCATGGCTGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCg                                                                                      <  1:365253/150‑1 (MQ=255)
                                                                 cTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAATCTGCAAATTGCGGTAGGTATTATTCg                                                                                      >  2:981251/1‑150 (MQ=255)
                                                                 cTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCg                                                                                      <  2:102648/150‑1 (MQ=255)
                                                                 cTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCg                                                                                      <  2:783930/150‑1 (MQ=255)
                                                                     cccGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAAc                                                                                  >  2:753191/1‑150 (MQ=255)
                                                                       cGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACga                                                                                <  2:1751346/150‑1 (MQ=255)
                                                                        gTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACgag                                                                               <  2:1578148/150‑1 (MQ=255)
                                                                                 tGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaa                                                                      >  2:335621/1‑150 (MQ=255)
                                                                                  ggTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGaaa                                                                     >  2:902681/1‑150 (MQ=255)
                                                                                   gTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAAt                                                                    >  1:1049558/1‑150 (MQ=255)
                                                                                       aaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                                <  1:1424041/150‑1 (MQ=255)
                                                                                       aaaaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCttt                                                                <  2:1572591/150‑1 (MQ=255)
                                                                                        aaaTACTAGCAGTGCCTTGGCGGCCTGCAATAAAATGTAACTGATGAAGTAATAGGCGCAGGATTCTGTGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCATCGAGAACAATGAAATCTTta                                                               <  1:344758/150‑1 (MQ=255)
                                                                                        aaaTACAAGCAGGGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTta                                                               <  1:464281/150‑1 (MQ=255)
                                                                                         aaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtat                                                              >  1:701004/1‑150 (MQ=255)
                                                                                         aaTACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTtat                                                              >  2:1007429/1‑150 (MQ=255)
                                                                                            aCAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAc                                                           <  2:1303494/150‑1 (MQ=255)
                                                                                             cAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcg                                                          >  2:1155279/1‑150 (MQ=255)
                                                                                             cAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcg                                                          >  1:1497207/1‑150 (MQ=255)
                                                                                              aaGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcgc                                                         <  2:1119146/150‑1 (MQ=255)
                                                                                               aGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcgcg                                                        >  2:342004/1‑150 (MQ=255)
                                                                                               aGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAAcgcg                                                        <  2:219921/150‑1 (MQ=255)
                                                                                               aGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTAGTCGCAACGAGAACAATGAAATCTTTATAAcgcg                                                        >  2:22363/1‑150 (MQ=255)
                                                                                                   gTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTcgc                                                    >  2:895159/1‑150 (MQ=255)
                                                                                                     gCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGAGTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCGTTATAACGCGTCgcgc                                                  >  1:221912/1‑150 (MQ=255)
                                                                                                       cATGGCCGCCTGAAATAAAAGCTAACTGTTTAAGTAAAAGGCTCAGGATTCAGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgcag                                                <  1:490651/150‑1 (MQ=255)
                                                                                                        aTGTCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCgcagc                                               <  1:1336887/150‑1 (MQ=255)
                                                                                                           gccgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGa                                            >  2:637406/1‑150 (MQ=255)
                                                                                                            ccgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGAt                                           >  1:479905/1‑150 (MQ=255)
                                                                                                             cgccTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATg                                          <  2:928618/150‑1 (MQ=255)
                                                                                                                 tGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGca                                      <  1:98529/150‑1 (MQ=255)
                                                                                                                  gCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGATCAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGcac                                     <  2:124003/150‑1 (MQ=255)
                                                                                                                  gCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGAGTCGCGCAGCAGATGCGcac                                     >  2:1851550/1‑150 (MQ=255)
                                                                                                                  gCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTACAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATgcgccc                                     >  2:1225585/1‑148 (MQ=255)
                                                                                                                     aTAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATg                                  <  1:573693/150‑1 (MQ=255)
                                                                                                                       aaaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGc                                >  2:589408/1‑150 (MQ=255)
                                                                                                                        aaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCg                               >  1:1916688/1‑150 (MQ=255)
                                                                                                                        aaaGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTAAAACGCGTCTCGCAGCAGATGCGCACATggag                               >  2:1414253/1‑148 (MQ=255)
                                                                                                                            cTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATa                           <  1:205533/150‑1 (MQ=255)
                                                                                                                            cTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATa                           <  1:198785/150‑1 (MQ=255)
                                                                                                                             tAACTGTTGAAGTAAAAGCCGCACGATTCTGCGCCTTTTTTATAGGATTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTAATAACGCGTCGCGCAGCAAATGCGCACATGACGAATaa                          >  2:470947/1‑150 (MQ=255)
                                                                                                                                   ttGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  1:512178/150‑1 (MQ=255)
                                                                                                                                   ttGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGa                    <  1:652404/150‑1 (MQ=255)
                                                                                                                                   ttGAAGTAAAAGGCGCAGGATTCTGCGCCGTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATAGCGGTAGGTATTTTTCGCAACGATAACAATGAATTATTTATAACTCGTCTCGCAGCAGATTCGCACATGGCGAATTATCTGGa                    >  2:978035/1‑150 (MQ=255)
                                                                                                                                    tGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAg                   <  1:930892/150‑1 (MQ=255)
                                                                                                                                     gAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGGAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACAGGAGt                  >  1:1263633/1‑150 (MQ=255)
                                                                                                                                      ataaaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGtt                 >  2:1684186/5‑150 (MQ=255)
                                                                                                                                       aGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGttt                >  2:694249/1‑150 (MQ=255)
                                                                                                                                        gTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTc               >  2:774765/1‑150 (MQ=255)
                                                                                                                                        aaaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGttta               >  2:1631381/3‑149 (MQ=255)
                                                                                                                                         aaaaaGGCGCAGAATCCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGATCAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTcc              >  1:171248/2‑150 (MQ=255)
                                                                                                                                          aaaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTccc             >  2:1707895/1‑150 (MQ=255)
                                                                                                                                          aaaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTccc             >  2:1758849/1‑150 (MQ=255)
                                                                                                                                           aaaGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCTCAACGAGAACAATGAAATCTTTATAACGCGTCGCGAATCAGATGCGCACATGGCGAATAAACAGGAGTTTCCcg            >  1:576616/1‑150 (MQ=255)
                                                                                                                                             aGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCcggc          <  2:1523027/150‑1 (MQ=255)
                                                                                                                                               gcgaAGGATTATTCGCATTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCTTTAGGTATTATTCTCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGAACATGGAGAATAAACTGGAGTTTCCcggcgg        <  1:787877/150‑1 (MQ=255)
                                                                                                                                                    ggATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTaaaa   >  1:1227412/1‑150 (MQ=255)
                                                                                                                                                    ggATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGGaaaa   >  2:1848336/1‑150 (MQ=255)
                                                                                                                                                     gATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAAt  <  2:235203/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACTCTGCAGCCGCAGCTGCACTGGCGGCGCAAACCGGAGAGCGCAAAGTAGGACGTAACGATCCTTGCCCGTGCGGTTCTGGTAAAAAATACAAGCAGTGCCATGGCCGCCTGCAATAAAAGCTAACTGTTGAAGTAAAAGGCGCAGGATTCTGCGCCTTTTTTATAGGTTTAAGACAATGAAAAAGCTGCAAATTGCGGTAGGTATTATTCGCAACGAGAACAATGAAATCTTTATAACGCGTCGCGCAGCAGATGCGCACATGGCGAATAAACTGGAGTTTCCCGGCGGTAAAAT  >  NC_000913/110864‑111162

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