Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,271	0	T	G	48.5%	94.1 / 16.3	68	G165G (GGT→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base T (6/29);  new base G (33/0);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.21e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAT  >  NC_000913/1157122‑1157418                                                                                                                                                       |                                                                                                                                                       catATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGt                                                                                                                                                        <  1:853606/150‑1 (MQ=255) catATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGt                                                                                                                                                        <  2:681592/150‑1 (MQ=255)  atatCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTaa                                                                                                                                                      <  1:1192437/151‑1 (MQ=255)      ccAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGaaa                                                                                                                                                     >  2:1298160/1‑148 (MQ=255)       cAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGaaaa                                                                                                                                                    >  1:754085/1‑147 (MQ=255)          aTTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTgg                                                                                                                                              <  1:827864/151‑1 (MQ=255)          aTTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTgg                                                                                                                                              <  2:193554/151‑1 (MQ=255)            tGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACttg                                                                                                                                              >  1:1067308/1‑147 (MQ=255)              gTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAATTAATGGGtc‑t                                                                                                                                         >  2:296801/1‑151 (MQ=255)                   gaactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGa                                                                                                                                     <  1:524332/150‑1 (MQ=255)                   gaactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATTCTGGATC‑TTgg                                                                                                                                      >  1:469294/1‑149 (MQ=255)                    aactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATATTGGATC‑TTGGCtt                                                                                                                                   >  1:191547/1‑151 (MQ=255)                     actgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGAttt                                                                                                                                  <  2:516610/151‑1 (MQ=255)                        gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑a                                                                                                                               <  2:777728/150‑1 (MQ=255)                        gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑Ac                                                                                                                              <  1:301645/151‑1 (MQ=255)                        gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATCGGTA‑TCGGGtttt                                                                                                                                 >  2:240947/1‑149 (MQ=255)                         aacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATAATGGATC‑TCGGGtttt                                                                                                                                 >  2:1302054/1‑148 (MQ=255)                             aGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATAATGGGTA‑TCGGAgttt                                                                                                                                 >  2:210468/1‑144 (MQ=255)                               ccGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCt                                                                                                                       <  2:1038489/151‑1 (MQ=255)                                cGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTTGATC‑TCGGGTTTT‑TTATTTCCtt                                                                                                                      >  1:1275334/1‑151 (MQ=255)                                     aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTGGGTC‑TCGGGTTTT‑TAATCTCCtttttc                                                                                                                  >  2:1234052/1‑148 (MQ=255)                                         tCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCa                                                                                                              <  1:954967/150‑1 (MQ=255)                                              aCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAAGTAATTGATA‑TCGGGTTTT‑ATATCTCCTTTTTCCGCAtttt                                                                                                          >  1:1057939/1‑147 (MQ=255)                                              aCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATAATTGGTC‑TCGGGTTTT‑TTATCTCCttttttgcatt                                                                                                             >  1:53545/1‑139 (MQ=255)                                                  gTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATGGATC‑TCGGGTTTT‑TTATCTCCTTTTCCCGCATTGTG‑ACCtt                                                                                                   >  1:443323/1‑151 (MQ=255)                                                     aCGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCg                                                                                                <  1:858649/151‑1 (MQ=255)                                                          cGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTTATTGATC‑TCGGGTTTT‑ATAACTCCTTTTTCGGCCTTGTG‑ACCCTCCGTaa                                                                                             >  1:522244/1‑149 (MQ=255)                                                                gATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGTAATTACCGGGTG‑TCGGGTTTT‑ATATCTTCTTTTTCGGGCCTGTG‑GCCtttcctgaagcgg                                                                                        >  2:306576/1‑137 (MQ=255)                                                                   aCGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATAATGGGTC‑TCGGGTTTT‑AAATTTCCCTTTTCCGCATTGTG‑ACCTTCCGTAATTCGGAGCaa                                                                                   >  2:1299100/1‑150 (MQ=255)                                                                    cGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACt                                                                                 <  2:185866/151‑1 (MQ=255)                                                                      cTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTGGGTC‑TCGGGTTTT‑TAATTTCCTTTTTCGGCCTTGTG‑ACCTTCCCTCATGCGGGGCaa                                                                                   >  1:1188542/1‑147 (MQ=255)                                                                        ggCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTgcg                                                                              <  1:1196628/150‑1 (MQ=255)                                                                        ggCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTgc                                                                               <  2:210390/149‑1 (MQ=255)                                                                          cAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAAATACTGGGTC‑TCGGATTTT‑ACATCTCCTTTTTCGGCCTTGTG‑ACCTTCCGTAATTCGGAGCAAATgcgc                                                                             >  1:482706/1‑149 (MQ=255)                                                                           aaTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTTGGTC‑TTGGATTTT‑AAATCTTCTTTTTCGGCATTTTG‑AACTTCCGTAATGCGGGGGAACCGCGCGAt                                                                          >  2:1039650/1‑151 (MQ=255)                                                                              tCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGAtgc                                                                        <  1:406913/150‑1 (MQ=255)                                                                              tCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATTGATCTTCGGATTTT‑‑TATCTCCTTTTTTGGCCTTTTG‑ACCTTCCTGAATTCGGAGCAACTgcg                                                                              >  2:953785/1‑144 (MQ=255)                                                                                tGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATAATGGGTC‑TCGGGTTTT‑AAATCTCCTTTTCCCGCATTGTG‑ACCTTCCGTAATGCGGAGGAAATGCGCGAt                                                                          >  2:327558/1‑146 (MQ=255)                                                                                 gcgcGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATAATGGGTC‑TCGGATTTT‑ACATCTCCTTTTTCGGCATTTTG‑AACTTCCGTAATGCGGAGGAAATGCGCGa                                                                           >  2:1052981/1‑144 (MQ=255)                                                                                          aaaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAAATTCTTGGTC‑TTGGGTTTT‑TTATTTTCCTTTTCGGCCTTTTT‑ACCTTCCGTAAATCGGGGGAACTTCGCGATGCCGCGCGCTAt                                                              >  2:598056/1‑148 (MQ=255)                                                                                          aaaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTGGGTC‑TCGGGGTTT‑AAATTTCCTTTTTCCGCATTGTG‑ACCTTTCCTAATGCGGAGCAACTGCGCGATGCTGCGCGtttt                                                              >  2:775735/1‑146 (MQ=255)                                                                                              gTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTccccc                                                       <  1:1162503/151‑1 (MQ=255)                                                                                              gTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTccccc                                                       <  2:834667/151‑1 (MQ=255)                                                                                               tGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTATTAGGTC‑TCGGGTTTT‑AAATAACCCTTTTCGGCCGTGTG‑GCCCTCCCTAAATCGGAGGGACAGCGCCATGCTGCGCGCTTTGTtcccc                                                       >  1:1266401/1‑150 (MQ=255)                                                                                                 aCGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTgg                                                    <  2:597951/151‑1 (MQ=255)                                                                                                 aCGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTgg                                                    <  1:715780/151‑1 (MQ=255)                                                                                                 aCGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGGTC‑TCGGGTTTT‑AAATCTCCTTTTTCGGCCTTGTG‑AGCTTCCGGTATGAGGGGGAACTGCGCGATGCTGCGCGTTTTGTtcccc                                                       >  1:68979/1‑148 (MQ=255)                                                                                                    gATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTTGGTT‑TCGGGTTTT‑AAAACTCCTTTTTCCGCCATGTG‑AACCTCCGTAATGCGGGGCAAATCCCCGAGGCTGCGCGGTTTGTTCCCCCGGATc                                                 >  1:1133624/1‑151 (MQ=255)                                                                                                           gTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACt                                          <  1:5257/151‑1 (MQ=255)                                                                                                              gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACtggt                                       <  1:155768/151‑1 (MQ=255)                                                                                                               cGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGTAAATAATTGG‑TC‑TTCGATTTT‑TTATCTCCTTTTTCGGCAATGTGAACCTTCC‑TTATACGGGGCAAATGCGCGGTGCCGCGCGGTATTTTCCCCCCGATCgg                                               >  2:837447/1‑142 (MQ=255)                                                                                                               cGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTAGATC‑TCCGGTTTT‑TAAACTCCTTTTTCGGGATTGTG‑AACCTCCGTCAAGCGGGGCCAACGCGCGATGCCGCGGGGTATGTCCCCCCCGGTCGGTTTCttg                                        >  2:859277/1‑147 (MQ=255)                                                                                                                gTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACtggtgg                                     <  1:777253/151‑1 (MQ=255)                                                                                                                gTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTAATTGGTC‑TCGGAGTTT‑ATAACTCCTTTTTCTGGATTTTG‑AACTTTCGTAATGAGGGGGAAATGCGCGGTGCTGCGCGGTATGTTCCCCCCGATTGGGTACtggggg                                     >  1:999797/1‑151 (MQ=255)                                                                                                                gTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATTGATC‑TCGGATTTT‑ATATCTCCTTTTTCGGCATTGTG‑ACCTTCCCTAATTCGGGGCAACCGCGCGATGCTGCGCGTTATGTCCCCCTCGGTCGGGTACtggtgg                                     >  1:1228424/1‑151 (MQ=255)                                                                                                                 tactacACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTTGATT‑TCGGGTTTTAAAACCTCCTTTT‑CGGCCATGTG‑GCCTTTCGTAATGCGGGGCAACTGCGCGATGCTGCGCGTTTTGTCCCCCCGGGTCgggttccggggg                                     >  1:1101792/1‑140 (MQ=255)                                                                                                                    tacACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATGGATC‑TCGGGTTTT‑AAAACTCCTTTTTCCGCATTGTG‑GCCTTCCGTAATGCGGAGGAACTGCGCGATGCTGCGCGTTATTGCCCCCCGGGTCGGTTTCTGGGGGaaa                                  >  1:751492/1‑150 (MQ=255)                                                                                                                     acacTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGaaa                                  <  2:443476/149‑1 (MQ=255)                                                                                                                     acacTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACt                                <  2:901668/151‑1 (MQ=255)                                                                                                                     acacTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGGAC‑TCGGGTTTT‑TTATCTCCTTTTTCGGCATTGTG‑ACCTTCCGCAATGGGGGGCAACTGCGCGATGCTGCGCGTTTTGTCCCCCCGGATCGGTTACTGGGGGgaa                                  >  2:1318384/1‑149 (MQ=255)                                                                                                                       acTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGATC‑TCGGGTTTT‑AAAACTCCTTTTTTGGCCTTGGG‑ACATTCCGTAAAGCGGGGGAAATGCGCGATGCTGCGCGGTTTGTTCCCCTGGATCGGGTACTGGTGGaaaat                                >  1:1032709/1‑147 (MQ=255)                                                                                                                       acTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATTGGAG‑TTGGGTTTT‑AAATCTCCTTTTTCCGCGTTGTG‑GGCTTCCCTTAATCGGGGCGACAGCGCGATGCTGCGCGGTATGTTCCCCCGGCTCgggttcttgggg                                     >  1:56125/1‑134 (MQ=255)                                                                                                                       acTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTGGGTC‑TCGGGTTTT‑TAATCTCCCTTTTCGGCATTGTG‑ACCTTTCGTAATGCGGAGCAAAAGCGCGATGCTGCGCGGTTTGTCCCCCCGGATCGGGTACTGGTGGAAAc                                 >  2:635109/1‑148 (MQ=255)                                                                                                                          gTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACt                            <  1:296798/150‑1 (MQ=255)                                                                                                                          gTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACt                            <  2:975349/150‑1 (MQ=255)                                                                                                                                     gACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTg                <  2:314859/151‑1 (MQ=255)                                                                                                                                      aCAGAGAAACGGC‑GGGGAAATTAATGGAGT‑TCGGGTTTT‑TAATCTCCTTTTTCGGCCTTGTG‑GCCTTCCGTAATGCGGGGCAACTGCGCGATGCTGCGCGGTTTGTCCCCCCGGGTTGGGTTCTGGGGGGAAATGACTCaccttacctt                 >  1:367383/1‑149 (MQ=255)                                                                                                                                            aaaCGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGt         <  1:536955/151‑1 (MQ=255)                                                                                                                                                ggc‑ggGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCg     <  1:240944/151‑1 (MQ=255)                                                                                                                                                    gggTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAt  <  2:918358/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATC‑TCGGATTTT‑ACATCTCCTTTTCCGGCATTGTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAT  >  NC_000913/1157122‑1157418

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