BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,821,489	T→A	5.9%	I21N (ATT→AAT)	gmk →	guanylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,821,489	0	T	A	5.9%	535.8 / 14.7	167	I21N (ATT→AAT)	gmk	guanylate kinase
Reads supporting (aligned to +/- strand): ref base T (79/78); new base A (4/6); total (83/84)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 168 total aligned reads displayed.

TGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATG  >  NC_000913/3821341‑3821635

      aGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCt                                                                                                                                             >  1:1149600/1‑150 (MQ=255)

        ggCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCttt                                                                                                                                           >  1:865537/1‑150 (MQ=255)

        ggCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCttt                                                                                                                                           >  1:27608/1‑150 (MQ=255)

        ggCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATCTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCGCAAGCCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtgg                                                                                                                                           <  2:460796/150‑3 (MQ=255)

         tcTACTTCACATGTCCGCGCCTTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTtta                                                                                                                                          <  1:381862/149‑1 (MQ=255)

          cTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtggat                                                                                                                                         <  2:90610/150‑5 (MQ=255)

           tACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGAATCAGGCtggatt                                                                                                                                        <  2:297185/150‑6 (MQ=255)

           tACGTCACATGGCCGCGCCCTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTttatt                                                                                                                                        >  2:1184274/1‑150 (MQ=255)

            aCGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTttatta                                                                                                                                       <  2:841596/150‑1 (MQ=255)

             cGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaa                                                                                                                                      <  2:1387915/150‑1 (MQ=255)

             cGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCAAATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaa                                                                                                                                      >  2:1388143/1‑150 (MQ=255)

              gTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaa                                                                                                                                     >  1:559350/1‑150 (MQ=255)

               tCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaa                                                                                                                                    >  1:612434/1‑150 (MQ=255)

               tCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaa                                                                                                                                    <  1:98419/150‑1 (MQ=255)

               tCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaa                                                                                                                                    >  1:438022/1‑150 (MQ=255)

               tCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTaaaa                                                                                                                                    >  1:1086354/1‑150 (MQ=255)

                 acaTGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAAc                                                                                                                                  <  2:374453/150‑1 (MQ=255)

                       ccGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAAcc                                                                                                                            >  1:279570/1‑150 (MQ=255)

                         gcgcCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGt                                                                                                                          >  2:1297512/1‑150 (MQ=255)

                           gcCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTg                                                                                                                        <  2:806270/150‑1 (MQ=255)

                           gcCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTg                                                                                                                        <  1:1266530/150‑1 (MQ=255)

                               tgtATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTAGTATg                                                                                                                    <  2:559574/150‑1 (MQ=255)

                                 tataATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCACCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGac                                                                                                                  >  1:940512/1‑150 (MQ=255)

                                  ataataAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGaca                                                                                                                 >  1:270761/1‑150 (MQ=255)

                                           tCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGttt                                                                                                        <  2:629093/150‑1 (MQ=255)

                                                  tAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTc                                                                                                 <  1:796376/150‑1 (MQ=255)

                                                   aGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTca                                                                                                <  2:1270135/150‑1 (MQ=255)

                                                    ggCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcac                                                                                               >  1:435064/1‑150 (MQ=255)

                                                    ggCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcac                                                                                               <  2:710429/150‑1 (MQ=255)

                                                    ggCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcac                                                                                               <  1:696839/150‑1 (MQ=255)

                                                    ggCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcac                                                                                               <  1:213315/150‑1 (MQ=255)

                                                     gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGt                                                                                                    <  1:1051714/144‑1 (MQ=255)

                                                     gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcaca                                                                                              >  2:113370/1‑150 (MQ=255)

                                                     gCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTcaca                                                                                              >  1:85376/1‑150 (MQ=255)

                                                             tCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCAcacaccacac                                                                                      <  1:1184274/150‑1 (MQ=255)

                                                             tCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCAcacaccacac                                                                                      >  1:717094/1‑150 (MQ=255)

                                                               gCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGc                                                                                    <  1:777703/150‑1 (MQ=255)

                                                                  aaTCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCaa                                                                                 <  1:113370/150‑1 (MQ=255)

                                                                   aTCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAAc                                                                                >  1:197063/1‑150 (MQ=255)

                                                                     cacaTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcg                                                                              <  2:197063/150‑1 (MQ=255)

                                                                       caTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACcgcg                                                                            <  2:27608/150‑1 (MQ=255)

                                                                        aTACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGt                                                                           >  1:621202/1‑150 (MQ=255)

                                                                         tACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTc                                                                          >  2:1391200/1‑150 (MQ=255)

                                                                           cGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCt                                                                        >  1:1004239/1‑150 (MQ=255)

                                                                             aaaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTgg                                                                      <  1:413139/150‑1 (MQ=255)

                                                                             aaaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTgg                                                                      <  2:1199180/150‑1 (MQ=255)

                                                                              aaGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCTGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGt                                                                     >  2:737529/1‑150 (MQ=255)

                                                                                gATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGt                                                                     >  2:379692/1‑148 (MQ=255)

                                                                                gATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGt                                                                     <  1:379692/148‑1 (MQ=255)

                                                                                  tACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAg                                                                 <  2:612434/150‑1 (MQ=255)

                                                                                  tACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAg                                                                 <  2:1028866/150‑1 (MQ=255)

                                                                                    cTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTc                                                               >  2:1278201/1‑150 (MQ=255)

                                                                                      cATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCa                                                             >  1:590119/1‑150 (MQ=255)

                                                                                      cATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCa                                                             >  2:802749/1‑150 (MQ=255)

                                                                                      cATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACACAACCGTTGTATGAAACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCATGGTGAAGTCCa                                                             >  1:747670/1‑150 (MQ=255)

                                                                                           cTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTg                                                        <  1:1270924/150‑1 (MQ=255)

                                                                                           cTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTg                                                        <  1:1388143/150‑1 (MQ=255)

                                                                                            tCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGa                                                       <  2:314308/150‑1 (MQ=255)

                                                                                              aaGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAAc                                                     <  2:1149600/150‑1 (MQ=255)

                                                                                                ggCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACat                                                   >  1:923586/1‑150 (MQ=255)

                                                                                                 gCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACatt                                                  <  2:1207800/150‑1 (MQ=255)

                                                                                                  cACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACatta                                                 >  2:302332/1‑150 (MQ=255)

                                                                                                   aCGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACattat                                                <  2:582297/150‑1 (MQ=255)

                                                                                                      cTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATAACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTc                                             >  1:317788/1‑150 (MQ=255)

                                                                                                      cTTTATATTGTTTCTGCCCCCAGTGGAGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTc                                             >  2:862731/1‑150 (MQ=255)

                                                                                                        ttATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCtt                                           >  2:206363/1‑150 (MQ=255)

                                                                                                        ttATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCtt                                           <  2:187299/150‑1 (MQ=255)

                                                                                                          atatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTg                                         <  1:468683/150‑1 (MQ=255)

                                                                                                          atatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTg                                         <  2:783793/150‑1 (MQ=255)

                                                                                                          atatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTg                                         <  2:433144/150‑1 (MQ=255)

                                                                                                           tatTGTTTCTGCCCCCAGTGGCGCGGGTTAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGt                                        >  1:217512/1‑150 (MQ=255)

                                                                                                           tatTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGt                                        >  1:91475/1‑150 (MQ=255)

                                                                                                            atTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGtt                                       >  1:1009657/1‑150 (MQ=255)

                                                                                                              tGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTaa                                     >  2:1123704/1‑150 (MQ=255)

                                                                                                                tttCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATc                                   <  2:1086354/150‑1 (MQ=255)

                                                                                                                   cTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACAGTTGTATTACACCCAGGTTTCTGTTTAACACACCACACGCCAACCGCGTCCTGGTGAAGTACAAGGTGAAAATTATTTCTTTGTTAATCata                                >  2:1104856/1‑149 (MQ=255)

                                                                                                                     gCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTGAATCatgat                              >  1:246567/1‑150 (MQ=255)

                                                                                                                     gCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGATTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgat                              >  2:1303788/1‑150 (MQ=255)

                                                                                                                      cccccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCatgatg                             >  2:1098307/1‑150 (MQ=255)

                                                                                                                       ccccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACgg                                                          >  1:859146/1‑120 (MQ=255)

                                                                                                                       ccccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACgg                                                          <  2:859146/120‑1 (MQ=255)

                                                                                                                        cccAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTCTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGaa                           <  2:1351905/150‑1 (MQ=255)

                                                                                                                          cAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAAtt                         <  1:1303708/150‑1 (MQ=255)

                                                                                                                           aGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACGCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAAttt                        <  2:126993/150‑1 (MQ=255)

                                                                                                                            gTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTa                       <  2:496741/150‑1 (MQ=255)

                                                                                                                              ggCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTaaa                     <  2:590119/150‑1 (MQ=255)

                                                                                                                                cgcgGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTca                                                                                                <  1:1204915/73‑1 (MQ=255)

                                                                                                                                   gggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAAt                <  1:1391200/150‑1 (MQ=255)

                                                                                                                                    ggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATg               >  2:622122/1‑150 (MQ=255)

                                                                                                                                    ggTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATg               <  2:85376/150‑1 (MQ=255)

                                                                                                                                       aaaTCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTg                                         <  1:571443/121‑1 (MQ=255)

                                                                                                                                       aaaTCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGAtt            >  1:1427239/1‑150 (MQ=255)

                                                                                                                                        aataaaGCCTGAATCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTa           >  1:532106/6‑150 (MQ=255)

                                                                                                                                        aataaaGCCTGAATCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTa           >  1:529988/6‑150 (MQ=255)

                                                                                                                                        aataaaGCCTGAATCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTa           >  1:545071/6‑150 (MQ=255)

                                                                                                                                         aTCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATCAg          >  1:1316887/1‑150 (MQ=255)

                                                                                                                                           ccAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCa        <  1:1278201/150‑1 (MQ=255)

                                                                                                                                            cAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCag       >  1:1088764/1‑150 (MQ=255)

                                                                                                                                                 tGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATg  <  1:784362/150‑1 (MQ=255)

                                                                                                                                                 tGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATg  <  2:940512/150‑1 (MQ=255)

TGAGCAAGGGCTACGTCACATGGCCGCGCCGTGTATAATAAGCTCGTATGTAGGCTTTATTTCGCTAATCACATACGAAAGATACTCATGGCTCAAGGCACGCTTTATATTGTTTCTGCCCCCAGTGGCGCGGGTAAATCCAGCCTGATTCAGGCTTTATTAAAAACCCAACCGTTGTATGACACCCAGGTTTCTGTTTCACACACCACACGCCAACCGCGTCCTGGTGAAGTCCACGGTGAACATTATTTCTTTGTTAATCATGATGAATTTAAAGAAATGATTAGCAGAGATG  >  NC_000913/3821341‑3821635

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