Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,324,839 G→T 12.2% E32* (GAA→TAA)  yciQ → enhancer of membrane protein expression; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,324,8390GT12.2% 243.9 / 19.6 77E32* (GAA→TAA) yciQenhancer of membrane protein expression; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base T (7/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGA  >  NC_000913/1324693‑1324977
                                                                                                                                                   |                                                                                                                                           
tAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaa                                                                                                                                           <  1:339699/149‑1 (MQ=255)
     gATGTGAAGGCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAgg                                                                                                                                      <  3:343697/149‑1 (MQ=255)
       tgtgAAGCCTTTCTTATAACAAATTGCGTTCCGCGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCCTTTTGCTGTTCATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAGGGa                                                                                                                                    <  4:93054/149‑1 (MQ=255)
          gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGAtct                                                                                                                                 >  2:409419/1‑149 (MQ=255)
           aaGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTATATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGAtctc                                                                                                                                <  1:510908/149‑1 (MQ=255)
             gCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGc                                                                                                                              >  4:43998/1‑149 (MQ=255)
               cTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGctc                                                                                                                            <  3:233012/149‑1 (MQ=255)
                tttCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGctct                                                                                                                           >  4:568960/1‑149 (MQ=255)
                tttCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGctct                                                                                                                           <  1:460374/149‑1 (MQ=255)
                       taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCAT‑AAGAAggg                                                                                                                                     <  2:339172/132‑1 (MQ=255)
                       taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAggg                                                                                                                                     >  1:339172/1‑132 (MQ=255)
                       taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                    >  3:38597/1‑149 (MQ=255)
                       taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                    <  2:412793/149‑1 (MQ=255)
                       taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                    >  2:364051/1‑149 (MQ=255)
                         aCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCttt                                                                                                                         <  3:357279/142‑1 (MQ=255)
                         aCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCttt                                                                                                                         >  4:357279/1‑142 (MQ=255)
                           aaaTTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGaa                                                                                                                <  3:314997/149‑1 (MQ=255)
                           aaaTTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGaa                                                                                                                >  4:300561/1‑149 (MQ=255)
                                gCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAc                                                                                                           <  2:550058/149‑1 (MQ=255)
                                   ttCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAc‑tct                                                                                                       >  4:461327/1‑149 (MQ=255)
                                       aCGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCaa                                                                                                   >  3:476365/1‑149 (MQ=255)
                                       aCGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCaa                                                                                                   <  2:171311/149‑1 (MQ=255)
                                         ggATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAGGGAt                                                                                                                                   >  2:283948/1‑116 (MQ=255)
                                         ggATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAGGGAt                                                                                                                                   <  1:283948/116‑1 (MQ=255)
                                            tGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGa                                                                                              <  3:300561/149‑1 (MQ=255)
                                             ggAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCCC‑TCTGCAACCAGAt                                                                                             <  3:43998/149‑1 (MQ=255)
                                                aGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATggg                                                                                          <  4:511419/149‑1 (MQ=255)
                                                 gACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCtt                                                                                                                          >  1:348352/1‑117 (MQ=255)
                                                 gACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCtt                                                                                                                          <  2:348352/117‑1 (MQ=255)
                                                  aCTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGt                                                                                        >  3:500632/1‑149 (MQ=255)
                                                      tGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTc                                                                                    <  2:54999/149‑1 (MQ=255)
                                                               aGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCaaa                                                                           >  1:300629/1‑149 (MQ=255)
                                                               aGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCaaa                                                                           <  3:461327/149‑1 (MQ=255)
                                                                gTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAg                                                                          >  2:82641/1‑149 (MQ=255)
                                                                  ttCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGaa                                                                        >  4:584016/1‑149 (MQ=255)
                                                                               gCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTc                                                           >  1:143014/1‑149 (MQ=255)
                                                                                 tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCAT‑AAGAAGGGAt                                                                                                                                   <  2:496514/76‑1 (MQ=38)
                                                                                 tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAGGGAt                                                                                                                                   >  1:496514/1‑76 (MQ=38)
                                                                                 tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAg                                                         >  3:303535/1‑149 (MQ=255)
                                                                                 tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAg                                                         >  3:524031/1‑149 (MQ=255)
                                                                                 tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAg                                                         >  2:382859/1‑149 (MQ=255)
                                                                                  gTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAgg                                                        >  1:359239/1‑149 (MQ=255)
                                                                                   ttGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGgc                                                       >  2:517715/1‑149 (MQ=255)
                                                                                   ttGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGgc                                                       >  1:440971/1‑149 (MQ=255)
                                                                                    tgaaataaGGGTTGTTTGTATCACCTCTAGGTTATGAAAAAACCCGGGAGACCCTTTTTT‑TGGAGGAAGGATTCCCCTTCTTGGATGTT‑AAGCCACTTCTCAAACAAAATGGGGTGC‑CGAAAACAAAAAAAATATTCATTTTCAGgcgc                                                     <  3:284865/142‑1 (MQ=255)
                                                                                             ggCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATcag                                             <  4:476365/149‑1 (MQ=255)
                                                                                                  tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGAt                                        >  1:158213/1‑149 (MQ=255)
                                                                                                  tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGAt                                        >  4:317250/1‑149 (MQ=255)
                                                                                                         catcTCTAAGTTATGCAGAAAACACGTAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCAc                                 <  4:98751/149‑1 (MQ=255)
                                                                                                          atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAAt                                                >  3:120359/1‑133 (MQ=255)
                                                                                                          atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAAt                                                <  4:120359/133‑1 (MQ=255)
                                                                                                               tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCa                                                                                                <  1:555025/80‑1 (MQ=255)
                                                                                                               tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCa                                                                                                >  2:555025/1‑80 (MQ=255)
                                                                                                               tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGa                                         <  1:224463/135‑1 (MQ=255)
                                                                                                               tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGa                                         >  2:224463/1‑135 (MQ=255)
                                                                                                               tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCttt                           >  2:196568/1‑149 (MQ=255)
                                                                                                               tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCttt                           >  3:581218/1‑149 (MQ=255)
                                                                                                                 aGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTa                         <  3:584016/149‑1 (MQ=255)
                                                                                                                 aGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTa                         <  4:243594/149‑1 (MQ=255)
                                                                                                                  gTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTAt                        >  2:388188/1‑149 (MQ=255)
                                                                                                                         aGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGAtt                 >  1:499053/1‑149 (MQ=255)
                                                                                                                           aaaaCACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTa               >  3:225461/1‑149 (MQ=255)
                                                                                                                               cacGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCAc           <  3:568960/149‑1 (MQ=255)
                                                                                                                                acGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACg          >  2:319168/1‑149 (MQ=255)
                                                                                                                                     gaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATACTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTAt     <  1:155925/149‑1 (MQ=255)
                                                                                                                                      aTCCCTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATg    >  4:426201/1‑149 (MQ=255)
                                                                                                                                      aTCCCTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATg    >  4:479697/1‑149 (MQ=255)
                                                                                                                                      aTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATg    <  2:417491/149‑1 (MQ=255)
                                                                                                                                       tCCCTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATgg   >  4:283903/1‑149 (MQ=255)
                                                                                                                                       tCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATgg   <  4:500632/149‑1 (MQ=255)
                                                                                                                                       tCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATgg   <  1:388188/149‑1 (MQ=255)
                                                                                                                                        cccTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGa  >  3:218714/1‑149 (MQ=255)
                                                                                                                                          cttctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCATCCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGa                   <  4:347467/130‑1 (MQ=255)
                                                                                                                                          cttctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATc                       <  4:215620/126‑1 (MQ=255)
                                                                                                                                          cttctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATc                       >  3:215620/1‑126 (MQ=255)
                                                                                                                                          cttctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGa                   >  3:347467/1‑130 (MQ=255)
                                                                                                                                          cttctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTAt     >  2:106651/1‑144 (MQ=255)
                                                                                                                                          cttctt‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTAt     <  1:106651/144‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGA  >  NC_000913/1324693‑1324977

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