Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,311,616 C→T 6.3% intergenic (‑82/+30) 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,6160CT 6.3% 249.1 / 5.5 79intergenic (‑82/+30)apbE/ompCputative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base T (2/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

GATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGAC  >  NC_000913/2311471‑2311757
                                                                                                                                                 |                                                                                                                                             
gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCa                                                                                                                                            >  4:264843/1‑149 (MQ=255)
gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCa                                                                                                                                            >  4:55118/1‑149 (MQ=255)
gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCa                                                                                                                                            >  1:365689/1‑149 (MQ=255)
    gCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGcc                                                                                                                                        <  2:61554/149‑1 (MQ=255)
    gCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGcc                                                                                                                                        <  4:566922/149‑1 (MQ=255)
      aaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCt                                                                                                                                      >  3:261859/1‑149 (MQ=255)
      aaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCt                                                                                                                                      >  3:353231/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                                                          >  4:53959/1‑149 (MQ=255)
                                                                                  gAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGt                                                          >  2:493708/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                                           <  3:136790/149‑1 (MQ=255)
                                                                                                 tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGt                                           >  1:448198/1‑149 (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)
                                                                                                                                                 |                                                                                                                                             
GATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGAC  >  NC_000913/2311471‑2311757

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