Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,410	0	T	G	41.2%	118.2 / 11.5	69	intergenic (+112/‑35)	dnaJ/insL1	chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base T (6/34);  new base G (28/0);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 2.70e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCT  >  NC_000913/15260‑15534                                                                                                                                                        |                                                                                                                              gATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGt                                                                                                                               <  1:363508/151‑1 (MQ=255)  aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGtt                                                                                                                              <  1:16609/151‑1 (MQ=255)    ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAGCGC‑TTACTTAAGGGGtg                                                                                                                             >  2:191042/1‑150 (MQ=255)    ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTtg                                                                                                                             <  2:1131275/150‑1 (MQ=255)      tgtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACTCCCCTAAACGC‑CAACTTTAGGGGTGTgg                                                                                                                          >  2:875073/1‑151 (MQ=255)        tgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCAATCCCCCAAACGC‑CAACTTTAGGGGTGTgg                                                                                                                          >  1:565737/1‑149 (MQ=255)         gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCATAAGCGC‑TAACTTAAGGGGTGTGGGAt                                                                                                                       >  2:1083183/1‑151 (MQ=255)           agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTa                                                                                                                     <  2:668568/151‑1 (MQ=255)            gaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTAc                                                                                                                    <  2:1191421/151‑1 (MQ=255)             aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACg                                                                                                                   <  1:988533/151‑1 (MQ=255)              agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGc                                                                                                                  <  2:780044/151‑1 (MQ=255)               gTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCGCCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGcc                                                                                                                 <  2:677166/151‑1 (MQ=255)               gTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGcc                                                                                                                 <  1:365156/151‑1 (MQ=255)                 tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCCCTCCCCCAAACGC‑TTAATTAATGGGTGTGGGAtttagcc                                                                                                                 >  2:391970/1‑144 (MQ=255)                 tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAACGC‑TAACTTTAGGGGTGTggggattagcc                                                                                                                 >  2:835327/1‑140 (MQ=255)                   tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGa                                                                                                              <  2:287714/150‑1 (MQ=255)                   tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAAGGGTGTGGGATTAAGcccg                                                                                                               >  2:783039/1‑147 (MQ=255)                    ttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACTCCCCTAAACGC‑TAAATTAAGGGGTGTGGGATTATGcc                                                                                                                 >  2:700067/1‑146 (MQ=255)                     tGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGata                                                                                                            <  1:1167099/150‑1 (MQ=255)                             tGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACCCCCCACAACCC‑GTACATAAAGGGTGTGGGATTACGCCTgcgttgtttt                                                                                                      >  2:985132/1‑139 (MQ=255)                             tGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAACGC‑TAAATTAAAGGGTGTGGGTTTACGCCCGATCTGAttt                                                                                                      >  1:1111903/1‑148 (MQ=255)                               aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                                                                                                 <  1:503033/151‑1 (MQ=255)                                  cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                                                                                              <  2:1068309/151‑1 (MQ=255)                                   gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCTAAGCGC‑TAACTTAAGGGGTGTGGGATTTCGCCTGATATGATTTTACCTGcc                                                                                              >  2:633593/1‑150 (MQ=255)                                      aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAAAGC‑TAACTTAAAGGTTGGGGGATTACGCCTTATATGATTTTACGTGcc                                                                                              >  2:211453/1‑147 (MQ=255)                                      aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAGAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGatat                                                                                                           <  2:766351/134‑1 (MQ=255)                                       aCCTCCCCAAAAGCTTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                                                          <  1:111045/150‑1 (MQ=255)                                           ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                       <  2:850054/149‑1 (MQ=255)                                            cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc                                                                                    <  2:82021/151‑1 (MQ=255)                                               aaaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACCCCCACAAACGC‑CAACTTAAAGGGTGTGGGATTACGCCTGCTATGATTTTACATGCCGATGaa                                                                                        >  1:263294/1‑144 (MQ=255)                                                 aaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa                                                                               <  1:972142/151‑1 (MQ=255)                                                   gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCATAAGCGC‑TAAATTAAGGGGTGTGGGATTTCGCCCGATATGATTTAAAGAGCCGGTGAATTAATCTCAc                                                                              >  2:559268/1‑150 (MQ=255)                                                   gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACg                                                                             <  2:440567/151‑1 (MQ=255)                                                      tGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAGCGC‑TAACTTAAGGGGTGTGGGGTTACGCCTGATATGTTTTTACGTGCCGATGAATTACTCTCAc                                                                              >  1:183281/1‑147 (MQ=25)                                                        cccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAc                                                                        <  1:242864/151‑1 (MQ=255)                                                         ccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGC‑CAACTTAAGGGGTGGGGGATTACGCCTGGTATGATTTTAAGTGCCCATGAATTACTCTCACGATaa                                                                         >  2:110819/1‑149 (MQ=255)                                                           gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCATAAGCGC‑TAACTTAAGGGGTGTGGGATTACGCCCGATATGATTTTAAGTGCCGATGAATTACTCTCACGATa                                                                          >  1:296075/1‑146 (MQ=25)                                                             gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACTCCCATAAGCGC‑CAAATTAAGGGGTGTGGTATTACGCCTGATATGATTTTAAGTGCCCATGAATTAATCTCACGATAACTgg                                                                     >  2:1060855/1‑149 (MQ=21)                                                                      tGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAGCGC‑TAACTTAAGGGGTGTGGGATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAATGGTCAGC‑AATTCt                                                         >  2:190532/1‑151 (MQ=255)                                                                      tGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGC‑TAACTTAAAGGGTGTGGGATTACGCCTGATATGATTTAACGTGCCGATGAATTAATCTCACGATAAATGGTCAGC‑AATTCt                                                         >  2:109546/1‑151 (MQ=21)                                                                         gTAAAAATAGGGGGCGGTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTc                                                          <  1:556476/147‑1 (MQ=255)                                                                            aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAGCGC‑TAAATTAAGGGGTGTGGGGTTACGCCTGATATGATTTAAAGTGCCGGTGAATTAATCTCCCGCTAAAAGGTCAGC‑AATTCTGGccc                                                    >  1:589794/1‑150 (MQ=255)                                                                            aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCa                                                   <  1:450242/151‑1 (MQ=255)                                                                              aaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGGACTCCCCCAAGCGC‑TAAATTTAGGGGTTTGGGATTATGCCTGGTATTATTTTACGTGCCGATGAAATACTCTCACGATAACAGGGCAGC‑AATTCTGGcccct                                                  >  1:140064/1‑148 (MQ=255)                                                                               aTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGC‑TAACTTTAGGGGTGTGGGGTTACGCCTGATATAATTTTAAGTGGCGATGAAATACTCTCCCGATAAATGGGCCGCAAATTTTGGcc                                                     >  2:166414/1‑147 (MQ=255)                                                                                 aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTg                                              <  2:1052317/151‑1 (MQ=255)                                                                                  gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGCCCCCCCCCAAGCGG‑TGACCTTTGGGCGGTGGTGGTATGCCTGGTATGATTTTACGGGGCGGTGAATTAGTCTCTCGACAACCGGTCATC‑GACTCTTGGCCATCtt                                               >  2:183976/1‑149 (MQ=255)                                                                                      gCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCTAAACGC‑TAAATTTAGGGGTGTGGGGTTAAGCCTGATATGATTTAACATGCCCCCGAATTAATCTCACGCAAAATAGTCAGC‑AATTTTGGCCCCTATTTGt                                            >  2:1025849/1‑148 (MQ=255)                                                                                      gCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTaa                                          <  1:510724/150‑1 (MQ=255)                                                                                      gCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAg                                         <  2:858058/151‑1 (MQ=255)                                                                                          tGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGGAACCCCCTAAACGCGTAAC‑TAAGGGGGGTGGGTTTTCGCCTGATATGATTTTACCTGCCGTTGAAATTCTCTTCCCAAAAATGGATAGC‑AACTCTTGcccccaattgttaaaccc                                     >  2:286667/1‑136 (MQ=255)                                                                                            aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTCCTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCgaa                                   <  2:679821/151‑1 (MQ=255)                                                                                             aGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAGCGC‑TAAATTAAGGGGTGTGGGGTTACGGCTGGTATGATTTTAAATGCCGCTGAAATACTCTCCCGATAAAAGGGCAGC‑AATTTTGGCCCCTATTTGGAAGccc                                      >  1:246266/1‑147 (MQ=255)                                                                                                 atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACt                              <  1:1159812/151‑1 (MQ=255)                                                                                                   gCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGC‑TAACTTAAGGGGTGTGGGATTACGCCTGATATGATTTAACGTGCCGGTGAATTAATCTCACCATAACTTGGCAGC‑AATTTTTGCCCCTATTGGTAAACCCCAAGAAATgg                            >  2:808545/1‑151 (MQ=17)                                                                                                   gCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAGCGC‑TAAATTAAAGGGTGTGGGATTACGCCTGGTATGATTTTAAGTGCCGATGAATTAATCTCACGCTAAATGGTCAGC‑AAATTTTGCCCATATTGGGAAACCCCAAGAACAgg                            >  2:906075/1‑151 (MQ=255)                                                                                                    cGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGC‑CCACTTTAAGGGGGTGGTGTTACGCCCGGTATGATTTTAAGTGCCGCTGAATTAATCTCCCGATAAATTGTCAGC‑AATTCTTGCCCCCATTTGTaaaaccccagaa                                >  1:495822/1‑136 (MQ=255)                                                                                                             gTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCg                    <  2:1194541/149‑1 (MQ=40)                                                                                                             gTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGc                  <  1:774056/151‑1 (MQ=40)                                                                                                                aaGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACg               <  2:816546/151‑1 (MQ=38)                                                                                                                  gTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa             <  1:1146385/151‑1 (MQ=37)                                                                                                                   tGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa             >  2:75615/1‑150 (MQ=35)                                                                                                                    ggCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTaa            <  1:448089/150‑1 (MQ=35)                                                                                                                    ggCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt           <  2:732317/151‑1 (MQ=35)                                                                                                                       ggggATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGcc        >  2:800185/1‑151 (MQ=33)                                                                                                                        gggATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg       <  1:384018/151‑1 (MQ=33)                                                                                                                        gggATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg       >  2:734393/1‑151 (MQ=33)                                                                                                                          gATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCggg     <  2:611493/151‑1 (MQ=33)                                                                                                                            tCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGc   >  2:342173/1‑151 (MQ=32)                                                                                                                             cACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct  >  2:841486/1‑151 (MQ=31)                                                                                                                              aCTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGTTATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATAGTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct  <  1:571849/150‑1 (MQ=14)                                                                                                                                                        |                                                                                                                              GATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGC‑TAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC‑AATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCT  >  NC_000913/15260‑15534

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