Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,391	0	A	C	43.5%	78.5 / 28.4	62	intergenic (+93/‑54)	dnaJ/insL1	chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base A (7/28);  new base C (27/0);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 2.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCC  >  NC_000913/15242‑15528                                                                                                                                                      |                                                                                                                                                cgcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAt                                                                                                                                                <  2:611845/151‑1 (MQ=255)     cAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAg                                                                                                                                             <  1:572982/150‑1 (MQ=255)        agagCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCAATCCCCTaaac                                                                                                                                            >  1:331390/1‑146 (MQ=255)         gagCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑a                                                                                                                                       <  1:188047/151‑1 (MQ=255)              tcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑t                                                                                                                                   <  1:137224/149‑1 (MQ=255)                  tGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAgg                                                                                                                              <  2:1056917/150‑1 (MQ=255)                   gATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAgg                                                                                                                              <  1:487183/149‑1 (MQ=255)                    aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGtt                                                                                                                           <  2:848778/151‑1 (MQ=255)                      ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCCAAACGCC‑AAC‑tttaagg                                                                                                                             >  1:773651/1‑142 (MQ=255)                       gtgtGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTtgtg                                                                                                                        <  2:1161124/151‑1 (MQ=255)                               aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACg                                                                                                                <  2:908317/151‑1 (MQ=255)                                   tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCCAAGCGGT‑AAC‑TTCAAGGGTGTGGGAtttgacctg                                                                                                            >  2:619252/1‑144 (MQ=255)                                        gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGatat                                                                                                        <  1:1165165/150‑1 (MQ=255)                                        gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑g                                                                                                      <  2:480296/151‑1 (MQ=255)                                         acgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACTCCCCTAAAAGCG‑AAC‑TTTAAGGGTGGGGGGTTTTGCCTGatat                                                                                                        >  2:192697/1‑149 (MQ=255)                                           gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCAATCCCCTAAACGCC‑AAC‑TTTAAGGGTGTGGGATTAAGCCTGata                                                                                                         >  1:641841/1‑146 (MQ=255)                                           gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑tt                                                                                                  <  1:1012868/151‑1 (MQ=255)                                           gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACCCCCCTAAACGCC‑AAC‑TTTAGGGGTGGGGGTTTAAGcccgg                                                                                                           >  1:526201/1‑141 (MQ=255)                                                  cccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACTCCCATAACCGCA‑AAA‑TTTAAGGGTATGGTATTataaccgg                                                                                                           >  1:766603/1‑130 (MQ=255)                                                     gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAACGCT‑AAC‑TTAAAGGGTGTGGGATTACGCCTGGTAT‑TA‑TTTTACGTGcc                                                                                         >  1:344369/1‑150 (MQ=255)                                                       tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATGACTCCCCTAAACGCT‑AAC‑TTAAGGGGTGGGGGATTAAGCCTGGTAT‑GA‑TTTaa                                                                                               >  2:262882/1‑142 (MQ=255)                                                        aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATg                                                                                     <  2:513270/151‑1 (MQ=255)                                                                cAAAAGCCTTCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTActc                                                                             <  2:771070/151‑1 (MQ=255)                                                                       cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCTAAGAGCT‑TAC‑TTTAGGGGTGTGGGATTTTGCCTGCTAT‑TA‑TTTTAAGTGCCGGAGAAATATTCTCA‑cca                                                                      >  1:997764/1‑148 (MQ=255)                                                                       cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTAACCCCCCTAAACGCT‑AAC‑TTAAAGGGTGTGGGATTTAGCCCGGTAT‑GA‑TTTTAAGTGCCGCAGAAATAATCTCA‑cca                                                                      >  1:1115959/1‑148 (MQ=255)                                                                         gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATa                                                                    <  2:1060354/150‑1 (MQ=255)                                                                            cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAAAGCT‑AAC‑TTAAGGGGTGTGGGGTTACGCCTGGTTT‑GA‑TTTTAAGTGCCGATGAATTATTCTCC‑CGATaa                                                                   >  2:955792/1‑148 (MQ=255)                                                                              tGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACTCCCCCAAAAGCG‑AAC‑TTTAGGGGTGGGGGATTACGCCCGGTAT‑TA‑TTTTAAGTGCCGACGGAATAATCTCA‑CGATaaattg                                                               >  2:441741/1‑146 (MQ=255)                                                                                  cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGTGGTGACTCCCCACAGCGCTGAAC‑TTATGGGGT‑TGGTGTTACCACTGATGT‑GA‑TTTTAAATGGCGATGAAGTAGTTTCA‑CGaaa                                                                    >  1:151485/1‑139 (MQ=255)                                                                                   aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACCCCCATAAACGCT‑AAC‑ATAAGGGGTGTGGTATTTAACCCGATAT‑GA‑TTTAAAGTGCCGCTGGATTTATCTCA‑CGCTAAATGGt                                                              >  2:741270/1‑146 (MQ=255)                                                                                   aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑a                                                        <  2:641823/151‑1 (MQ=255)                                                                                      ccTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCAAAACGCT‑AAC‑TTAAGGGGTGTGGGGTTTTGCCTGATAT‑GA‑TTTTAAGTGCCGGTGAATTATTCTCA‑CGATAACTAGTCAGC‑AAtt                                                     >  2:278284/1‑151 (MQ=255)                                                                                       cTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCC‑AAA‑TTTAGGGTTGTGGTGTTTCGCCTGGTCT‑TC‑TTTTACCTGCCGATGAATTACTCTCA‑CGGTAACTGGTCCGC‑aaat                                                     >  1:970402/1‑148 (MQ=21)                                                                                        tGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATGACTCCCCTAAACGCC‑AAA‑TTTAGGGGTTTGGGATTACGCCTGATAT‑GA‑TTTAAAATGACCGTGAATTAATCTCA‑CAAGAAATTGTCAGC‑AATTCt                                                   >  1:384316/1‑151 (MQ=21)                                                                                          ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTgg                                                 <  2:992057/151‑1 (MQ=255)                                                                                            tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAGCGCT‑AAC‑TTAAGGGGTGTGGGATTTCGCCTGGTAT‑GA‑TTTAACGTGCCGCTGAATTAATTTCA‑CGATAACCGGTCAGC‑AAATCTTGcc                                               >  1:989941/1‑151 (MQ=255)                                                                                            tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGCACTCCCCCAAAAGCT‑AAC‑CTAAAGGGTGTGGGATTAAGCCTGGTAT‑GA‑TTTTAAGTGGCGATGAAATAATCTCA‑CGATAAATGGTCAGC‑AAtt                                                     >  2:246107/1‑145 (MQ=255)                                                                                             aaaaaTAGGGTGCGTTGAAGATATGCGAGCGCCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGcc                                               <  2:289663/150‑1 (MQ=255)                                                                                             aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATGACTCCCCTAAAAGCTAAAC‑TTAAGGG‑TGTGGGGTTTAGCCCGGTAT‑TA‑TTTTAACTGCCGGGTAAATACTCTCA‑CGATAAATGGGTAGCAAATTTTGGcc                                               >  2:188847/1‑151 (MQ=255)                                                                                                aaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCata                                           <  2:875599/151‑1 (MQ=255)                                                                                                 aTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCAAAACGCC‑AAC‑TTAAAGGGTGGGGGATTTCGCCCGGTAT‑TA‑TTTTAAGTGCCGAAGAAATACTCTCA‑CCATAAATTGGCAGC‑CATTCTTGGCCctat                                          >  1:1050302/1‑151 (MQ=255)                                                                                                   aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAACGCC‑AAC‑CTAAAGGGTGGGGTATTACGCCCGATAT‑TA‑TTTTACGTGCCGGTGAAATAATCTCC‑CGATAACTTGTCAGC‑AAATCTGGCCCCTAtt                                         >  1:985197/1‑150 (MQ=255)                                                                                                   aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTg                                        <  1:1080852/151‑1 (MQ=255)                                                                                                   aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCC‑AAC‑TTAAGGGGTGTGGGATTACACCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTTGCCCCTAtt                                         >  1:771113/1‑150 (MQ=18)                                                                                                    gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACCCCCACAAACGCT‑AAC‑TTAAGGGGTGTGGGGTTACGCCTGGTAT‑TA‑ATTTACATCCCGATGAATTAATCTCC‑CGCTAAATGGTCAGC‑AATTTTTGCCCCTAtt                                         >  1:629256/1‑149 (MQ=255)                                                                                                       tGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACACCCCAAAAAGCC‑AAC‑TTTAGGGGGGGGGGGTTTTGCCTTGTAT‑GATTTTAACG‑GCCGATGAAATAATCTCACCGAAAAC‑GGTCAGG‑AAATTTGGCCCCTAtt                                         >  2:620190/1‑146 (MQ=255)                                                                                                        gCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑CAC‑TTTAAGGGTGTGGGGTTACGCCCGGTAT‑GA‑TTTTACGTGCCGGTGAATTAATCTCA‑CGATAAATGGTCAGC‑AATTCTGGCCCCTAtt                                         >  1:446227/1‑145 (MQ=255)                                                                                                         cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCCAAACGAC‑AAC‑TTTAGGGGTGTGGTGTTATGCCCGGTAT‑GA‑TTTTACGTGGCGGTGAATTAATCTCA‑CCATAACTGGGCAGC‑AATTCTGGGCCCTATTTGt                                      >  2:864765/1‑147 (MQ=255)                                                                                                         cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAACGCT‑AAC‑TTAAGGGGTGTGGTATTAAGCCCGGTAT‑TA‑TTTTAAGTGCCGGTGAATTAATCTCA‑CCAAAAATTGTCAGC‑AATTTTTGCCCCTATTGGTaa                                    >  1:831600/1‑149 (MQ=255)                                                                                                             gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAACGCC‑AAA‑TTAAAGGGTGTGGGATTATGCCTGGTAT‑GA‑TTTTAAGTGCCGGTGAATTACTCTCA‑CGATAAATGGGCAGC‑AATTTTTGCCCATATTGGTAAGccc                                >  1:94098/1‑149 (MQ=255)                                                                                                               aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCgaag                            <  1:637503/151‑1 (MQ=255)                                                                                                                gATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCgaaga                           <  2:341430/151‑1 (MQ=255)                                                                                                                 atatGCGAGCACCTGTAAAGTGGCGGGGGGTACCCCCCAAAAAGCC‑AACTTTAAGGGTT‑TGGGTTTAAGCCTGGTATGGA‑TTTAA‑GTTCCGAAGAAATTATCTCA‑CGCTAAAAGGGCAGC‑AAATCTTGCCCCTATTGGTAAAcccccaga                           >  1:38605/1‑145 (MQ=255)                                                                                                                  tatGCGAGCACCTGTAAAGTGGCGGGGATCACCCCCCCAAACGCTAAAT‑TTAAGGGTTGGGGTTTTA‑GCCCTATAT‑GA‑TTTAAAAGTCCCATTAAATACTCTCA‑AGATAACATGTCAGC‑CAATTTGTCCCCAATTTGTTAACCCcaagaa                          >  1:735881/1‑150 (MQ=255)                                                                                                                   atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACt                        <  2:1006729/151‑1 (MQ=255)                                                                                                                     gCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg                      <  2:394302/151‑1 (MQ=255)                                                                                                                       gAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAAAGCC‑AAC‑TTTAAGGGGGTGGGGTTAAGCCTGGTAT‑GA‑TTTAAACTGGCGGAGAAATAAACTCA‑CCATAAATTGTGCAC‑AAATCTTGCCCATAATTGGAAACCCgaagaaaag                       >  1:1085190/1‑145 (MQ=255)                                                                                                                       gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                    <  2:1061144/151‑1 (MQ=255)                                                                                                                       gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                    <  1:826471/151‑1 (MQ=255)                                                                                                                          cACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACt                 <  2:947658/151‑1 (MQ=255)                                                                                                                                   aGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGt        <  2:230969/151‑1 (MQ=37)                                                                                                                                         ggggATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGcc  >  1:457172/1‑151 (MQ=33)                                                                                                                                                      |                                                                                                                                                CGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCC  >  NC_000913/15242‑15528

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