Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	869,107	0	A	G	43.1%	81.2 / 15.5	59	I530V (ATT→GTT)	gsiA	glutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (5/28);  new base G (25/0);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 8.35e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGC  >  NC_000913/868960‑869251                                                                                                                                                    |                                                                                                                                                              aaTCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTa                                                                                                                                                            <  2:762626/151‑1 (MQ=255) aaTCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGAGGTAGAGAGGagt                                                                                                                                                             >  2:330552/1‑148 (MQ=255)        aaGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCg                                                                                                                                                     <  1:1171621/151‑1 (MQ=255)            gATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGGGGGGTAGTCATCgtgt                                                                                                                                                  >  1:607019/1‑150 (MQ=255)            gATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCgtgt                                                                                                                                                  <  1:322497/150‑1 (MQ=255)                  tGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGGTAATCATCGCGTTGC‑gg                                                                                                                                            >  1:92570/1‑149 (MQ=255)                    ccGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAt                                                                                                                                        <  1:723151/151‑1 (MQ=255)                    ccGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGTGTGGC‑gggg                                                                                                                                          >  1:48570/1‑147 (MQ=255)                     cgacgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑g                                                                                                                                     <  2:168162/151‑1 (MQ=255)                      gacgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGTGTGgg‑ggggg                                                                                                                                         >  1:346191/1‑145 (MQ=255)                         gaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTTGTCATCGTGTGGC‑GGTgggt                                                                                                                                       >  1:1132546/1‑144 (MQ=255)                            gCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑g                                                                                                                         <  2:25592/151‑1 (MQ=255)                            gCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑g                                                                                                                         <  2:684340/151‑1 (MQ=255)                             ccGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑gg                                                                                                                        <  2:21603/151‑1 (MQ=255)                               gTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGGGAGGGGGGTAGTTATTGTGTGGC‑GGGGgt‑‑gttgctctg                                                                                                                             >  1:715533/1‑142 (MQ=255)                                  tCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGCCATCGTGTGGC‑ggggggggttttcggg                                                                                                                              >  2:311428/1‑133 (MQ=255)                                   cGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGc                                                                                                                      <  2:1244682/147‑1 (MQ=255)                                   cGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGAGTTATCATCGGGTGGC‑GGGTGT‑‑GT‑A‑TT‑T‑T‑GGGGaaa                                                                                                                   >  1:135129/1‑150 (MQ=255)                                        cTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGGGGGGGGGTTTGTCATCGTGTGgg‑gggggtgggttttggg                                                                                                                              >  2:525932/1‑127 (MQ=255)                                         tGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTTGTCATTGGTTGGG‑GGGGAT‑‑GTGT‑TT‑T‑C‑GGGCAAAtatt                                                                                                               >  2:896067/1‑146 (MQ=255)                                            aTGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTTGTCATCGTGTGGC‑GGGGGT‑‑GT‑T‑TT‑TCC‑GGGGAAAtttt                                                                                                               >  2:891179/1‑144 (MQ=255)                                            aTGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAACGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGTGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGa                                                                                                            <  2:1129873/148‑1 (MQ=255)                                             tGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGC‑GGGGGT‑‑GT‑A‑TT‑TCC‑GGGCAAAtgttttaa                                                                                                           >  1:945030/1‑142 (MQ=255)                                              gTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGATTTGTTATCGTGTGGC‑GGGGGT‑‑GT‑A‑TC‑T‑CGGGGGAAATTTTTGA‑AAt                                                                                                        >  1:1169005/1‑150 (MQ=255)                                                  cTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTTGTCCTCCTGTGgg‑gggggtgggttt                                                                                                                                  >  2:867682/1‑117 (MQ=255)                                                    aTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTc                                                                                                  <  2:8268/148‑1 (MQ=255)                                                    aTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGAGGTAGGGGGGGTTTGTCATCTGGTGGG‑GGGGATAGGT‑A‑TC‑T‑‑‑GGGGAAATTTTTTA‑AAtttttccc                                                                                                 >  2:835303/1‑145 (MQ=255)                                                       cgcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGgc                                                                                            <  1:455423/151‑1 (MQ=255)                                                         cgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTTGTAGAGCGGGTTTGTCATCGTGTGGG‑GGGGGG‑‑GT‑ATTC‑T‑C‑GGGAAAAtttttttatttgtgccc                                                                                                  >  2:96151/1‑134 (MQ=255)                                                         cgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGGGTGGG‑GGGGGG‑‑GG‑A‑TT‑TCC‑GGGGAAATTTTTGA‑AATTgtgccc                                                                                                 >  1:708845/1‑141 (MQ=255)                                                          gggCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGC‑GGGGGT‑‑TT‑A‑TT‑T‑CGGGGCAAATTTTTTA‑AATTgtgccc                                                                                                 >  1:544234/1‑140 (MQ=255)                                                                aTTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTTATCGGGTGGG‑GGGGGT‑‑GT‑A‑TT‑TCC‑GGGGAAATTTTTGA‑AATTG‑TTCCCCg                                                                                              >  1:1150088/1‑141 (MQ=255)                                                                   aTCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTc                                                                              <  1:1002145/150‑1 (MQ=255)                                                                            ttgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGCATAGGCATCGCGTGGC‑GGGGTT‑‑GT‑G‑TG‑T‑CGGGGGAAATGGTTTA‑AA‑TG‑GTCCCCGGCCCC‑C‑gggttttcggaa                                                                          >  1:25780/1‑136 (MQ=255)                                                                             tgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGCGGGGGTT‑‑GT‑T‑TT‑C‑G‑GGGCAAAATTTTGA‑AATTG‑GTCCCCGGGGGG‑G‑GGGT‑TCCGAAAAccc                                                                      >  2:50527/1‑149 (MQ=255)                                                                                tCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTTGAGCGGGTTAGTTATCGTGTTGC‑GGGGGG‑‑GT‑T‑TCTT‑T‑GGGGAAATT‑TTGAGAATTGTGTACAGGGCGCG‑G‑GGTT‑TTCAAAAcccc                                                                      >  2:463219/1‑147 (MQ=255)                                                                                 cGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCgcagca                                                                <  1:697515/150‑1 (MQ=255)                                                                                   aTCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGGAGAGCGGGTTAGTCATCGTTTGGC‑GGGGGG‑‑GT‑A‑TT‑TCC‑GGGGAAATTTTTGA‑AATTGTGCCCCCGGCGCG‑G‑GGTT‑TTGAAAAAcccaccccac                                                               >  2:477980/1‑144 (MQ=255)                                                                                       ccAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCtt                                                           <  1:465824/149‑1 (MQ=255)                                                                                       ccAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCt                                                            <  1:687314/148‑1 (MQ=255)                                                                                       ccAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGC‑GGGGGT‑‑GG‑A‑TC‑TCC‑GGGGAAATTTTTTA‑AATT‑‑GTCCCCCGCGCG‑C‑GGGCTTTCGAAAAaccccgcaatcct                                                            >  1:445930/1‑139 (MQ=255)                                                                                       ccAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTAGAGCGGGTTTGTTATCCTGTGGG‑GGGGGT‑‑GT‑T‑TT‑TCC‑GGGGAAATTTTTGA‑AATT‑‑GTCCACCCCGCG‑CGGGGT‑TTCGAAAAccccccccctc                                                              >  2:1105468/1‑139 (MQ=255)                                                                                        cAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGGGTGGC‑GGGGGG‑‑GG‑A‑TC‑TCC‑GGGCAAATGTTTGA‑AATTG‑GTCCCCGG‑GCG‑C‑GGGCTTTCGAAAAcccccgcaacccttt                                                          >  1:629519/1‑138 (MQ=255)                                                                                        cAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGGGTGGC‑GGGGAG‑‑GT‑A‑TT‑C‑T‑GGGCCAATTTTTTT‑AAATA‑GTCCACCGCCCGCC‑GGCT‑TTGGAAAAcccccgccccc                                                              >  1:1112053/1‑138 (MQ=255)                                                                                               atgtCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTaa                                                 <  1:321200/148‑1 (MQ=255)                                                                                               aTTTCGGCATTGCTTATCTGTTTAACTCCCACGCGATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑G‑GGTC‑TTCGAAAACCCGCAGCTTCCTTATACGCGTa                                                  <  2:127463/150‑1 (MQ=255)                                                                                                   cGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAAtt                                              <  1:730369/150‑1 (MQ=255)                                                                                                    ggCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTAc                                            <  2:729860/151‑1 (MQ=255)                                                                                                         tGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTAGAGAGGGTTAATCATCGTGTGGC‑GGCGGG‑‑TT‑G‑TT‑T‑CGGGGAAAATTTTTTA‑ATTTG‑TTCCAC‑GCGCG‑CGGGCC‑TTCG‑AAACCCCCAGCATTCCtataagacggaaataatt                                          >  2:971989/1‑134 (MQ=255)                                                                                                          gCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTggcgg                                      <  2:798958/151‑1 (MQ=255)                                                                                                          gCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTggcg                                       <  2:603262/150‑1 (MQ=255)                                                                                                                  gTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTcc                               <  2:346176/150‑1 (MQ=255)                                                                                                                     tATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGt                            <  1:910199/150‑1 (MQ=255)                                                                                                                      aTCTCCCACGATATGGCGGTGGGAGAGCGGGGTAGTCATCGTGTGGG‑GGGGGT‑‑GT‑T‑TT‑T‑C‑CGGGGAAATGTTTT‑AAATG‑GTGCACCGCGCG‑C‑GGGG‑TTTGTAAAACCCCCGCAACCTTTTTAGCGGAAAATACTtgcggccggttccg                             >  2:166009/1‑140 (MQ=255)                                                                                                                          cccACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTg                       <  2:771904/150‑1 (MQ=255)                                                                                                                           ccACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGaa                     <  2:1132568/151‑1 (MQ=255)                                                                                                                               gATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGt                 <  1:525953/151‑1 (MQ=255)                                                                                                                                atatGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTc                <  2:204624/151‑1 (MQ=255)                                                                                                                                               agCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGc  >  1:490606/1‑150 (MQ=255)                                                                                                                                                gCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGGAACGTCCCGACAACGAcc    >  1:472275/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                              AATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGC‑GGTGAT‑‑GT‑A‑TC‑T‑C‑GGGCAAATTGTTGA‑AATTG‑GTCCACGGCGCG‑C‑GGTC‑TTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGC  >  NC_000913/868960‑869251

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