Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,311,619 A→G 6.3% intergenic (‑85/+27) apbE ← / ← ompC putative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,311,6190AG 6.3% 247.1 / 6.6 79intergenic (‑85/+27)apbE/ompCputative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base G (2/3);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

GATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTG  >  NC_000913/2311471‑2311766
                                                                                                                                                    |                                                                                                                                                   
gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCa                                                                                                                                                     >  4:55118/1‑149 (MQ=255)
gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCa                                                                                                                                                     >  4:264843/1‑149 (MQ=255)
gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCa                                                                                                                                                     >  1:365689/1‑149 (MQ=255)
    gCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGcc                                                                                                                                                 <  4:566922/149‑1 (MQ=255)
    gCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGcc                                                                                                                                                 <  2:61554/149‑1 (MQ=255)
      aaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCt                                                                                                                                               >  3:353231/1‑149 (MQ=255)
      aaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCt                                                                                                                                               >  3:261859/1‑149 (MQ=255)
        ccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCttt                                                                                                                                             >  4:208238/1‑149 (MQ=255)
          aacaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGt                                                                                                                                           >  2:300821/1‑149 (MQ=255)
           acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGtt                                                                                                                                          >  1:116850/1‑149 (MQ=255)
           acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGtt                                                                                                                                          >  3:560713/1‑149 (MQ=255)
           acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGtt                                                                                                                                          <  3:397859/149‑1 (MQ=255)
            caaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTc                                                                                                                                         >  3:317671/1‑149 (MQ=255)
             aaaGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCg                                                                                                                                        >  3:559054/1‑149 (MQ=255)
             aaaGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCg                                                                                                                                        >  4:352943/1‑149 (MQ=255)
              aagaagAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAACAAAGGGCCTGCGGGCCCTTTTTTcat                                                                                                                                       <  3:488407/149‑3 (MQ=255)
              aagaagAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGa                                                                                                                                       <  3:103199/149‑1 (MQ=255)
               agaagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCTTTTTTcatt                                                                                                                                      <  1:51724/149‑4 (MQ=255)
               agaagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGat                                                                                                                                      <  3:352943/149‑1 (MQ=255)
                   gagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCa                                                                                                                                  <  3:115897/149‑1 (MQ=255)
                    agCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCaa                                                                                                                                 <  1:443190/149‑1 (MQ=255)
                     gcgcGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAAt                                                                                                                                <  4:559054/149‑1 (MQ=255)
                      cgcgGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATc                                                                                                                               <  4:513356/149‑1 (MQ=255)
                      agcgGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATc                                                                                                                               <  2:276307/148‑1 (MQ=255)
                       gcgGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACACTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCg                                                                                                                              >  3:68313/1‑149 (MQ=255)
                             gccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTa                                                                                                                        >  1:219642/1‑149 (MQ=255)
                                 gcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAAc                                                                                                                    <  1:225973/149‑1 (MQ=255)
                                        cACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTATAACTGGTaaa                                                                                                             <  2:453850/149‑1 (MQ=255)
                                        cACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTaaa                                                                                                             <  2:489571/149‑1 (MQ=255)
                                         aCGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAAc                                                                                                            >  1:206343/1‑149 (MQ=255)
                                         aCGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAAc                                                                                                            <  4:500371/149‑1 (MQ=255)
                                                   aGCTTATTTCCATACGTGATTATCCCCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCaga                                                                                                  >  2:542964/1‑149 (MQ=255)
                                                     cTTATTTCCATTCGTGATTATCCTCATGCGAACGGTCGCACGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGc                                                                                                >  3:491184/1‑149 (MQ=255)
                                                     cTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGc                                                                                                >  4:140152/1‑149 (MQ=255)
                                                     cTTAGTTCCAGACGTTATTATACTCATGCGAACGGTCGAAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTTGTAAGCCAGACCCAGAGc                                                                                                <  1:329483/149‑1 (MQ=255)
                                                        aTTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTAc                                                                                             <  3:208238/149‑1 (MQ=255)
                                                                    gATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgt                                                                                 <  2:545072/149‑1 (MQ=255)
                                                                      ttATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgtgt                                                                               <  4:196831/149‑1 (MQ=255)
                                                                       tataCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                              <  4:448278/149‑1 (MQ=255)
                                                                       tATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                              <  3:385908/149‑1 (MQ=255)
                                                                        atcatcATGCAAACGGTAGCAAGAGTTCACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTg                                                                             <  3:410694/149‑1 (MQ=255)
                                                                        aTCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTg                                                                             >  3:91093/1‑149 (MQ=255)
                                                                            tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                         >  3:15/1‑149 (MQ=255)
                                                                               tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAg                                                                      >  2:382266/1‑149 (MQ=255)
                                                                               tGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAg                                                                      >  3:507687/1‑149 (MQ=255)
                                                                                  gAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGt                                                                   >  2:493708/1‑149 (MQ=255)
                                                                                  gAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGt                                                                   >  4:53959/1‑149 (MQ=255)
                                                                                       gTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGa                                                              <  3:53959/149‑1 (MQ=255)
                                                                                       gTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGAGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGa                                                              >  2:175062/1‑149 (MQ=255)
                                                                                         cGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGt                                                            <  3:55118/149‑1 (MQ=255)
                                                                                            aaGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGaa                                                         >  4:136790/1‑149 (MQ=255)
                                                                                            aaGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGaa                                                         <  1:175062/149‑1 (MQ=255)
                                                                                            aaGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAAc                                                        <  3:445092/149‑1 (MQ=255)
                                                                                                gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                     <  4:261859/149‑1 (MQ=255)
                                                                                                 tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGt                                                    >  1:448198/1‑149 (MQ=255)
                                                                                                 tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGt                                                    <  3:136790/149‑1 (MQ=255)
                                                                                                  acacCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                                   <  4:353231/149‑1 (MQ=255)
                                                                                                    acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                                 >  3:117715/1‑149 (MQ=255)
                                                                                                    acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                                 >  2:223775/1‑149 (MQ=255)
                                                                                                     ccAAAAAACTGGGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAATTGGTTgtc                                                <  4:5682/149‑1 (MQ=255)
                                                                                                       aaaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAACCCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgt                                              <  4:388211/149‑1 (MQ=255)
                                                                                                       aaaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgt                                              >  1:295491/1‑149 (MQ=255)
                                                                                                       aaaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgt                                              >  3:148928/1‑149 (MQ=255)
                                                                                                        aaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgtc                                             <  3:473037/149‑1 (MQ=255)
                                                                                                             cTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagca                                        >  2:561380/1‑149 (MQ=255)
                                                                                                                aGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGt                                     <  4:109536/149‑1 (MQ=255)
                                                                                                                       aCGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAtttt                              <  4:317671/149‑1 (MQ=255)
                                                                                                                        cGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTg                             >  3:67100/1‑149 (MQ=255)
                                                                                                                         gCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTGGTCCAGCAGGTTGATTTTGt                            >  3:465454/1‑149 (MQ=255)
                                                                                                                         gCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGt                            >  1:87996/1‑149 (MQ=255)
                                                                                                                          cTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                           <  2:365689/149‑1 (MQ=255)
                                                                                                                           tGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAg                          >  1:209863/1‑149 (MQ=255)
                                                                                                                                   aTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTa                  >  1:290655/1‑149 (MQ=255)
                                                                                                                                    tGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAg                 >  3:155307/1‑149 (MQ=255)
                                                                                                                                      aacaaaGGGCCTGCGGGCCCTTTTTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt               >  3:384516/4‑149 (MQ=255)
                                                                                                                                        aaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGGAGGTTGATTTTGTAGTCAACGTAGGTgg             <  4:15/149‑1 (MQ=255)
                                                                                                                                          aaGGGCCTGCGGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt               >  4:371937/1‑145 (MQ=255)
                                                                                                                                          aaGGGCCTGCGGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt               <  3:371937/145‑1 (MQ=255)
                                                                                                                                          aaGGGCCCGCAGGGCCTTTGTTCGATATCAATCGAGATTAGAACGGGTAAACCAGACCCAGAGTTGCGATGTTATCAGTGTTTATGCCAGCGTCAAGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGAc           <  4:335917/149‑1 (MQ=255)
                                                                                                                                                   cAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttg  <  4:507687/149‑1 (MQ=255)
                                                                                                                                                   cAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttg  >  3:438040/1‑149 (MQ=255)
                                                                                                                                                   cAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttg  <  4:67100/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTG  >  NC_000913/2311471‑2311766

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