Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,245,234	0	G	T	44.1%	104.1 / 35.2	68	L25F (TTG→TTT)	ycgY	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (9/29);  new base T (30/0);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGAT  >  NC_000913/1245086‑1245383                                                                                                                                                     |                                                                                                                                                         tCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTg                                                                                                                                                          <  1:896913/149‑1 (MQ=255)    atatCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTg                                                                                                                                                      >  2:1088129/1‑150 (MQ=255)      atCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTgcg                                                                                                                                                    <  2:600775/150‑1 (MQ=255)          gATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTAt                                                                                                                                               <  2:122682/151‑1 (MQ=255)             aaaGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTAtt                                                                                                                                              >  1:954189/1‑149 (MQ=255)                gTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACa                                                                                                                                         <  1:557365/151‑1 (MQ=255)                  cGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTAttttaa                                                                                                                                          >  1:1133733/1‑144 (MQ=255)                       cccGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCGGCGCTTATTGAATAttt                                                                                                                                      >  1:12681/1‑147 (MQ=255)                        ccGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAAt                                                                                                                                 <  1:281659/151‑1 (MQ=255)                           gTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAACATTTAGAAtt                                                                                                                                >  2:970282/1‑149 (MQ=255)                             aacaacaacAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTTAAAATTTTGaa                                                                                                                                  >  1:858237/1‑145 (MQ=255)                              acaacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCa                                                                                                                           <  1:755915/151‑1 (MQ=255)                              acaacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCa                                                                                                                           <  1:685003/151‑1 (MQ=255)                                 acaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTGAACATTTTGAATTGATCaa                                                                                                                          >  2:1081728/1‑149 (MQ=255)                                 acaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCCTTTTGAAAATTTAGAATTTATCaa                                                                                                                          >  2:466722/1‑149 (MQ=255)                                 acaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAAAATTTAGAATTTATTaa                                                                                                                          >  1:1123179/1‑149 (MQ=255)                                  caacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAAAAATTAGAATTGATCAATcc                                                                                                                       >  1:326600/1‑151 (MQ=255)                                   aacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑t                                                                                                                     <  1:318527/151‑1 (MQ=255)                                        aaTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTTAACATTTAGAATTTATCAATCC‑TTCt                                                                                                                  >  1:903637/1‑149 (MQ=255)                                               ttttAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAAAATTTAGAATTTATTAATTC‑TTCTGGGAGAtt                                                                                                          >  1:65776/1‑150 (MQ=255)                                                   aaTTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑c                                                                                                     <  1:880588/150‑1 (MQ=255)                                                    attattGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAACATTTAGAATTGATCAAACC‑TTCTGGCAGATTAAG‑cc                                                                                                    >  1:948212/1‑150 (MQ=255)                                                       attGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAACATTTAGAATTTATTAATCC‑TTCTGGCAGATTAAG‑CCaaa                                                                                                 >  1:585807/1‑150 (MQ=255)                                                          gTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATAc                                                                                             <  2:234808/151‑1 (MQ=255)                                                            ttCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAAAATTTTGAATTTATCAAACC‑TTCTTGGCGATTTAG‑CCaaaaaaac                                                                                            >  1:590577/1‑146 (MQ=255)                                                              cGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGt                                                                                         <  2:533417/151‑1 (MQ=255)                                                              cGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGt                                                                                         <  2:897023/151‑1 (MQ=255)                                                               gttgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAACATTTAGAATTGATTAATCC‑TTCTGGGAGATTAAA‑CCCAAATACCCGtt                                                                                        >  2:1142346/1‑151 (MQ=255)                                                               gttgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGCTTTTTGAAAATTTTGAAATGATTAAACC‑TTTTTGGGGATTAAA‑CCaaaaaaaccctt                                                                                        >  1:607272/1‑143 (MQ=255)                                                                 tgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAAAAATTAGAATTTATCAATAC‑TTCTTGCGGAATTAG‑CCAAAATAacc                                                                                           >  1:915728/1‑146 (MQ=255)                                                                  gttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTTAAATTTATGATTTATTAATCC‑TTCTTGGAGATTAAG‑CCAAAAAAccc                                                                                           >  1:418697/1‑145 (MQ=255)                                                                     tttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCa                                                                                   <  1:648895/150‑1 (MQ=255)                                                                                tCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCCTTTTGAACATTTAGAATTTATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATAACCGTTCGGCAATGCTT‑CAGtt                                                                       >  2:439317/1‑150 (MQ=255)                                                                                 cAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTGAACAATTAGAAATTATCAATCC‑TTCTGGCCGATTTAG‑CCAAAATACCCGTTCGGCAATGCtttcggt                                                                        >  2:850278/1‑144 (MQ=255)                                                                                  aGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTAcg                                                                    <  2:648301/151‑1 (MQ=255)                                                                                  aGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTAcg                                                                    <  2:852547/151‑1 (MQ=255)                                                                                  aGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTAcg                                                                    <  2:931739/151‑1 (MQ=255)                                                                                  aGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGGTTTTGCTTCGCTTATTTAACATTTTAAAATTATCAATCC‑CTCTGGCAGATTAAG‑CCCAAAAAACCGTTCGGGAAAGCTT‑CAGtt                                                                       >  1:841531/1‑148 (MQ=255)                                                                                   gAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTGAACATTTTAAATTGATCAAATC‑TTCTGGCAGATTAAG‑CCAAAATAACCCTTTCGCTATGTTT‑TAGtt                                                                       >  2:403654/1‑147 (MQ=255)                                                                                    aaaaTAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTGAAAATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATAACCGTTTGGCAATGCTT‑CAGTTTcgcg                                                                  >  2:479045/1‑151 (MQ=255)                                                                                     aaataaatATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTGAACATTTAGAATTTATCAATCC‑TTCTGGCAGATTAAG‑ACAAAATACCCCTTCGGCAATGCTT‑CCGTTAcgc                                                                   >  2:279794/1‑149 (MQ=255)                                                                                      aataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGCTTTTTTAACATTTAGAAATGATCAAACC‑TTCTGGCAGATTAAG‑CCAAAATAACCGTTCGGGAATGCTT‑CCGTTaag                                                                    >  1:76880/1‑145 (MQ=255)                                                                                      aataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGATTTGGCTGCGCTTTTTCAACATTTAGAAATGTTTAAACC‑CTCTTGCAGCTTAAG‑CCCAAAGAACCCTTCCGCCATGCTT‑CTGGTAcggg                                                                  >  2:152625/1‑147 (MQ=255)                                                                                            tatTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGGGGTTATTTAAAATTTAGAATTGATGAATAC‑TTCTTGCAGATTAAA‑GCCAAATACCCCTTTGGCAAAGCCT‑TAGTTTCCCGAAGAgg                                                            >  1:154670/1‑149 (MQ=255)                                                                                               tttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATg                                                        <  1:688623/150‑1 (MQ=255)                                                                                                ttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTGTACATTTAGAATTGATTAATCC‑TTCTGGCAGATTAAG‑CCAAAAAACCCCTTCGGCAATGCTT‑CATTTACGCGGAGAgg                                                            >  1:1234549/1‑145 (MQ=255)                                                                                                 tAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGGTTTTGTGGCGCCTATTTAAAATTTAGAAATTATCAATCC‑TTTTAGCCGATTAAG‑CCAAAATAACCGTTTGGCCAAGCTT‑TCGTTTCGCGGagag                                                             >  1:79491/1‑143 (MQ=255)                                                                                                  aGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCt                                                    <  1:970273/151‑1 (MQ=255)                                                                                                    aCAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGa                                                  <  1:1138220/151‑1 (MQ=255)                                                                                                         aGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTTTTGAACATTTTGAATTGGTCAAACC‑TTCTGGCAGATTAAG‑CCAAAATAACCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGGCTGAATTCt                                             >  1:381764/1‑151 (MQ=255)                                                                                                            cAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCtttt                                          <  2:1158169/151‑1 (MQ=255)                                                                                                             aGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCttttt                                         <  1:738629/151‑1 (MQ=255)                                                                                                              gTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGCTTTTGAAGCGCCTTTTGAAAAATTAGAATAGATCAAACCTTTCGGGAACATTACA‑CCAAAAAACCGGTCCGGC‑ATGATTCCAGTT‑CGCGGCGAGGAATGGTCTGAATTTcttttt                                        >  2:395024/1‑151 (MQ=255)                                                                                                                cAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGATTTGGGTGCGCTTTTTGGAAATTTTGAATTTTTCAAAAC‑TTCTTGCAGATTAAG‑CCAAAATAACCGTTCGGCAATGCCT‑CTGTTTAGCGAAGAGGGATGGTTTTAAtttttttttt                                       >  1:937449/1‑150 (MQ=255)                                                                                                                 aGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTGAACATTTAGAATTTATCAATCC‑TTCTTGCAGATTAAG‑CCAAAAAACCCGCTCGGCAATGCTT‑CAGTTACGCGAAGAGGGATGGTCTGAATTCttttttt                                       >  2:579586/1‑149 (MQ=255)                                                                                                                 aGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTgg                                     <  1:133523/151‑1 (MQ=255)                                                                                                                  ggCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGt                                    <  1:1038517/151‑1 (MQ=255)                                                                                                                   gCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGtt                                   <  2:326603/151‑1 (MQ=255)                                                                                                                   gCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGtt                                   <  1:1228400/151‑1 (MQ=255)                                                                                                                          aaaTGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAAt                            <  2:806225/151‑1 (MQ=255)                                                                                                                            aTGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATtc                          <  2:478602/151‑1 (MQ=255)                                                                                                                                ggCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTg                      <  1:319182/151‑1 (MQ=255)                                                                                                                                     aTAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAgtg                 <  2:1016790/151‑1 (MQ=255)                                                                                                                                      tAGTGAGGCGGTTTTGGCTGCGCTTTTTGAACATTTTGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAg                   >  1:1082919/1‑148 (MQ=255)                                                                                                                                       aGTGAGGCGGTTTTTGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAGCCCAAAAAACCCGTTGGGAAATGTTT‑CGGTTCGGGGAGGGGGATGGGTTTAATTTTTTTTTTTGGTTCTTAAAt                            >  1:19802/1‑139 (MQ=255)                                                                                                                                       aGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAgtgt                <  2:12681/150‑1 (MQ=255)                                                                                                                                          gAGGCGGTTTTTGCTGCGCTTTTTGAACATTTAGAATTGATCAATCC‑TTCTTGCAGATTAAG‑CCCAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTTTCTGGAGTGTCTgg            >  1:1190154/1‑151 (MQ=255)                                                                                                                                             gCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGAtt         >  2:1141699/1‑151 (MQ=255)                                                                                                                                               ggTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTAt       >  2:323022/1‑151 (MQ=255)                                                                                                                                                     ggCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGat  <  1:1225458/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                         TCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCC‑TTCTGGCAGATTAAG‑CCAAAATACCCGTTCGGCAATGCTT‑CAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTGAT  >  NC_000913/1245086‑1245383

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