BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	40.6%	85.0 / 34.7	71	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (34/7); new base A (0/28); total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.81e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCACC  >  NC_000913/1222447‑1222745

gCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑Cg                                                                                                                                                         >  2:62368/1‑150 (MQ=255)

    cGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAaaac                                                                                                                                                    <  2:1446887/151‑1 (MQ=255)

              acccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑cc                                                                                                                                         <  1:668204/150‑1 (MQ=255)

                 cattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGcc                                                                                                                                     >  2:1131890/1‑151 (MQ=255)

                  attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGcc                                                                                                                                     >  1:526874/1‑150 (MQ=255)

                    taccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGt                                                                                                                                  >  1:830449/1‑151 (MQ=255)

                      ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTa                                                                                                                                 >  2:935486/1‑150 (MQ=255)

                       cGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTaa                                                                                                                                >  1:1299159/1‑150 (MQ=255)

                       cGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTaa                                                                                                                                >  1:844326/1‑150 (MQ=255)

                         caacaaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGc                                                                                                                             >  2:708680/1‑151 (MQ=255)

                          aacaaTTTCTGCACCTTGTGCGCTGTCGTTTTTTGCTTTTTCGTTGTCTTCTGCGCAATAAATATGTG‑TTTCTGACGGAAACCT‑T‑TTTTT‑AA‑TTTTTTACTCCTGCGCCCAGG‑TCTTT‑GG‑T‑G‑CGAAAAAAAACAG‑G‑CCCGCCAGTAAAg                                                                                                                              <  1:926058/149‑1 (MQ=255)

                          aacaaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGc                                                                                                                             >  2:1390963/1‑150 (MQ=255)

                          aacaaTTTCTGAAACCTGTGCGCTGTCATTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGAGGGCAACCC‑T‑TTTTT‑AA‑TATTGTAGTCCGGCGCCCAGG‑TCATT‑GG‑T‑G‑AAAAAAAAAACAG‑G‑CCCGCCAGTAAAGc                                                                                                                             <  1:724045/150‑1 (MQ=255)

                                ttCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAAtttt                                                                                                                      >  2:1339807/1‑151 (MQ=255)

                                     cAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTtcat                                                                                                                 >  2:733257/1‑151 (MQ=255)

                                      acctggtgCGTTGTAATTTTTGCCCTTTTCTTTGGTTTCTGCCCAAGAAATAGAGT‑TTCCTGGCGGAAACCC‑T‑TTTTT‑AA‑TATTGTATCCCTGCCCCCAGG‑TCTTT‑GG‑G‑G‑AGAAAACAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTtca                                                                                                                  <  1:187460/144‑1 (MQ=255)

                                      acctggtgCGCTGTTATTTTTTGCTTTTTCTTTGGTTTCTGCAAAAGAAATATAGT‑TTTCTGACGGAAACCC‑G‑TTTTT‑AA‑TTTTGTACCCCGGCCCCCATT‑TTTTG‑GT‑T‑G‑CAAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTtca                                                                                                                  <  1:1361561/144‑1 (MQ=255)

                                      aaCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                                >  1:1049961/1‑151 (MQ=255)

                                          tgtgCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTtcatcatc                                                                                                             >  1:946283/1‑150 (MQ=255)

                                           gtgcgttttaattttttCCCATTTCTTTGGCTTCTGCGAAAGTAATATAGG‑TTTCTGCGGACAACCC‑T‑TTTTA‑AA‑TATTGTAGCCCTCCGCCCAGT‑TTATG‑GG‑T‑G‑CAAAAACAACCAG‑G‑CCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                           <  1:332598/139‑1 (MQ=255)

                                           gggcgttttCAGTTTTTGCCATTTTTTTGGTTTCTGCGAAAGAAATATAGT‑TTCCTGACGGCAACCC‑T‑TTTTA‑AA‑TTTTGTACCCCTCCCCCCAGG‑TATTT‑GG‑T‑G‑CAAAAACAACCAG‑G‑CCCGCCAGTAAAGCAAATTTTTtcatcat                                                                                                              <  2:1208142/140‑1 (MQ=255)

                                              cgttgtctttttttGCTTTTTCTTGGGCTTCTGCAAAAAAAATATGGT‑TTTCGGGCGGCAACCC‑T‑TTTTA‑AA‑TTTTGTAGCCCT‑CCCCCATT‑TTTTT‑GG‑G‑G‑GCAAAAAAAACCG‑GCCCCGCCAGTAAAGCAAATTTTTtcatcatc                                                                                                             <  1:587754/137‑1 (MQ=255)

                                               gttttCATTTTTTGCTTTTTTGTTGGTTTTTGCGAAATTAATAT‑GG‑TTTATTGCCGACAACC‑T‑TTTTT‑AATTTTTGTGTTCCTGCGCCAAGT‑TCTTG‑GT‑T‑G‑CAAAAACAACCAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                       <  2:280761/147‑1 (MQ=255)

                                                 tttCAGTTTTTCCTATTTCTTTGTCTTCTCCGAAAAAAATATAGG‑TTCCTGAAAACAACCT‑T‑TTTTT‑AA‑TATTATCCTCCCCCCCCCA‑G‑TTTTT‑GG‑TCA‑CAAAAACAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                       <  1:473378/147‑1 (MQ=255)

                                                   tCATTTTTTGTCATTGCGTTGTCTTCTGCGAAATTAATATAGT‑TTTTTGAGGACAACCC‑G‑TTTTT‑AA‑TAGTGGATTCCGGAGCCTAGT‑TTTTT‑GT‑G‑G‑CCAAAAAAATATG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTcc                                                                                                        <  1:1007562/146‑1 (MQ=255)

                                                   tCAGTTTTTGCCTTTTCGTTGTTTTCTGCGCAAAAAATATGGG‑TTTCGGCCGG‑AAACC‑T‑TTTTTAAA‑TTTTGTAGCCCGGCGCCCAGG‑TCTTG‑GT‑T‑G‑CAAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                                    <  1:1244931/150‑1 (MQ=255)

                                                      ttttttGCTTTTTGTTTTGCTTCTGCGAAAGTAATATGGG‑TTTCTGACGCCAACCC‑T‑TTTTT‑AA‑TTTTGTAGCCCTGC‑CCCAAT‑TTCTT‑GG‑T‑GCCAAAAACAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                                      <  1:902734/144‑1 (MQ=255)

                                                      gTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                                 >  2:542872/1‑150 (MQ=255)

                                                       tttttGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAACCC‑T‑TTTTT‑AA‑TATTGTAGCCCGGCCCCCAGG‑TCATT‑GG‑T‑G‑CAAAAACAAACGG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                                 <  2:839121/149‑1 (MQ=255)

                                                       tttttGCCATTTCGTTGGCTTCGGCGCAAGAAATATAGT‑TTTCTGACGGAAAACC‑T‑TTTTT‑AA‑TATTGTAGCCCGGCGCCCAGT‑TCTTT‑GG‑T‑G‑CAAAAACAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCa                                                                                                <  2:1205562/150‑1 (MQ=255)

                                                          ttGCCTTTTCGTTGGTTTCGGCGCAAGTAATATAGG‑TTTCTGACGGAAAACC‑T‑TTTTA‑AA‑TATTGTACTCCTGCACAAAGG‑TAATT‑GG‑T‑G‑AAAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                              <  1:937015/149‑1 (MQ=255)

                                                           tGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtttt                                                                                            >  1:912655/1‑150 (MQ=255)

                                                            gCCATTTCGTTGGTTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAACCC‑T‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CAAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtttt                                                                                            <  2:211367/149‑1 (MQ=255)

                                                             ccttttCGTTGGCTTCGCCACAATTAATATAGT‑TTTCTGCCGACAACCC‑T‑TTTTT‑AA‑TATTGTAGTCCTGCCCCCAGG‑TCTTT‑GG‑T‑G‑CAAAAACAACCAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAc                                                                                          <  2:1349452/147‑1 (MQ=255)

                                                             ccATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                         >  2:396550/1‑151 (MQ=255)

                                                             ccATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                         >  1:650285/1‑151 (MQ=255)

                                                               aTTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                        >  1:1284024/1‑150 (MQ=255)

                                                               aTTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCACCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                        >  1:1330757/1‑150 (MQ=255)

                                                                  tCGTTGTTTTCTGCGAAAGTAATATAGT‑TTTCTGACGGCAAACC‑T‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGT‑TCATT‑GT‑T‑G‑CAAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGa                                                                                     <  1:784108/150‑1 (MQ=255)

                                                                  tCGTTGGGTTCTGCGAAGAAAATATAGG‑TTTCTGACGGCAACCC‑T‑TTTTT‑AA‑TATTTTAGTCCT‑CCCCCAGGTTTTTT‑GG‑T‑G‑AAAAAAAAAACGG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGa                                                                                     <  1:832831/150‑1 (MQ=255)

                                                                     tggacTTCTGCCAAAGAAAT‑TAGGATTTCTGGGGGCAATCC‑T‑TTTTA‑AA‑TTTTGTAATCCT‑CGCCCAAG‑TTTTT‑GGTT‑G‑CAAAACCAAAAAG‑A‑ACCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGaa                                                                                    <  1:357671/144‑1 (MQ=255)

                                                                      tGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                                >  2:631776/1‑151 (MQ=255)

                                                                       ggtttCTGCGCAAGTAATATAGG‑TTTCTGCCGGAAACCC‑T‑TTTTT‑AA‑TATTGTAGTCCTGCCCCAAGT‑TCTTG‑GG‑T‑G‑CAAAAACAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGca                                                                               <  1:185202/148‑1 (MQ=255)

                                                                          ttCGGCGGAAGTAATATGGT‑TTTCGATCAGCAACCC‑T‑TTTTA‑AA‑TATTGGAGTCCTGGGCCCAGG‑TAATG‑GT‑A‑G‑CGAAAAC‑AACAGTC‑CCCGCCAGTAAAGCACATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGca                                                                               <  1:210410/148‑1 (MQ=255)

                                                                           tttGCGCAATAATTATATG‑TTTCTGACGGAAAGCT‑T‑TTTTT‑AA‑TTTTGTCTCCCTGCCCCCGGT‑TTTTT‑GT‑T‑G‑AGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacca                                                                            <  2:287543/148‑1 (MQ=255)

                                                                             tGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCAccc                                                                         >  2:1193087/1‑151 (MQ=255)

                                                                               cgaaattAATATAGG‑TTTCTGACGGCAACCT‑T‑TTTTT‑AA‑TATTGTATTCCTGCCCCCAGG‑TCATT‑GG‑T‑G‑AAAAAAAAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCg                                                                        <  2:547726/144‑1 (MQ=255)

                                                                                  aattAATATAGT‑TTTC‑GAACGCAAACCTT‑TTTTT‑AA‑TTTTGTAGTCCTGCCCCCAGT‑TTTTT‑GG‑T‑G‑AAAAAAAAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAcc                                                                    <  1:856043/148‑1 (MQ=255)

                                                                                  aaGTAATATAGT‑TTTCTGACGGCAACCC‑T‑TTTTT‑AA‑TATTGTAGCCCAGCCCCCAGG‑TCATT‑GG‑T‑G‑CAAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTa                                                                      <  2:186754/149‑1 (MQ=255)

                                                                                            gg‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACg                                                          >  2:181546/1‑151 (MQ=255)

                                                                                                ttCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                                        >  1:612439/1‑150 (MQ=255)

                                                                                                ttCTGA‑GGCCAACC‑GTTTTTT‑AT‑TGTTTTTGCCCTCCCCCCAGT‑TCTTT‑GG‑T‑‑‑AGAAAAAAAACCG‑GCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCAc                                                           <  1:322250/147‑1 (MQ=255)

                                                                                                    gACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCtt                                                   >  1:1285548/1‑151 (MQ=255)

                                                                                                            aaacgtttttt‑AA‑TTTGTTAGCCCGGCGCCCAGG‑TTTGT‑GG‑T‑G‑AAAAAAAAAACGG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCAcc                                                     <  1:977888/137‑1 (MQ=255)

                                                                                                             cc‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGccc                                          >  1:1515391/1‑151 (MQ=255)

                                                                                                              c‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGcccc                                         >  1:1053900/1‑151 (MQ=255)

                                                                                                                  ttttt‑AA‑TTTTGTAGTCCTCCGCCAAGT‑TCTTG‑GG‑G‑G‑AAAAAAAAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCt                                        <  2:1024352/150‑1 (MQ=255)

                                                                                                                  tttta‑AA‑TTTTGT‑GTCCCGCCCCCAGG‑TTATTGGG‑G‑G‑AAAAAAAAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGcccc                                         <  1:708637/149‑1 (MQ=255)

                                                                                                                    ttt‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                     >  2:1262070/1‑151 (MQ=255)

                                                                                                                         aatttggTATTCCTGCCCCCGGG‑TC‑TT‑GG‑TGG‑AAAAAAAAAACGG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAgg                                  <  1:269650/145‑1 (MQ=255)

                                                                                                                         a‑tatTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGc                                 >  2:77947/1‑151 (MQ=255)

                                                                                                                             ttGTAGCCCTGCCCCCGGG‑TATTG‑GG‑T‑G‑CAAAAAAAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGc                                 <  2:830496/148‑1 (MQ=255)

                                                                                                                             ttGTAGCCCTGCCCCCAGG‑TCTTT‑GG‑T‑G‑CAAAAAAAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                                <  2:859717/149‑1 (MQ=255)

                                                                                                                              tGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAtt                             >  2:1039117/1‑151 (MQ=255)

                                                                                                                               gTAGCCCGGCCCCCAGT‑TCTTT‑GG‑T‑G‑AAAAAAAAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAtt                             <  1:802595/150‑1 (MQ=255)

                                                                                                                                  gTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCg                         >  1:107682/1‑151 (MQ=255)

                                                                                                                                    ccTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGt                       >  1:877690/1‑151 (MQ=255)

                                                                                                                                          cccAGG‑TCTTT‑GG‑T‑G‑CAAAAAAAAACAG‑C‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGaa                  <  1:625993/150‑1 (MQ=255)

                                                                                                                                          cccAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAAc                 >  1:195575/1‑151 (MQ=255)

                                                                                                                                                  cATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg          >  1:695952/1‑151 (MQ=255)

                                                                                                                                                    tt‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc        >  2:1251925/1‑151 (MQ=255)

                                                                                                                                                     t‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGAc        >  2:844471/1‑150 (MQ=255)

                                                                                                                                                               gAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAcc  >  2:243056/1‑150 (MQ=255)

GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGG‑TTTCTGACGGCAAACC‑G‑TTTTT‑AA‑TATTGTAGTCCTGCGCCCAGG‑TCATT‑GG‑T‑G‑CGAAAACAAACAG‑G‑CCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCACC  >  NC_000913/1222447‑1222745

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