Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,587,825 C→A 9.6% A158S (GCG→TCG)  ydeR ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,587,8250CA 9.6% 267.5 / 5.3 88A158S (GCG→TCG) ydeRputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (29/49);  new base A (7/3);  total (36/52)
Fisher's exact test for biased strand distribution p-value = 8.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATCCCTCAGCTCTA  >  NC_000913/1587679‑1587971
                                                                                                                                                  |                                                                                                                                                  
gccATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCATGTGTAGATCACATTGATTAGCGcc                                                                                                                                                   <  4:51746/148‑1 (MQ=255)
    ttAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGa                                                                                                                                               <  1:147036/148‑1 (MQ=255)
    ttAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGa                                                                                                                                               >  7:380877/1‑148 (MQ=255)
    ttAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGa                                                                                                                                               >  1:43292/1‑148 (MQ=255)
     tAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGAt                                                                                                                                              >  5:92814/1‑148 (MQ=255)
      aaTAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGACCACATTGATTAGCGCCTCGAt                                                                                                                                              >  8:38079/1‑147 (MQ=255)
      aaTAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGACCACATTGATTAGCGCCTCGAt                                                                                                                                              <  7:38079/147‑1 (MQ=255)
          gaagaTAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTc                                                                                                                                         >  5:76085/1‑148 (MQ=255)
          gaagaTAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTc                                                                                                                                         <  8:115008/148‑1 (MQ=255)
          gaagaTAAATACCGAACAGGACCTTTATAATGATTGTTTTACCCATTATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGTTGTAGATCACATTGATTAGCGACTCGATCGtt                                                                                                                                          <  2:197793/147‑1 (MQ=255)
          gaagaTAAATACAGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTc                                                                                                                                         <  1:157328/148‑1 (MQ=255)
           aataTAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTcc                                                                                                                                        <  4:310941/145‑1 (MQ=255)
           aagaTAAATACCGAACAGGACCTTTATACGGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTcc                                                                                                                                        >  7:145224/1‑148 (MQ=255)
             gaTAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGACTCTATcataaact                                                                                                                                      >  8:171058/1‑141 (MQ=255)
              aTAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGACTCAATCGATCACt                                                                                                                                      >  2:309732/1‑147 (MQ=255)
                        aaCAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTTATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGttt                                                                                                                           <  3:347443/148‑1 (MQ=255)
                            ggACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCAt                                                                                                                       <  2:115178/148‑1 (MQ=255)
                            ggACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCAt                                                                                                                       <  7:348013/148‑1 (MQ=255)
                                  ttATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGAGTTTCCATTCACCg                                                                                                                 >  7:316273/1‑148 (MQ=255)
                                   tataCTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTGCCCGt                                                                                                                >  1:155119/1‑148 (MQ=255)
                                     taCTGATTGTTTTACCCATGCTATCTCCTAAGGGTTAATATTGATTTAAATTGTAGGCTTAGAAAATCGAATCGACGTTATTGCCAGGTGGAGATCACATTGATTAGCGAGTCGATGGTTCCCTGGCTTGGGTTTCCATTCACCGTGa                                                                                                              <  6:124016/148‑1 (MQ=255)
                                     taCTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTTAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACTTTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTg                                                                                                               <  1:111836/147‑1 (MQ=255)
                                     taCTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTg                                                                                                               <  7:155459/147‑1 (MQ=255)
                                     taCTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGa                                                                                                              <  8:211268/148‑1 (MQ=255)
                                     taCTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGa                                                                                                              <  3:81331/148‑1 (MQ=255)
                                      catGATTGTTTTCCCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGATGTCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGa                                                                                                              <  6:118718/145‑1 (MQ=255)
                                         gATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTcc                                                                                                                         <  4:171874/133‑1 (MQ=255)
                                         gATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTcc                                                                                                                         >  3:171874/1‑133 (MQ=255)
                                         gATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGc                                                                                                          >  3:207765/1‑148 (MQ=255)
                                          aTTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATa                                                                                                            <  7:179071/145‑1 (MQ=255)
                                             gTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGctct                                                                                                       <  1:309732/147‑1 (MQ=255)
                                              ttttACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTg                                                                                                     <  3:279923/148‑1 (MQ=255)
                                               tttACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAATTCGAATAGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGc                                                                                                    >  7:419911/1‑148 (MQ=255)
                                                    ccATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTaa                                                                                               <  6:353547/148‑1 (MQ=255)
                                                       tGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCAc                                                                                                                   <  3:32268/125‑1 (MQ=255)
                                                       tGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTg                                                                                            <  1:48717/148‑1 (MQ=255)
                                                       tGATATATACTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTg                                                                                            <  8:145224/148‑1 (MQ=255)
                                                        gATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTg                                                                                            >  7:58665/1‑147 (MQ=255)
                                                        gATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATAACCGAATCAATATATAAATACCTTGCGTTTCCATTCAcc                                                                                                                  >  4:32268/1‑125 (MQ=255)
                                                         atatatCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTcc                                                                                                                         <  5:74529/117‑1 (MQ=255)
                                                           atttCCTAAGGTTAAAATTTGATTTAAATAGAAGGCTAAGAAAATCGATTGGAGGTTATTGCCAGGTGTAGATAACATTGATTAGATCCTCGATCGTTCGCTTGCTTGCGGTTCCATTAACCCTGATAGCTCTTGCCTTAAGTGGaa                                                                                         <  2:72778/144‑1 (MQ=255)
                                                           atatCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGAATCCATATATCACTGTCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGaaa                                                                                        >  2:368589/1‑148 (MQ=255)
                                                            tatCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGAGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGaaa                                                                                        <  8:316273/147‑1 (MQ=255)
                                                              tCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAGCtg                                                                                     <  4:207765/148‑1 (MQ=255)
                                                              tCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGACAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACtg                                                                                     >  1:134361/1‑148 (MQ=255)
                                                                 tAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTgg                                                                                           >  4:258118/1‑139 (MQ=255)
                                                                 tAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTgg                                                                                           <  3:258118/139‑1 (MQ=255)
                                                                 tAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGaa                                                                                         <  8:400664/141‑1 (MQ=255)
                                                                 tAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGaa                                                                                         >  7:400664/1‑141 (MQ=255)
                                                                   atgTTAAAATTTGATTTAAAAAGATGGCTAAGTAAATAGATTAGACGTTATTGCCAGTTTTATATAAATTTGATTTGCGACTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTTCCTTATGTGGAAACTGTGc                                                                                  <  8:104737/144‑1 (MQ=255)
                                                                     gTTAAAAATTGATTTAAAAAGTAGTCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCtt                                                                                                      <  5:19141/124‑1 (MQ=255)
                                                                      ttAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATAAGCGACTCAATCTTTACCTGACTTGCGTTTCCATTCACCGTGATAGCTCTTg                                                                                                     >  6:19141/1‑124 (MQ=255)
                                                                       tAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTc                                                                                                                                         <  8:138309/87‑1 (MQ=255)
                                                                       tAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATGAGCGCCTCGATCGTTc                                                                                                                                         >  7:138309/1‑87 (MQ=255)
                                                                        aaaaaTTGATTTAAAAAGAAGGCTAAGTAAATCGAATAGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACg                                                                            <  7:332556/147‑1 (MQ=255)
                                                                         aaaaTTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGa                                                                           <  5:319453/147‑1 (MQ=255)
                                                                           aaTTGATTTAAAAAGAAGGCTAAGAAACTCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTg                                                                        <  1:306942/148‑1 (MQ=255)
                                                                            aTTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCTTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGa                                                                       <  1:368589/148‑1 (MQ=255)
                                                                            aTTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGa                                                                       >  5:292625/1‑148 (MQ=255)
                                                                            aTTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGa                                                                       <  2:198641/148‑1 (MQ=255)
                                                                             ttGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATGACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGAt                                                                      >  1:363813/1‑148 (MQ=255)
                                                                                 tttAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGAtctc                                                                   <  8:58665/147‑1 (MQ=255)
                                                                                 tttAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCtca                                                                  <  6:292625/148‑1 (MQ=255)
                                                                                   tAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCtcat                                                                 <  6:92814/147‑1 (MQ=255)
                                                                                   tAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGAAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCtcat                                                                 <  4:64382/147‑1 (MQ=255)
                                                                                      aaaGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATAACCGAATCAATAAAAAAAAAACTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGGGATCTCATcaa                                                              >  4:191850/1‑147 (MQ=255)
                                                                                                 aGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGttt                                                  <  6:76085/148‑1 (MQ=255)
                                                                                                 aGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGttt                                                  <  2:152441/148‑1 (MQ=255)
                                                                                                 aGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGtt                                                   <  2:43292/147‑1 (MQ=255)
                                                                                                  gAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGtttt                                                 <  8:380877/148‑1 (MQ=255)
                                                                                                     aatcgaatcgACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCt                                              <  8:88499/148‑1 (MQ=255)
                                                                                                      atcgaatcgaCGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTa                                             <  7:171058/148‑1 (MQ=255)
                                                                                                           atcgaCGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGc                                         >  7:277892/1‑147 (MQ=255)
                                                                                                            tcgaCGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCa                                       <  7:251201/148‑1 (MQ=255)
                                                                                                                         cAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACgcag                          >  3:124672/1‑148 (MQ=255)
                                                                                                                           ggTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTTATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACgcagca                        >  1:374389/1‑148 (MQ=255)
                                                                                                                            gtgtAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGa                                                                                                              <  7:236931/61‑1 (MQ=255)
                                                                                                                             tgtAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTCGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATc                      >  1:377811/1‑148 (MQ=255)
                                                                                                                             tgtAGATCACATTGATTACCGACTCAATAAATAAATGTCTTGCGTTTCCATTCACCGTGAt                                                                                                             >  8:236931/1‑61 (MQ=255)
                                                                                                                              gtAGATCACATTGATTAGCGCCTCGATCGATCCCTGGCTTGCGTTTCCATTCACCGTGATAGATCTAGCCTTAAGTGGAAACTGTGCATTCCGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACTCAGCATcc                     >  1:202672/1‑148 (MQ=255)
                                                                                                                               tataTCACATTGATTAGCTCATAGATCGTTCCATGGCTTGCGTTTCCATACACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATcc                     <  4:122719/144‑1 (MQ=255)
                                                                                                                                   tCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGt                 >  2:190687/1‑147 (MQ=255)
                                                                                                                                   tCACATTGATTAGCGCCTCCATCGTTCCCTGGCTTGCGTTTCCTTTCACCGTGATAGCTCTTGCCTTAAGTGGAAACAGTGCATTACGTGTGATAACATCAACAATACCCGTTTTCCTATCGCCATTTTTTAACGCATCATCCTTGt                 >  3:276912/1‑147 (MQ=255)
                                                                                                                                   tCACATTGATTAGCGCATCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTc                <  8:419911/148‑1 (MQ=255)
                                                                                                                                      caTTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATc             >  1:143587/1‑148 (MQ=255)
                                                                                                                                        ttGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATACCCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATccc           >  2:319579/1‑148 (MQ=255)
                                                                                                                                        ttGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATccc           >  7:348651/1‑148 (MQ=255)
                                                                                                                                           aTTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTTCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATCCCTc         <  3:71382/147‑1 (MQ=255)
                                                                                                                                             tAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATCCCTCAGc      >  1:72109/1‑148 (MQ=255)
                                                                                                                                               gcgcCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTTATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATCCCTCAGct     <  2:374389/147‑1 (MQ=255)
                                                                                                                                                 gcCTCGATCGTTCCCTGGCTTGCTTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATCCCTCAGCTCTa  <  4:307048/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATTGATTAGCGCCTCGATCGTTCCCTGGCTTGCGTTTCCATTCACCGTGATAGCTCTTGCCTTAAGTGGAAACTGTGCATTACGAGTGATCTCATCAACAATAACCGTTTTGCTATCGCCATTTTTTAACGCAGCATCCTGGTCATCCCTCAGCTCTA  >  NC_000913/1587679‑1587971

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