Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,324,837 A→C 11.4% Y31S (TAT→TCT) ‡ 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,8370AC11.4% 245.8 / 23.2 80Y31S (TAT→TCT) ‡yciQenhancer of membrane protein expression; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base C (7/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGA  >  NC_000913/1324690‑1324977
                                                                                                                                                    |                                                                                                                                             
gCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑At                                                                                                                                              >  2:155925/1‑149 (MQ=255)
gCGTAGGTGATGTGAAGCCTTTCTGATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑At                                                                                                                                              <  3:446437/149‑1 (MQ=255)
   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                                                                                                                    <  2:412793/149‑1 (MQ=255)
                          taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                    >  2:364051/1‑149 (MQ=255)
                          taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                    >  3:38597/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                                                                                                                          <  2:348352/117‑1 (MQ=255)
                                                    gACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCtt                                                                                                                          >  1:348352/1‑117 (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:524031/1‑149 (MQ=255)
                                                                                    tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAg                                                         >  2:382859/1‑149 (MQ=255)
                                                                                    tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAg                                                         >  3:303535/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)
                                                                                                                                                    |                                                                                                                                             
GCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAC‑TCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGA  >  NC_000913/1324690‑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: ‑