Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,395	0	G	A	41.4%	66.1 / 26.2	59	intergenic (+97/‑50)	dnaJ/insL1	chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base G (8/26);  new base A (24/0);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.40e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCC  >  NC_000913/15246‑15528                                                                                                                                                      |                                                                                                                                            cAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAg                                                                                                                                             <  1:572982/150‑1 (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)                                                                                                                                                      |                                                                                                                                            CAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCT‑AAC‑TTAAGGGTTGTGGTATTACGCCTGATAT‑GA‑TTTAACGTGCCGATGAATTACTCTCA‑CGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCC  >  NC_000913/15246‑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: ‑