Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,395	0	G	A	42.9%	84.0 / 34.2	70	intergenic (+97/‑50)	dnaJ/insL1	chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base G (6/34);  new base A (30/0);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 3.60e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTC  >  NC_000913/15250‑15535                                                                                                                                                  |                                                                                                                                             gagCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAacg                                                                                                                                            >  1:309115/1‑148 (MQ=255) gagCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTa                                                                                                                                         <  1:469466/151‑1 (MQ=255)    cttcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCCaa                                                                                                                                        >  2:713175/1‑149 (MQ=255)      tcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                                                                                                                                      <  2:228616/149‑1 (MQ=255)      tcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                                                                                                                                      <  2:391595/149‑1 (MQ=255)       cttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACACCCCTAAGCGCCAAcatt                                                                                                                                    >  1:991748/1‑147 (MQ=255)        tttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACACCCCTAAACGCTAACttta                                                                                                                                   >  2:785555/1‑148 (MQ=255)         ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAACGCTTACTTaa                                                                                                                                   >  2:1132044/1‑149 (MQ=255)         ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAg                                                                                                                                  <  1:638091/150‑1 (MQ=255)          tGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCATAAACGCTAACTTAAggg                                                                                                                                >  2:340124/1‑151 (MQ=255)            aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGGACTCCCCTAAACGCGAAATTAAggg                                                                                                                                >  1:885589/1‑149 (MQ=255)             tGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGCACTCCCCTAAACGCTAACCTAAagg                                                                                                                                >  2:236277/1‑148 (MQ=255)              ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGTGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgt                                                                                                                            <  2:1134059/151‑1 (MQ=255)              ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACACCCACAAACGCTAACTTAAggggt                                                                                                                              >  2:974526/1‑147 (MQ=255)                      gaagTGTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                                                                                    <  1:341446/151‑1 (MQ=255)                       aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCTCTTCCATAAACGCTAAACTAAAGGGTGTGGTGtt                                                                                                                      >  2:1136891/1‑148 (MQ=255)                       aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                                                                                                   <  2:885392/151‑1 (MQ=255)                        agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACTCCCCCAAAAGCTCACCTAAAGGGTGTggggtttcg                                                                                                                   >  2:973852/1‑143 (MQ=255)                           tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                                                                                                <  1:463585/150‑1 (MQ=255)                                 acgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAACGCCAACATAAggg                                                                                                                                >  1:98542/1‑128 (MQ=255)                                 acgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                                                                                                         <  1:62152/151‑1 (MQ=255)                                   gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCTAAACGCTAACTTTAGGGGTGTGGGATTTAGCCTGatat                                                                                                           >  2:1020610/1‑147 (MQ=255)                                     ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                     <  2:627210/151‑1 (MQ=255)                                      cTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCCAAACGCCAACTTAAAGGGTGTGGGGTTAAGGCCGATATTAttt                                                                                                      >  1:324726/1‑149 (MQ=255)                                         aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                                                                                                 <  2:811806/151‑1 (MQ=255)                                           ccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCTAAATTTAAGGGTGTGGTATTTTGCCTGATTTGATTTTACGt                                                                                                 >  1:561824/1‑149 (MQ=255)                                            cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGAACTCCCCCAAGAGCTAACTTAAGGGGTGTGGGATTTCGCCCTatatttttt                                                                                                      >  1:490434/1‑138 (MQ=255)                                             gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGGGGCGGGGGTGCCTCCCCCAAAAGCAAACATAAATGTTGTGGTAACACTCCCGATATGAATTaaaattacg                                                                                             >  1:708259/1‑144 (MQ=255)                                               tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGAACTCCCCTAAACGCGAAATTAAGGGGTGTGGTGTTTCGCCTGATATGAttt                                                                                                      >  1:267893/1‑140 (MQ=255)                                                aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGTACTCCCCTAAACGCTAAATTAAAGGGTGTGGGATTTCGCCTGCTATGATTTTACATGCCGATg                                                                                          >  1:422301/1‑151 (MQ=255)                                                    tCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAGCGCTCACTTTAGGGGTGTGGGATTATGCCCGATATTATTTTAAGTGCCGGTGaa                                                                                        >  1:320315/1‑149 (MQ=255)                                                     ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                     <  1:345985/151‑1 (MQ=255)                                                      cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt                                                                                      <  1:978468/149‑1 (MQ=255)                                                       ccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                                   <  1:812790/151‑1 (MQ=255)                                                         aaaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct                                                                                 <  2:693180/151‑1 (MQ=255)                                                           aaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa                                                                               <  2:627179/151‑1 (MQ=255)                                                            aGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGAGTCTCCCCTAAGGGGTACCCTAAAGGTTTCGGGGTTAAGCCCGTTATGATTTAACGTACCGATGACTTACtctc                                                                                >  1:253963/1‑149 (MQ=255)                                                            aGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGCACTCCCCTAAAAGCTAAATTTAAGGGTGTGGGATTTAGGCCGATATGATTTTAAGGGTCGATGAAATAAtctc                                                                                >  1:385087/1‑149 (MQ=255)                                                             gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAATACACCCATAAACGGTCAATTTAGGGTTGTGGGAATTCACCCTATATTATTTAAAATGCTCATGAAtat                                                                                     >  2:207345/1‑141 (MQ=255)                                                                 gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCTAAATTTAAGGGTGTGGGGTTTAGCCCGATATGATTTAAAGTGCCGGTGGATTAATCTCCCCATa                                                                          >  2:755803/1‑150 (MQ=255)                                                                 gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAt                                                                           <  2:990264/149‑1 (MQ=255)                                                                 gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATaa                                                                         <  1:434551/151‑1 (MQ=255)                                                                 gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa                                                                          <  1:85762/150‑1 (MQ=255)                                                                  cccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAAc                                                                        <  2:576068/151‑1 (MQ=255)                                                                   ccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCACAAACGCCAAATTAAGGGGTGTGGGGTTTAGCCTGATATGATTTTAAGTGCCGCTGAATTAAtctccccaaaa                                                                         >  2:90517/1‑142 (MQ=255)                                                                    cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCCAAACGCCAACTTTAGGGGTGTGGGATTTCGCCTGATATGATTTTAAGTGGCGATGAAATTATCTCACGATaa                                                                         >  2:762044/1‑148 (MQ=255)                                                                      tGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGCACTCCCCTAAAAGCTAACTTAAGGGGTGTGGGGTTACGCCTGATATGATTTTAAATGCCGATGAATTATTCTCACGATAACTgg                                                                     >  2:776833/1‑150 (MQ=255)                                                                        ggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCCAACTTAAAGGGTGTGGGGTTAACCCCGATATGATTTTAAGTGCCGATGAATTAATCTCACGATaaaatg                                                                     >  2:1030419/1‑144 (MQ=255)                                                                        ggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                                                                  <  2:111551/151‑1 (MQ=255)                                                                         gCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAACACTCCCATCAAGGGAAACCTTAGGGGTGGGGGGTTACAACCGGTAATTTTTAACGGTCCGGTGAAATATTCTCCCGATaaatcggc                                                                   >  1:438737/1‑143 (MQ=255)                                                                                     aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACTCCCCCAAACGCTAACTTAAAGGGTGTGGGGTTTCGCCTGGTATTATTTTACGTGCCGACGAATTAATCTCACGCTAACTTGTTAGCAAATTTTGccc                                                     >  1:523663/1‑151 (MQ=255)                                                                                     aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGAGGGCGCCCCCACAAAAGCTAACTTAAGGGGTGTGGTATTTCGCCTGGAATGATTTTAAGTGCCGATGGATTTTTCTCACCATAACTGGTCAGCAAATTTGGgcc                                                     >  1:286324/1‑151 (MQ=255)                                                                                           aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCCTACTTAAGGGGTGTGGGATTATGCCTGGTATTATTTTACGTGCCGATGGATTAATCTCACGATAACTGGTCAGCAATTCTTGCCCCTAtt                                                >  2:477953/1‑150 (MQ=18)                                                                                           aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                                               <  2:44429/151‑1 (MQ=255)                                                                                            gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAACGCGAACTTCAAGGGTGTGGGATTTAGCCCGGTATGATTTTACCTGCCGAAGAATTATTCCCAACATAAATTGTCAGCAATTCTGGGCCAAAtttg                                              >  1:1066823/1‑149 (MQ=255)                                                                                            gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTgg                                              <  2:394100/151‑1 (MQ=255)                                                                                                gCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTaa                                           <  1:69464/150‑1 (MQ=255)                                                                                                   ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAAAGCTCAATTTAGGGGTGTGGGGTTACGCCCGATATGATTTAAAGTGCCGAGGAAATATTCTTTCCATAAATGGGCAGCAActctttgcccaatttg                                              >  2:409067/1‑132 (MQ=255)                                                                                                      aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaa                                    <  2:1080723/151‑1 (MQ=255)                                                                                                           atGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACt                               <  2:561845/151‑1 (MQ=255)                                                                                                              cGAGCACCTGTAAAGTTGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGa                            <  2:495870/151‑1 (MQ=255)                                                                                                              cGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGa                            <  1:1020574/151‑1 (MQ=255)                                                                                                               gAGCACCTGTAAAGTGGCGGTGTTGACTCCCCACAACGCCAACTTAAGGGTTGTGGGATTTCCCCCGTTATGACTTAACGTGCCCTTGGATTTCTTTTACGATAACTGGTCAGCAATTCTTGCCCCTTTTGGTTAGCCCGAAAAACTg                              >  1:188524/1‑148 (MQ=9)                                                                                                               gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                           <  1:40672/151‑1 (MQ=255)                                                                                                               gAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                           <  1:1030418/151‑1 (MQ=255)                                                                                                                   aCCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACtt                       <  1:336982/151‑1 (MQ=255)                                                                                                                         aaaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCAc                 <  1:124516/151‑1 (MQ=39)                                                                                                                          aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACg                <  1:631263/151‑1 (MQ=38)                                                                                                                            gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa              <  1:963722/151‑1 (MQ=37)                                                                                                                                       cACTCGCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct   >  2:455178/1‑151 (MQ=14)                                                                                                                                        aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGctc  >  1:299948/1‑151 (MQ=31)                                                                                                                                                  |                                                                                                                                             GAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTC  >  NC_000913/15250‑15535

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