Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,311,616 C→T 9.6% 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 9.6% 270.7 / 12.4 94intergenic (‑82/+30)apbE/ompCputative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (50/34);  new base T (5/4);  total (55/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTT  >  NC_000913/2311469‑2311763
                                                                                                                                                   |                                                                                                                                                   
ttGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTAAACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGccc                                                                                                                                                     <  2:301699/148‑1 (MQ=255)
  gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCg                                                                                                                                                    >  2:275583/1‑147 (MQ=255)
    tCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAgg                                                                                                                                                >  7:282261/1‑149 (MQ=255)
        aaccaacAAAGAAGAGCGCTGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAACAAAGGGCCTGCGGGCCCt                                                                                                                                            <  5:219794/149‑1 (MQ=255)
        aaccaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCt                                                                                                                                            >  4:56365/1‑149 (MQ=255)
         accaacAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAACAAAGGGCCTGCGGGCCCtt                                                                                                                                           >  2:146347/1‑149 (MQ=255)
           caacaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGCTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTg                                                                                                                                         >  8:238354/1‑149 (MQ=255)
            aacaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGTTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCACTTTGt                                                                                                                                        >  6:128482/1‑149 (MQ=255)
            aacaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGt                                                                                                                                        >  3:155764/1‑149 (MQ=255)
               aaaGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGGGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTc                                                                                                                                      >  7:365008/1‑148 (MQ=255)
                 agaagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTGATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGa                                                                                                                                    >  5:295293/1‑148 (MQ=255)
                  gaagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGata                                                                                                                                  >  7:75114/1‑149 (MQ=255)
                   aagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGatat                                                                                                                                 <  8:282261/149‑1 (MQ=255)
                    agagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATc                                                                                                                                <  3:56365/149‑1 (MQ=255)
                     gagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATc                                                                                                                                >  3:323461/1‑148 (MQ=255)
                       gcgcGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATGTCCATACGGGAGTATCCGCATGCGACCGGGCGCCAGAGTACCCCAAAAACCTGCGGGTGTACGCTGAAACCAATGCAAAAAGGGCCCGCCGGCCCTGGGTGCGAGATCact                                                                                                                             >  1:45706/1‑147 (MQ=255)
                                     agTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTg                                                                                                               <  8:75114/149‑1 (MQ=255)
                                      gTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTgg                                                                                                              <  4:184888/149‑1 (MQ=255)
                                         ccACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCt                                                                                                                                            <  6:368862/116‑1 (MQ=255)
                                          cACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTaaa                                                                                                          >  8:222276/1‑149 (MQ=255)
                                             gcgGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCa                                                                                                       <  7:270219/149‑1 (MQ=255)
                                              cgGGTAAAGCTTATTTCCATACGTGATTAACCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCAGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAg                                                                                                      >  3:37459/1‑149 (MQ=255)
                                                   aaaGCTTATTTCCATACGTGATTATCCTCATGCTAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAAc                                                                                                                 <  3:292175/133‑1 (MQ=255)
                                                   aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCTTTTTTcattg                                                                                                                                  <  5:145394/116‑5 (MQ=39)
                                                   aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCTTTTTTcattg                                                                                                                                  >  6:145394/1‑112 (MQ=39)
                                                   aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCTTTTTTcattg                                                                                                                                  <  7:311907/116‑5 (MQ=39)
                                                   aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCTTTTTTcattg                                                                                                                                  >  8:311907/1‑112 (MQ=39)
                                                   aaaGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAAc                                                                                                                 >  4:292175/1‑133 (MQ=255)
                                                     aGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCaga                                                                                               <  4:155764/149‑1 (MQ=255)
                                                      gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCGAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGat                                                                                                                                   >  5:54144/1‑112 (MQ=255)
                                                      gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCGAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGat                                                                                                                                   <  6:54144/112‑1 (MQ=255)
                                                        ttATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGc                                                                                             >  8:348645/1‑148 (MQ=255)
                                                               cATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTa                                                                                           <  8:47362/143‑1 (MQ=255)
                                                               cATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTa                                                                                           >  7:47362/1‑143 (MQ=255)
                                                               cATACGTGATTAACCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCCGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGt                                                                                     <  4:37459/149‑1 (MQ=255)
                                                                 tACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTa                                                                                   >  2:301921/1‑149 (MQ=255)
                                                                    gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCa                                                                                <  3:329693/149‑1 (MQ=255)
                                                                      gATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgt                                                                              >  5:26580/1‑149 (MQ=255)
                                                                        ttATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgtg                                                                             <  7:222276/148‑1 (MQ=255)
                                                                          aTCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGAGtt                                                                           >  5:326167/1‑148 (MQ=255)
                                                                             cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgtg                                                                             >  7:10144/1‑143 (MQ=255)
                                                                             cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgtg                                                                             <  8:10144/143‑1 (MQ=255)
                                                                             cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                                        <  7:348645/148‑1 (MQ=255)
                                                                             cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATg                                                                       >  2:204358/1‑149 (MQ=255)
                                                                              tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                      >  5:92498/1‑149 (MQ=255)
                                                                              tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACACAGAGCTACGATGTTATAAGTGTTGATGc                                                                      >  1:67697/1‑149 (MQ=255)
                                                                              tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTCTACGCTGAAAACAATGAAAAAAGGGACCGAAGGCCCTTAGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTCTTGATGc                                                                      >  4:169576/1‑149 (MQ=255)
                                                                                aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                           >  6:162381/1‑142 (MQ=255)
                                                                                aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                           <  5:162381/142‑1 (MQ=255)
                                                                                  gCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGc                                                                  >  5:15261/1‑149 (MQ=255)
                                                                                    gAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGCGCTACGATGTTATCAGTGTTGATGCCAGCGt                                                                >  6:109011/1‑149 (MQ=255)
                                                                                         gTCGCAAGTGTCCCCCAAAATCCGGAGGTTGTACGCTGAAAACAATGCAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCACTCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGCTGTTCGCAGTGTTGATGCCCGCGTCCCGa                                                           <  2:45706/149‑1 (MQ=255)
                                                                                              aaGAGTACCCCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGaa                                                      >  3:135686/1‑149 (MQ=255)
                                                                                               agagTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGTGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAAc                                                     <  1:275583/149‑1 (MQ=255)
                                                                                                  gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGAt                                                                                       >  7:157666/1‑112 (MQ=255)
                                                                                                  gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGAt                                                                                       <  8:157666/112‑1 (MQ=255)
                                                                                                  gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                  <  4:323461/149‑1 (MQ=255)
                                                                                                      acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                              >  6:135747/1‑149 (MQ=255)
                                                                                                      acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                              >  7:297792/1‑149 (MQ=255)
                                                                                                      acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                              >  8:232732/1‑149 (MQ=255)
                                                                                                        cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCTTTTTTcat                                                                                                                                    >  8:182906/1‑59 (MQ=25)
                                                                                                        cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCTTTTTTcat                                                                                                                                    <  7:182906/61‑3 (MQ=25)
                                                                                                        cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcg                                            >  4:240393/1‑149 (MQ=255)
                                                                                                        cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCACTTTGTTCGATATCAATCGAGATTAGCACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGGGTTGATGCCCGCGTCACCAGTGAACTTGTTgtcg                                            >  6:213477/1‑149 (MQ=255)
                                                                                                         aaaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgt                                           >  3:261352/1‑149 (MQ=255)
                                                                                                          aaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                           <  8:259511/116‑1 (MQ=255)
                                                                                                          aaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGtt                                                                           >  7:259511/1‑116 (MQ=255)
                                                                                                          aaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGAAGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgacgtc                                          >  7:284678/1‑149 (MQ=255)
                                                                                                             aaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcag                                       >  5:146239/1‑149 (MQ=255)
                                                                                                              aCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagc                                      <  4:261352/149‑1 (MQ=255)
                                                                                                                      tGTACGCTTAAAACAATTAAAAAAGGTCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTTCGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGa                               <  8:284678/148‑1 (MQ=255)
                                                                                                                        tACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAtt                             <  1:301921/148‑1 (MQ=255)
                                                                                                                         aCGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgtc                                          <  3:352672/134‑1 (MQ=255)
                                                                                                                         aCGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcgtc                                          >  4:352672/1‑134 (MQ=255)
                                                                                                                          cGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGAGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTg                          <  6:326167/149‑1 (MQ=255)
                                                                                                                            cTGAAAACAATGAATAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCGAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                        <  5:213477/149‑1 (MQ=255)
                                                                                                                            cTGAAAACAATGAACACAGGCCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGACCTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                        <  2:67697/149‑1 (MQ=255)
                                                                                                                             tGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACAGGTTGTCGTCCAGCAGGTTGATTTTg                          >  7:261220/1‑146 (MQ=255)
                                                                                                                               aaaaCATTGAAAACAGGGCCCGCTGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAg                       <  6:15261/147‑1 (MQ=255)
                                                                                                                               aaaaCAATGAAAAACGTGCCCGCTGTCCCTTTGTTCGATATCTATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTa                        <  8:261220/146‑1 (MQ=255)
                                                                                                                                aaaCAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTc                     >  2:36826/1‑148 (MQ=255)
                                                                                                                                    aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                                        >  6:220961/1‑93 (MQ=255)
                                                                                                                                    aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                                        <  5:220961/93‑1 (MQ=255)
                                                                                                                                    aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACg                 >  1:380759/1‑148 (MQ=255)
                                                                                                                                      tGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAAc                                                                                                                 >  7:386967/1‑50 (MQ=255)
                                                                                                                                      tGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAAc                                                                                                                 <  8:386967/50‑1 (MQ=255)
                                                                                                                                        aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAgg             <  5:135747/148‑1 (MQ=255)
                                                                                                                                        aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGt            >  2:92157/1‑149 (MQ=255)
                                                                                                                                         aaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg           >  1:37887/1‑149 (MQ=255)
                                                                                                                                          aaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg           >  4:317100/1‑148 (MQ=255)
                                                                                                                                           aaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa         >  4:145370/1‑149 (MQ=255)
                                                                                                                                           aaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa         >  8:35924/1‑149 (MQ=255)
                                                                                                                                            aaGGGCCTGCGGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa         >  4:258382/1‑148 (MQ=255)
                                                                                                                                              gggCCAGCAGGCCCGTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACAt      <  6:295293/149‑1 (MQ=255)
                                                                                                                                                gCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGt    <  3:240393/149‑1 (MQ=255)
                                                                                                                                                  ccgccgGCCCTTTGTTCGATAGCAAGCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGCTGCCAGCGTCACGCGTGAACTGGTTGTCGTCCAGCAGGTGGATTTTGTAGTCAACGTAGGTGGACATGttt  <  4:135686/149‑1 (MQ=255)
                                                                                                                                                  ccGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGa                                                           <  3:206978/92‑1 (MQ=255)
                                                                                                                                                  ccGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGa                                                           >  4:206978/1‑92 (MQ=255)
                                                                                                                                                  ccGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttt  >  3:92373/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTT  >  NC_000913/2311469‑2311763

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: ‑