Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,578,301	0	A	C	41.3%	111.7 / 12.0	63	G142G (GGT→GGG)	gntR	d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (33/4);  new base C (0/26);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 7.67e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATAACC  >  NC_000913/3578153‑3578450                                                                                                                                                        |                                                                                                                                                      cgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATAcc                                                                                                                                                     >  2:815297/1‑151 (MQ=255) cgtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATAcc                                                                                                                                                     >  1:409701/1‑151 (MQ=255)       gaCGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCAc                                                                                                                                               >  1:709739/1‑151 (MQ=255)            gcgcCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCa                                                                                                                                           >  2:1318906/1‑150 (MQ=255)                cGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAatcat                                                                                                                                      >  1:860887/1‑151 (MQ=255)                 gagaTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAatcatc                                                                                                                                     >  2:1136757/1‑151 (MQ=255)                    aTAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCttt                                                                                                                                  >  1:675399/1‑151 (MQ=255)                          atgggggCGAGCCCCGCGAGAAAAAGGGGGAGTGGTCTTCGGGGGTGCTTCTAAAAATTTAAAAAAC‑CCGGGGCCA‑ATAAAAAGGCAGGGCGTTTTGTTGCCCACCTGTCCCCCCCCGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTg                                                                                                                            <  2:1272401/146‑1 (MQ=255)                              gggCGCTCCCCGCGAACAATAAGGGCTGTGGCCATCGGGCTTGCTTCTAAAAATTTAACAAAC‑CCGGCGGCT‑ATATCAAGGCGTGGCCACTTGCTCTCCATCATTCCCCCCGGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                                                                                        <  1:592794/150‑1 (MQ=255)                               ggCGAGCCCCGGGAAAAAAAAGGGCAGTGGTCTTCTGGCGTCCTCCTAAAAATTAAAAAAAA‑CCGACGGCA‑ATAACAAGGCAGGCCGATTTGCTCTCCACCATTCCCCCCAGGGA‑AAAACCCGCCACTTCAATCATCTTTAAGGTgcgc                                                                                                                         <  2:1036797/149‑1 (MQ=255)                               ggCGACGCCCGGGAGAAATAATGGCAGTGGTCATCGGGCGTGCTGTTTAAAATTTAAAAAAC‑CCGACGGCG‑ATATCAAGGCAGGGCGATTTGCTGCCCATCATTCCCCCCACGGG‑AATCCCCGCCACTTCAATCATCTTTAAgg                                                                                                                              <  1:1169898/144‑1 (MQ=255)                                  gATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgt                                                                                                                    >  1:500500/1‑151 (MQ=255)                                   aTGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtg                                                                                                                   >  2:648005/1‑151 (MQ=255)                                         gcgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGt                                                                                                              >  2:801344/1‑150 (MQ=255)                                          cgAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTAc                                                                                                            >  2:307317/1‑151 (MQ=255)                                               aaaaagggggggggTCTTCGGGCGTGTTGCTT‑AAAGTTTACAAAAGCCGACGGCG‑ATACAAAGGGATGGCAATTGGTTCTCCATCAGTCCCACCGGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGtt                                                                                                         <  1:930091/137‑1 (MQ=255)                                                 taatGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                                     >  2:265910/1‑151 (MQ=255)                                                     ggagggggtattcgggCGTGCTACTAAAAATTAAAAAAAC‑CCGAGGGGT‑ATATAAGGGCAGGGCGACTTGTTGTCCAGCATTCCCCCCGCGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                                                                                  <  1:785133/136‑1 (MQ=255)                                                     ggCGGTGTTCATTTGGCGTGCTCCTAAAAATTTACAAAAC‑CCGCCGGCA‑ATATCAGGGCAGGGCGCCTTGCTGCCCACCTTTCCCCCCAGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                                                                                    <  2:710949/148‑1 (MQ=255)                                                     ggCGGTGGTCATCGGGCGTGCTGCTTAAAGGTAAAAAAAC‑CCGGCGGCT‑ATACAAAGGAAGGGCAACTTGCTGTCCATCAGTCCCCCCCCGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGt                                                                                                     <  1:859479/147‑1 (MQ=255)                                                      ggggggttaatgtGGCGTTCTGTTAAAAATTTAAAAAAC‑CCGAGGGGA‑AAATCAGGGCAGGGCGCTTTGCGGCCCCTCATTCCCCCCGGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGa                                                                                                <  2:1090637/139‑1 (MQ=255)                                                         ggggtatttgggCGTGTTGTTTAAAATTTACAAAAC‑CCGAGGGCG‑ATATAAGGGCAGGGCGATTTGCTGTCCATCAGTCCCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                                             <  1:1188800/141‑1 (MQ=255)                                                          tgctCATCTGGCGTGCTTCTAAAAATTTAAAAAAC‑CCGACGGCG‑ATAACAAGGCATGGCGATTTGTTGTCCATCATTCCCCCCGGGGA‑AACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa                                                                                             <  2:210122/147‑1 (MQ=255)                                                          tGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATc                                                                                              >  2:510804/1‑149 (MQ=255)                                                           ggTCTTCTGGGGGGCTGCTTAAATTTTAAAAAAC‑CCGACGGCA‑ATATAAGGGCAGGGCGCTTTGCTTTCCATCAGTCCCCCCAGGGA‑AACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAgg                                                                                           <  1:961055/151‑1 (MQ=255)                                                                  gggcttgttcttaaaaaGTTATAAAAC‑CCGCCGGCG‑ATAAAAAGGCAGGGCTACTTCTTGTCCATCATTCCCCCCAGGGA‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGa                                                                                    <  2:155707/137‑1 (MQ=255)                                                                   ggCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                                                                                   >  1:971048/1‑151 (MQ=255)                                                                      gTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACTCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATAtt                                                                                <  2:348589/151‑1 (MQ=255)                                                                       tgctgcTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATAtt                                                                                >  1:130081/1‑150 (MQ=255)                                                                       tgctgcTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCACGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTc                                                                               >  2:672312/1‑151 (MQ=255)                                                                       tgctcttaaaaaGTTATAAAAC‑CCGGCGGCG‑AAATCAAGGCAGGGCGCCTTGTTGTCCATCAGTCCCACCGGGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATAtt                                                                                <  1:80031/141‑1 (MQ=255)                                                                          tgctTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTcc                                                                              >  1:995507/1‑149 (MQ=255)                                                                               cAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgag                                                                         >  1:550470/1‑149 (MQ=255)                                                                                aaaGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGc                                                                      >  1:960131/1‑151 (MQ=255)                                                                                 aaGTTTACAAAA‑CCAGCGGGGAATATCAGGGCCGGGCGCCTTCCTGCCCACCATTCCCCCCACGGA‑AACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg                                                                           <  2:698968/146‑1 (MQ=255)                                                                                 aaGCTATCTAAA‑CCGACGGCC‑ATATCAAGGCAGGGCGACCTGCTGTCCATCAGTTCCACCAGGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGCCCGGTGAGGCTCAGGCCGTCGAGCTTCCAGGcgag                                                                       >  2:142475/1‑149 (MQ=255)                                                                                   gTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATg                                                                   >  1:1276078/1‑151 (MQ=255)                                                                                    taaacAAAC‑CCGACGGCG‑ATATAAGGGCATGGCGTTTTTCTCTCCACCATTCCCCCCGGGGG‑AACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCAt                                                                    <  2:868863/145‑1 (MQ=255)                                                                                      accAAAC‑CCGGCGGCA‑ATATCAGGGCGGGGCGCTTTGCTTCCCATCAGTCCCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                                <  2:640877/149‑1 (MQ=255)                                                                                         aaaa‑CCGCCGGCA‑AAAAAAAGGCATGGCGCCTTGCTGTCCATCA‑TTCCCCCCCGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt                                                                <  1:40406/148‑1 (MQ=255)                                                                                          aaa‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCg                                                             >  2:71962/1‑150 (MQ=255)                                                                                            a‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGCTCAGGCCGTCGATCTTCCAGGAGAGCATGGATTCGAgg                                                          >  1:18484/1‑151 (MQ=255)                                                                                               cGACGGCG‑ATATCAGGGATTGGCGTCTTGCGGTCCATCAGTCCCCCCACGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCg                                                        <  1:932352/151‑1 (MQ=255)                                                                                                gACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGt                                                       >  2:63149/1‑151 (MQ=255)                                                                                                        atacaaaGGCATGGCTACTTGCTGGCCATAAGTCCCCCCCGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCt                                                <  1:1158689/146‑1 (MQ=255)                                                                                                        ataaCAAGGCGGGGCGTTTTGCTTCCCATCAGTCCCCCCGGGGG‑AACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCt                                                <  2:925489/151‑1 (MQ=255)                                                                                                               ggcggggCGCCTTGTTGCCCATCATTCCCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAttt                                          <  1:979321/145‑1 (MQ=255)                                                                                                                   tGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                                     >  2:675240/1‑151 (MQ=255)                                                                                                                    ggCGCCTTGCTCCCCACCAGCCCCCCCGGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGt                                     <  2:346340/150‑1 (MQ=255)                                                                                                                       gACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtat                                 >  1:234570/1‑151 (MQ=255)                                                                                                                          ttGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATAATCTTTAAGGTGCGCGGCGTGTGGGAACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTtata                                >  1:1019022/1‑149 (MQ=255)                                                                                                                            gtttcccATCAGTCCCCCCGCGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAAc                              <  2:108680/144‑1 (MQ=255)                                                                                                                               gTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCg                            >  1:449863/1‑148 (MQ=255)                                                                                                                                tCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGt                        >  2:368028/1‑151 (MQ=255)                                                                                                                                 ccATCATTCCCCCCGCGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAg                         <  1:251512/149‑1 (MQ=255)                                                                                                                                 ccATCAGTCCCCCCAGGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTa                          <  1:1335795/148‑1 (MQ=255)                                                                                                                                  caacattccccccccGGT‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAg                         <  2:1363000/134‑1 (MQ=255)                                                                                                                                          cccccagGGG‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGc                <  2:38216/146‑1 (MQ=255)                                                                                                                                          cccccagGGA‑AATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGc                <  1:326700/146‑1 (MQ=255)                                                                                                                                                 ggg‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTg        >  1:1277753/1‑150 (MQ=255)                                                                                                                                                  gg‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTATATAACCGTAGTGCGCCAGCATGGTCTGAt      >  1:790444/1‑151 (MQ=255)                                                                                                                                                     aaTACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATa     >  1:242795/1‑150 (MQ=255)                                                                                                                                                       tACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGCGTGGGAACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATAAcc  >  2:1346792/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      CGTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAA‑CCGACGGCG‑ATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGG‑AATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCGTAGTGCGCCAGCATGGTCTGATAACC  >  NC_000913/3578153‑3578450

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