Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,211,187 A→T 5.3% intergenic (+87/‑38) yehQ → / → yehR pseudogene/lipoprotein, DUF1307 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,211,1870AT 5.3% 272.9 / 4.3 75intergenic (+87/‑38)yehQ/yehRpseudogene/lipoprotein, DUF1307 family
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base T (2/2);  total (36/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 = 9.76e-01

ATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGACAAAGA  >  NC_000913/2211041‑2211304
                                                                                                                                                  |                                                                                                                     
aTGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGagag                                                                                                                      <  1:1119850/148‑1 (MQ=255)
        cGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGt                                                                                                              <  1:471844/148‑1 (MQ=255)
        cGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGt                                                                                                              <  3:791149/148‑1 (MQ=255)
         ggTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGtt                                                                                                             <  3:357951/148‑1 (MQ=255)
                          aCAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          >  3:306580/1‑134 (MQ=255)
                          aCAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          <  4:306580/134‑1 (MQ=255)
                          aCAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAAt                                                                                            >  3:909826/1‑148 (MQ=255)
                            aaTTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATgg                                                                                          <  4:909826/148‑1 (MQ=255)
                                  tCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGACTCTCACCAg                                                                                                               >  4:770543/1‑121 (MQ=255)
                                  tCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGACTCTCACCAg                                                                                                               <  3:770543/121‑1 (MQ=255)
                                     ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTcc                                                                                                           >  3:895036/1‑122 (MQ=255)
                                     ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTcc                                                                                                           <  4:895036/122‑1 (MQ=255)
                                     ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          >  1:339774/1‑123 (MQ=255)
                                     ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          <  2:339774/123‑1 (MQ=255)
                                     ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGtt                                                                                      >  1:107124/1‑143 (MQ=255)
                                     ccATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGtt                                                                                      <  2:107124/143‑1 (MQ=255)
                                       aTGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTc                                                                                                                    >  2:65373/1‑111 (MQ=255)
                                       aTGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTc                                                                                                                    <  1:1034644/111‑1 (MQ=255)
                                       aTGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTc                                                                                                                    >  2:1034644/1‑111 (MQ=255)
                                       aTGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTc                                                                                                                    <  1:65373/111‑1 (MQ=255)
                                                   caaGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAAt                                                                                            >  1:16794/1‑123 (MQ=255)
                                                   caaGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAAt                                                                                            <  2:16794/123‑1 (MQ=255)
                                                             aTAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGt                                                                                                              >  2:462883/1‑95 (MQ=255)
                                                             aTAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGt                                                                                                              <  1:462883/95‑1 (MQ=255)
                                                               aGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCa                                                                                                                   <  1:120741/88‑1 (MQ=255)
                                                               aGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCa                                                                                                                   >  2:120741/1‑88 (MQ=255)
                                                                  gCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGtttt                                                         >  3:238516/1‑143 (MQ=255)
                                                                  gCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGtttt                                                         <  4:238516/143‑1 (MQ=255)
                                                                     ttttACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGtttt                                                                                    >  3:751821/1‑113 (MQ=255)
                                                                     ttttACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGtttt                                                                                    <  4:751821/113‑1 (MQ=255)
                                                                            aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              >  2:1092202/1‑96 (MQ=255)
                                                                            aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              <  1:1092202/96‑1 (MQ=255)
                                                                            aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCtt                                                                        <  2:107456/118‑1 (MQ=255)
                                                                            aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCtt                                                                        >  1:107456/1‑118 (MQ=255)
                                                                            aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGACTCTCACCAg                                                                                                               >  4:488279/1‑79 (MQ=255)
                                                                            aaGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGACTCTCACCAg                                                                                                               <  3:488279/79‑1 (MQ=255)
                                                                              gagTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCa                                                                                                                   >  2:102635/1‑73 (MQ=255)
                                                                              gagTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCa                                                                                                                   <  1:102635/73‑1 (MQ=255)
                                                                                gTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGt                                      >  1:788283/1‑148 (MQ=255)
                                                                                 ttCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCgt                                                  <  2:1116136/135‑1 (MQ=255)
                                                                                 ttCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCgt                                                  >  1:1116136/1‑135 (MQ=255)
                                                                                  tCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              <  2:874089/90‑1 (MQ=255)
                                                                                  tCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              >  1:874089/1‑90 (MQ=255)
                                                                                  tCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTg                                                             <  1:725403/123‑1 (MQ=255)
                                                                                  tCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTg                                                             >  2:725403/1‑123 (MQ=255)
                                                                                       tCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          <  4:322568/73‑1 (MQ=255)
                                                                                       tCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCt                                                                                                          >  3:322568/1‑73 (MQ=255)
                                                                                          aGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGaa                                                                             <  1:446819/99‑1 (MQ=255)
                                                                                          aGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGaa                                                                             >  2:446819/1‑99 (MQ=255)
                                                                                           gTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGtttt                             <  1:4434/146‑1 (MQ=255)
                                                                                           gTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGtttt                             >  2:4434/1‑146 (MQ=255)
                                                                                             tCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATgg                                                                                          >  1:92641/1‑83 (MQ=255)
                                                                                             tCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATgg                                                                                          <  2:92641/83‑1 (MQ=255)
                                                                                                 gACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCtt                      <  4:248633/147‑1 (MQ=255)
                                                                                                 gACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCtt                      >  3:248633/1‑147 (MQ=255)
                                                                                                       tCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              <  1:84440/69‑1 (MQ=255)
                                                                                                       tCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATtaaataaa                                                                                              >  2:84440/1‑69 (MQ=255)
                                                                                                              aCGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGa        <  2:788283/148‑1 (MQ=255)
                                                                                                                     tAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCtt                                                                        >  3:417702/1‑77 (MQ=255)
                                                                                                                     tAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCtt                                                                        <  4:417702/77‑1 (MQ=255)
                                                                                                                       ggCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGt                                      >  3:147625/1‑109 (MQ=255)
                                                                                                                       ggCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGt                                      <  4:147625/109‑1 (MQ=255)
                                                                                                                        gCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgc                    <  1:218917/126‑1 (MQ=255)
                                                                                                                        gCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgc                    >  2:218917/1‑126 (MQ=255)
                                                                                                                               cTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         >  1:153262/1‑114 (MQ=255)
                                                                                                                               cTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         <  2:153262/114‑1 (MQ=255)
                                                                                                                                  tataGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         >  1:552745/1‑111 (MQ=255)
                                                                                                                                  tataGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         <  2:552745/111‑1 (MQ=255)
                                                                                                                                   ataGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAGTAAGCTGTTTTCCCTGGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgctggct               >  2:758617/1‑120 (MQ=255)
                                                                                                                                   ataGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgctggct               <  1:758617/120‑1 (MQ=255)
                                                                                                                                       ccccTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgctg                  <  3:611661/113‑1 (MQ=255)
                                                                                                                                       ccccTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTgctg                  >  4:611661/1‑113 (MQ=255)
                                                                                                                                           tGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGttgttg                                             <  4:354498/82‑1 (MQ=255)
                                                                                                                                           tGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGttgttg                                             >  3:354498/1‑82 (MQ=255)
                                                                                                                                             gTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGACAaaga  <  2:26000/123‑1 (MQ=255)
                                                                                                                                             gTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGACAaaga  >  1:26000/1‑123 (MQ=255)
                                                                                                                                               gagaGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         <  1:184035/98‑1 (MQ=255)
                                                                                                                                               gagaGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCCGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTtctct                         >  2:184035/1‑98 (MQ=255)
                                                                                                                                                  |                                                                                                                     
ATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGAGTTCTACTCCAGTTCCGGACTGCTCACGCCACGGTATTAGGCATATCCTATATAGCCCCTGGTGAGAGTCACCAGTTCCTTGATTAAATAAAATGGAGTTTTACATGAAGGCTTTCAATAAGCTGTTTTCCCTCGTTGTTGCATCTGTTCTGGTTTTCTCTCTTGCTGGCTGCGGTGACAAAGA  >  NC_000913/2211041‑2211304

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