Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,578,2870TC41.9% 94.7 / 13.9 62E147G (GAA→GGA) gntRd‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (28/8);  new base C (0/26);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGC  >  NC_000913/3578148‑3578436
                                                                                                                                             |                                                                                                                                                      
aCGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACggg                                                                                                                                              <  2:972554/150‑1 (MQ=255)
      gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAccc                                                                                                                                      <  2:882892/151‑1 (MQ=255)
      gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATAccc                                                                                                                                      <  1:961858/151‑1 (MQ=255)
        cgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGc                                                                                                                                    >  2:1011806/1‑151 (MQ=255)
         gagaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGACGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGcc                                                                                                                                   >  1:237063/1‑151 (MQ=255)
          agaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCACAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGcc                                                                                                                                   >  1:639208/1‑150 (MQ=255)
          agaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGcc                                                                                                                                   >  1:88078/1‑150 (MQ=255)
           gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCa                                                                                                                                  <  2:869650/150‑1 (MQ=255)
             cGTGCGCCGAGATAGACAATGTGGCGATGCCCGCGAGCAATAATGGCAGGGGTCATTTGGCGTGCTGCTCAAAAGTTATAAAAACCGACGGCGATATCAAGGCATGGCGTCTTGCTGTCCATCAGTTCCCCCACGGGAATTCCCGCCACtt                                                                                                                               <  2:776001/151‑1 (MQ=255)
                  gcCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatc                                                                                                                          >  2:935008/1‑151 (MQ=255)
                        atgggaatgggggCGAGCCCGGCGAAAAATAAGGGCGGGGGTCCTCGGGGGTGCTCCTTAAAATTTAAAAAAACCGAGGGGAATATCAAGGCATGGCGCTTTTCTTTCCATCAGTCCCCCCACGGGAATACCCGCCACTTCAATCATCttt                                                                                                                    <  1:977736/140‑1 (MQ=255)
                        aggggCAATGGGGCGATGCCCGCGAGAAATAATGGCAGTGGTCACCGGGCGTGGTCCTCAAAATTGAAAAAACCCGACGGCAATAACAAGGCAGGGCGCTTTGCTGCCCAGCTGTCCCACCAGGGGAATCCCCGCCACTTCAATCATCtt                                                                                                                     <  1:832813/147‑1 (MQ=255)
                          aGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTa                                                                                                                   >  2:630312/1‑150 (MQ=255)
                           ggatgtggggCGGTGCCCGAGAAAAAAAATGGCGGGGTTCTTCTGGCGTGCTGCTAAAAATTTATAAAACCCGACGGCGATAACAAGGCAGGGCGTTTTGCTCTCCATCAGTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTaa                                                                                                                  <  1:544881/142‑1 (MQ=255)
                           gacAATGTGGCGATGCCCGCGAGCAATAAGGGCGGTGGTCTTCGGGCGTGCTGCTACAAAGTTATCAAACCCGACGGCGATATCAAGGCATGGCGACTTGCTCTCCATCAGTCCCACCACGGAAATCCCCGCCACTTCAATCATCTTTa                                                                                                                   <  1:267398/147‑1 (MQ=255)
                                   ggCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGACGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                                                                         >  1:844771/1‑151 (MQ=255)
                                   ggCGAGGCGCGCGAAAAAAAAGGGAAGGGGTCATCGGGCGTGCTGCTTAAAATTTAAAAAACCCGACAGCGATATCAAGGCAGGGCGACTTGCTGTCCAACATTCCCCCCCCGGGAATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                                <  1:827751/144‑1 (MQ=255)
                                   ggCGAGGCCCGGGAGAAATAAGGGCGGTGGTCTTCTGGCGTGCTGCTTAAAATTTATAAAACCCGACGGCGATATCAAGGCATGGCGTTTTGCTGCCCACCAGTCCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTgcg                                                                                                            <  2:828751/148‑1 (MQ=255)
                                      gtcccccGGGAGAAATAGGGGGGGTGTTCTTCGGGGGTGCTCCTAAAAATTTAAAAAACCCGCCGGCGATAACAGGGCGGGGCGACTTGCTGTCCATCATTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgt                                                                                                      <  2:30880/147‑1 (MQ=255)
                                      gATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgt                                                                                                      >  2:293214/1‑151 (MQ=255)
                                        tGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtgt                                                                                                    >  1:838138/1‑151 (MQ=255)
                                         gCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGtgtg                                                                                                   >  2:715453/1‑151 (MQ=255)
                                          cccGCGAGAAATAAGGGCAGTGGTCATCTGGGGTGCTGCTTCAAAGTTAAAAAACCCGACGGCGATATCAAGGCATGGCGATTTGCTGTCCATCAGTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTgg                                                                                                  <  2:346118/151‑1 (MQ=255)
                                           ccccgAGAAAAAAGGGGAGGGGTCTTCTGGCGGGCTGCTAAAAATTAACAAAACCCGACGGCAATACAAGGGCTGGGCTATTGGCTGCCCATCGTTCCCCCCAGGAAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTggg                                                                                                 <  1:608912/148‑1 (MQ=255)
                                             gcgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGt                                                                                                >  2:385770/1‑150 (MQ=255)
                                              cgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTAc                                                                                              >  2:10145/1‑151 (MQ=255)
                                                agaaataagGGCAGTGGTCATCGGGCGTGCTGCTAAAAATTTAAAAAACCCGGCGGCAATACCAGGCCAGGGCGTCTTGCTGTCCATCATTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTAc                                                                                              <  1:696936/146‑1 (MQ=255)
                                                aGCAATAAGGGCGTTGGTCCTCGGGCGGGCTGCTAAAAAGTTATAAAAACCGACGGCAATATCAAGGCATGGCGTCTTGCTGTCCATCAGTCCCCCCGCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTgggc                                                                                                <  1:781226/147‑2 (MQ=255)
                                                 gCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGtt                                                                                           >  2:1001743/1‑151 (MQ=255)
                                                  cAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                                                                          >  1:589725/1‑151 (MQ=255)
                                                   aatattGGAAGTGTTCTTCGGGCGTGCTTCTTAAATTTTTTAAAACCCGACGTCAATACCGCGGAGTGGCGACTTTCCGCCCCTCATTCCCCCCGGGGGAGTACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                                                                          <  2:2021/150‑1 (MQ=255)
                                                   aaaaagggCAGTGTCCATGTGGCGTGCTGCTTAAAAGTTAAAAAACCCGACGGCGATAAAAAGGCTTGGCCACTTGCTGCCCATCATTCCCCCCCGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                                                                          <  1:528208/144‑1 (MQ=255)
                                                         ggCAGTGTTCATCTGGGGTGCTGCTACAAATTTAAAAAACCCGACGGGGATACAAGGCAAGGGCGCTTTGTTTCCCCTCGTCTCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGa                                                                                     <  1:17277/149‑1 (MQ=255)
                                                          gCAGTGTTCATCTGGCGTGCTGCTTCAAAGTTATCAAACCCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGa                                                                                  <  2:689958/151‑1 (MQ=255)
                                                              tGGTCTTGTGGCGTGCTGCTTAAAAGTTATCAAACCCGACGGCAATATAAAGGCGGGGCGATTTGCTGTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                               <  2:510098/150‑1 (MQ=255)
                                                              tGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                >  2:184026/1‑149 (MQ=255)
                                                               gttCATCGGGGGTGCTGCTTAAAATTTACAAAACCCGACGGCGAAATAAAGGCATGGCGCCTTGCTGTCCATCATTCCCCCCAGGGAAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAg                                                                              <  2:462285/148‑1 (MQ=255)
                                                                     cTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                       >  2:80960/1‑151 (MQ=255)
                                                                       ggggggctgctAAAAATTTAAAAAACCCGCCGGCGATATCAAGGAATGGCGTTTTGTTCCCCACCATTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                                                                     <  2:724105/146‑1 (MQ=255)
                                                                       ggCGTGCTCTTTAAAATTTAAAAAACCCGCGGGCAATATAAGGCCAGGGCGCTTTGCTCTCCATCATTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                                                                     <  2:782945/151‑1 (MQ=255)
                                                                       ggCGTGCTCCTACAAATTTAAAAAACCCGACGGCAATAACAGGGCATGCGCACTTCCTCCCCATCACTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                                                                       <  1:434540/149‑1 (MQ=255)
                                                                        gCGTGCTGTTAAAAATTTAAAAAACCCGACGGCTATACCAAGGCATGGCGCTTTGTTGTCCATCATTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                <  2:705292/139‑1 (MQ=255)
                                                                          gTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATAtt                                                                  >  1:972279/1‑151 (MQ=255)
                                                                                  tCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                             >  2:1933/1‑148 (MQ=255)
                                                                                    aaaTTTATCAAACCCGAGGGCAATACCAAGGCATGGCGACTTGTTGTCCATCATTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGagag                                                         <  1:532264/150‑1 (MQ=255)
                                                                                    aaaGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGc                                                        >  2:577107/1‑151 (MQ=255)
                                                                                       gTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                     >  2:297739/1‑151 (MQ=255)
                                                                                       gTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                     >  2:1032155/1‑151 (MQ=255)
                                                                                       gTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                     >  1:296104/1‑151 (MQ=255)
                                                                                          accAAACCCGACGGGGATACCAGGGCATGGCGATTTGCTGTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg                                                    <  1:140021/147‑1 (MQ=255)
                                                                                          aaaaaacCCGACGGAAATATAAAGGCATGGCGCCTTGCTGTCCATCATTTCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGa                                                   <  2:945060/147‑1 (MQ=255)
                                                                                                  ccGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGc                                           >  1:41274/1‑151 (MQ=255)
                                                                                                   cGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGAGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCg                                          >  1:260776/1‑151 (MQ=255)
                                                                                                       ggCAATATAAAGGCAGGGCGATTTGCTGTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTc                                       <  1:165178/150‑1 (MQ=255)
                                                                                                           atatCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCt                                  >  2:367211/1‑151 (MQ=255)
                                                                                                             atCACGGCACGGCGACTTGCTGCCCATCAGTTCCACCACGGGAATACCCGACACTTCAATCATCTTTAAGGTGCGCGGCGCGTGGGTACGTTCGGTGAGGATCAGGCCGTCGCTATTCCAGGAGAGCATGGATTCGAGGCGTTCTCGCTcc                                >  1:153702/1‑151 (MQ=255)
                                                                                                              agaaggcctcggCGATTTGCTTCCCTCCAGTCCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc                                   <  2:955792/138‑1 (MQ=255)
                                                                                                                  ggCATGGCGACTTGTTGTCCACCTGCCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAttt                            <  1:837616/150‑1 (MQ=255)
                                                                                                                  ggCAGGGCTATTTGCTCCCCTTGAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAt                              <  1:730355/148‑1 (MQ=255)
                                                                                                                        ggcgttttGCTGCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCgg                        <  2:569612/143‑1 (MQ=255)
                                                                                                                              ttGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATaa                 >  1:630974/1‑150 (MQ=255)
                                                                                                                                     ccATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTg         >  1:894697/1‑151 (MQ=255)
                                                                                                                                            ttCCCCCGCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCa    <  2:221891/149‑1 (MQ=255)
                                                                                                                                            ttCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGc  >  2:726649/1‑151 (MQ=255)
                                                                                                                                             tCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGc  >  2:982879/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
ACGTTCGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGC  >  NC_000913/3578148‑3578436

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