Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,174,626	0	G	A	41.1%	83.2 / 32.7	74	F180F (TTC→TTT)	ycfT	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (29/14);  new base A (0/30);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 5.06e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACC  >  NC_000913/1174478‑1174775                                                                                                                                                        |                                                                                                                                                        gAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaa                                                                                                                                                       <  2:851539/151‑1 (MQ=255) gAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGa                                                                                                                                                        >  2:612394/1‑150 (MQ=255)  aaCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaa                                                                                                                                                       >  2:874165/1‑150 (MQ=255)   aCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCACACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑t                                                                                                                                                    <  2:787140/151‑1 (MQ=255)         aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCg                                                                                                                                             <  2:125094/151‑1 (MQ=255)          gCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGc                                                                                                                                            >  2:9159/1‑151 (MQ=255)            aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGcc                                                                                                                                           >  1:1117284/1‑150 (MQ=255)                   ccaccaGCAGATGGCGGCGCAACGGCCCCTCTTAAACACAGGTCATTATTGTTGCGCCAAACAAT‑GCGCAAAGGCTGTAAAAAGGCAAATG‑GC‑GGATCACACTGTTCTTCCCCCCCCACGGCGTGGAAACGAAA‑TTAA‑CCGCCACACTca                                                                                                                                    <  1:190797/150‑1 (MQ=255)                     aTCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGCCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTcagc                                                                                                                                  >  2:111812/1‑150 (MQ=255)                        agcggATGGCGGCGCACCGGCCCCTCTTAAACACAGATCTTTATTTTTCCGCAAAACCAT‑GCGCAAAGGCTGTAAAAAGGAAATTT‑GG‑GGATCCCCCTGTTCATCCCCCCCCGCGGGGGGGAAACAAAA‑TAAA‑CCGCCACACTCAGCAGTa                                                                                                                              <  1:546903/151‑1 (MQ=255)                         gcagATGGCGGCGCAACGGCCCCTCTTTAACACAGGTTATTATTGTTGCGCAAAACCAG‑GCGCCAAGGTTGTAAAAAGCCAATTG‑GC‑GGATCACACTTTTCATCCCCCCCCACGGCGTGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTa                                                                                                                              <  2:743101/150‑1 (MQ=255)                               ggcggcGAAACGCCACCTTTTTAACACAGGTAATTATTGTGGCGCCAAACCAT‑GCGCAAAGGCTGAAAAAAGGCAAATG‑GC‑GGATCACATTGTTCATCCCCCCCCACGGCGTGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTa                                                                                                                              <  2:1057091/144‑1 (MQ=255)                                     gcAAGGGCACTTTTTTAACAGAGGTTATTATTGTGGC‑CCAAACCCTGGCGCCACGTATGTAA‑AAGGGAAATTGGG‑GGATCCCATTGTTTATCCCCCCCCGCGGCG‑GGGAAGTAAATTTAC‑CCGCCACACTCAGCAGTACAAACAAGGCg                                                                                                                   <  2:110933/149‑1 (MQ=255)                                       aaCGCCCCCTTTTAAACACAGGTTATTATTTTTCCCCCAACCCAT‑GCCCAAGGGGTGTAAAAAGAAAAATT‑GG‑GGATCCCCCTTTTCATCCCCCCCCACGGGGGGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATa                                                                                                               <  1:253175/151‑1 (MQ=255)                                        aCGGCACCTCTTAAACACAGGTCATTTTTGTTGCGCCAACCAGT‑GCGCCAAGGCTAAAAGAAGGAAAAGT‑GC‑GGATCACCCTGTTATTCCCCCCCCACGGCG‑GGGAAAGAAA‑ATAAGCCGCCACACTCAGCAGTACAAACAAGGCGAATa                                                                                                               <  1:838878/150‑1 (MQ=255)                                         cGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGt                                                                                                             >  2:807476/1‑151 (MQ=255)                                            cACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                                          >  2:627421/1‑151 (MQ=255)                                            cACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                                          >  1:1099891/1‑151 (MQ=255)                                              cctcttaaaacaaGGTTATTATTGTGGCGCAAAACAAG‑GCGCCAGGGTTGAAAAAAGGAAAATG‑GC‑GAACCACTCTGTTCATCCCCCCCCACGGGGGGGGAAAAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                                                         <  1:197845/138‑1 (MQ=255)                                              ccTTTTAAACACGGGTCATTATTGTTGCGCAAACCCAG‑GCGCCAGGGCTTAAAAAAGGAAAATT‑GG‑GAACCACCCTTTTCTTCCCCCCCCACGGCGGGGAAAAGAAA‑TAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                                                         <  1:1181736/150‑1 (MQ=255)                                              ccTCTTAAACCAAGTTAATTATTGTGGCGCCAACCCAG‑GCGCAAAGGCTGTAATAAGGAAAATG‑GC‑GGATCACATTGTTCATCCCCCCCCACGGCGTGGGAAAAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGt                                                                                                             <  2:106968/146‑1 (MQ=255)                                                     aaCACGGTTCTTTTTTTTTCCGCAAAACCTG‑GCGCCAGGGTTATAATAAGGAAATTT‑GC‑TGATCACATTTTTCATCCCCCCCCCGGGGGTGGGAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCa                                                                                                   <  2:711308/149‑1 (MQ=255)                                                     aaCACGGCTCATTATGTTTGCGAAAACCCCG‑GCGCCAAGGCTGTAA‑AAGGGAAATT‑ACAGGATCCTGCTGTCGTTACCCCACCAGGGGGTGGAGACAAAA‑TTAG‑CCGCCACCCTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGACAGa                                                                                                 <  1:334975/151‑1 (MQ=255)                                                     aaCACAGGTCATTTTTGTTGCGCAAAACCAT‑GCGCCAAGGCTGAAAAAAGGAAAATT‑GC‑GGATCACCCTTTTCATCCCCCCCCCCGGCGGGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGa                                                                                                 <  2:858452/151‑1 (MQ=255)                                                     aaCAAAGTTCATTTTTTTGGCGCAAACCAAG‑GCCCCAGGGTTTAAAAAAGGAAAATT‑GC‑GAATCACTCTTTTCATCCCCCCCCCCGGGGTGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAg                                                                                                  <  2:710733/150‑1 (MQ=255)                                                      acacAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGAGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGa                                                                                                 >  1:548011/1‑150 (MQ=255)                                                               attattGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                                                        >  1:453756/1‑150 (MQ=255)                                                                tttttgtggCCCAAAACCAG‑GCGCCAAGGCTGTAATAAGAAAAATT‑GG‑GGATCCCCCTTTTCATCCCCCACCACGGCGGGGGAACAAAA‑TTAC‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                                            <  1:915627/142‑1 (MQ=255)                                                                   tttttgCGCAAACCCTT‑GCGCAAGGGCTTAAATAAGGCAATTG‑GC‑GATCCCCATTGTTCTTCCCCCCCCAGGGCGGGGGAACAAAT‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                                                      <  2:390192/145‑1 (MQ=255)                                                                   ttgttgCCCAAAACCAT‑GCCCCAGGGCTTAAAAAAGGAAATTT‑GG‑GGATCACCCTTTTCATCCCCCCCCCCGGCGAGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                                                     <  2:379405/149‑1 (MQ=255)                                                                      tgccgacAACCCAT‑CCCCCAGGGCTTAAATAAGAAAAATG‑GC‑GCATCACACTGTTCATCCCCCCCCCCGGGGGGAAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaa                                                                                         <  1:176822/136‑1 (MQ=255)                                                                          gcCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                                            >  2:448327/1‑151 (MQ=255)                                                                           ccAACCCAT‑GCGCCAGGGCTGTAATAGGGCAAATT‑GC‑GGATCACACTTTTCTTCCCCCCCCACGGGGTGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCa                                                                           <  2:586179/151‑1 (MQ=255)                                                                             aaaCCAT‑GCCCCAGGGTTGAAAAAAGGCAAATG‑GC‑GGACCACCCTGTTCCTCCCCCCCCCCGGCGGGGGAACAAAA‑TAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAc                                                                          <  1:158387/150‑1 (MQ=255)                                                                             aaaCCAT‑GCCCCAAGGATCAAATAGGAAAATTT‑GC‑GTATCCCTCTGTTCATCCCCCCCCGCGGCGTGGGAAAGAAA‑GTAT‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                                         <  2:408458/151‑1 (MQ=255)                                                                               aCCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGa                                                                        >  2:408346/1‑150 (MQ=255)                                                                                ccAG‑GCCCCAGGGCTTAAATAGGGAAATTT‑GC‑GGACCCCACTTTTCATCCCCCACCAGGGCGGGGAAACAAAA‑TAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCa                                                                           <  1:878515/146‑1 (MQ=255)                                                                                 cAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                                                      >  1:1188765/1‑150 (MQ=255)                                                                                   aggcgcAAGGGCTGTAAAAAGGAAAATG‑GC‑GGATCCCACTGTTCTCCCCCCACCACGGGGTGGGAAAAAAA‑TAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                                                     <  1:136901/148‑1 (MQ=255)                                                                                        caaaGGCTGAAATAAGGCAAATG‑GC‑GGATCACACTGTTCATCCCCCACCCGGGGGGGGGAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACa                                                                <  2:418701/148‑1 (MQ=255)                                                                                         cAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAt                                                              >  1:663800/1‑151 (MQ=255)                                                                                          aaGGCTGAAAAAAGCAAAATT‑GC‑GAACCACACTTTTCATCCCCCCCCACGGCGTGGAAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGATGGCTAAAAATTTTACATACCACGAAATAGACaa                                                               <  2:1063722/149‑1 (MQ=255)                                                                                            ggCTGTACTAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCGGTACAAACAAGGCGAACAGTGGCAGGGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTa                                                            >  1:1089325/1‑150 (MQ=255)                                                                                            gcttgaaaaagggaaaaTG‑GC‑GGATCCCCCTGTTCTTCCCCCCCCCCGGGGGGGGAAGAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                                           <  2:166352/137‑1 (MQ=255)                                                                                                  aataaGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                                      >  2:502390/1‑150 (MQ=255)                                                                                                  aataaGGAAAATT‑GC‑GGATCACACTTTTCTTCCCCCCCCACGGCGGGGAAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCa                                                       <  2:845147/149‑1 (MQ=255)                                                                                                     aaGGAAAATT‑GC‑GGATCACCCTTTTCATCCCCCCCCCGGGCGTGGGAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAg                                                  <  2:191703/151‑1 (MQ=255)                                                                                                      aGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                    >  1:138380/1‑148 (MQ=255)                                                                                                       ggaaaaTT‑GC‑GGATCACACTTTTCATCCCCCACCACGGCGGGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGa                                                 <  2:1075532/147‑1 (MQ=255)                                                                                                       ggaaaaTG‑GC‑GGTTCACACTGTCCATCCCCCACCCCGGCGGGGAAACAAAA‑TTAC‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGa                                                 <  2:3147/147‑1 (MQ=255)                                                                                                        gCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                                                >  1:1075129/1‑150 (MQ=255)                                                                                                        gCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                                                >  2:576092/1‑150 (MQ=255)                                                                                                          aaaTT‑GC‑GGCTCACCCTGTTCATCCCCCCCCAGGGCGTGGGAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAc                                              <  1:688309/150‑1 (MQ=255)                                                                                                                 gcgaataccacTTTTCTTCCCCCCCCCGGGGGTGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa                                         <  2:1193899/142‑1 (MQ=255)                                                                                                                    gATCACCCTTTTCATCCCCCCCCACGGCGTGGAAAAAAAA‑TAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg                                     <  1:425493/151‑1 (MQ=255)                                                                                                                       cccctttttcttcccccccCACGGGGGGGAAAAAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTg                                  <  1:924983/138‑1 (MQ=255)                                                                                                                       cacaCTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTg                                  >  2:861342/1‑151 (MQ=255)                                                                                                                            tGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTgg                              >  1:608583/1‑150 (MQ=255)                                                                                                                              ttCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCgg                           >  1:144957/1‑151 (MQ=255)                                                                                                                              ttCATCCCCCACCACGGGGGGGAAACAAAA‑TAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCg                            <  1:996765/150‑1 (MQ=255)                                                                                                                                cATCCCCCCCCACGGGGGGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCg                            <  2:152071/148‑1 (MQ=255)                                                                                                                                cATCCCCCACCGCGGGGGGGAAACGAAT‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt                          <  2:516496/150‑1 (MQ=255)                                                                                                                                    ccccccCACGGCGGGGGAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatc                      <  2:1012726/150‑1 (MQ=255)                                                                                                                                    ccccACCGCGGTGTGGAAACAATT‑TTAT‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGat                       <  1:275025/149‑1 (MQ=255)                                                                                                                                     ccccccAGGGGGGGGGAAAGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatcat                    <  1:17320/151‑1 (MQ=255)                                                                                                                                     cccACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatca                     >  2:666295/1‑150 (MQ=255)                                                                                                                                     cccACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatc                      >  2:944608/1‑149 (MQ=255)                                                                                                                                       caccacGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATcc                  >  2:1183920/1‑151 (MQ=255)                                                                                                                                         ccacGGCGGGGAAACAAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATcc                  <  1:863293/149‑1 (MQ=255)                                                                                                                                          cacGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGGGATCATccc                 >  1:1083380/1‑149 (MQ=255)                                                                                                                                           acGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTAAATACCACGAAAAAGACAATTAACGCATACAGATACCACAAGCTGGGGCTGGCGGCGATTATCCCg                >  2:828597/1‑149 (MQ=255)                                                                                                                                                gTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATTCCGTGCAgg          >  1:730061/1‑150 (MQ=255)                                                                                                                                                gTGGAAACGAAA‑TAAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAgg          <  2:623771/150‑1 (MQ=255)                                                                                                                                                       cGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCAcc  >  2:720005/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        GAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATT‑GC‑GGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAA‑TTAA‑CCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACC  >  NC_000913/1174478‑1174775

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