BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,484,221:1	+T	100%	pseudogene (1/252 nt) pseudogene (750/750 nt)	gapC ← gapC ←	pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment) pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,484,221	1	.	T	99.4%	633.7 / ‑3.6	161	pseudogene (1/252 nt) pseudogene (750/750 nt)	gapC gapC	pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment) pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment)
Reads supporting (aligned to +/- strand): ref base . (0/1); new base T (94/66); total (94/67)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 173 total aligned reads displayed.

GCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑TTTTTT‑CCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTT  >  CP009273/1484073‑1484365

 cggTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tt                                                                                                                                                 >  1:349827/1‑150 (MQ=255)

   gTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttt                                                                                                                                               <  2:1098051/150‑1 (MQ=255)

   gTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttt                                                                                                                                               >  1:1388989/1‑150 (MQ=255)

    tAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑ttttt                                                                                                                                              <  1:107452/150‑1 (MQ=255)

     aaTTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttt                                                                                                                                             >  1:132347/1‑150 (MQ=255)

      aTTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttt                                                                                                                                             <  2:210775/149‑1 (MQ=255)

       tttCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttttc                                                                                                                                           >  2:1233205/1‑148 (MQ=255)

         tCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttttcc                                                                                                                                          <  1:326335/149‑4 (MQ=255)

         tCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttttcc                                                                                                                                          <  1:871726/149‑4 (MQ=255)

             tttGCGTGGCATCAAACACCGAACCGAAATCGCTGCCAATGATCTCGGAAGAGACTATTTCTTCATCGGTATAACCAAATCACTCGTTATTGGCGGTTGCTTGTTTAAGTGCGTTATTTACCTCTTCGGCAGTCACTTTTTTT‑CCGAGaa                                                                                                                                     >  2:1205182/1‑150 (MQ=255)

                 cGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGa                                                                                                                                 <  1:649057/150‑1 (MQ=255)

                  gTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGTGGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGa                                                                                                                                 >  2:786983/1‑149 (MQ=255)

                  gTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGac                                                                                                                                <  1:923626/150‑2 (MQ=255)

                     gCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGacac                                                                                                                              <  2:208812/149‑1 (MQ=255)

                     gCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACAcc                                                                                                                             <  1:239386/150‑1 (MQ=255)

                       aTCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAg                                                                                                                           >  2:421821/1‑150 (MQ=255)

                            acacCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAg                                                                                                                      <  2:481486/150‑1 (MQ=255)

                            acacCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAg                                                                                                                      <  2:31994/150‑1 (MQ=255)

                            acacCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑TTTTTTGCCGAGAATGGACACCAGCTCAg                                                                                                                      <  1:419897/150‑1 (MQ=255)

                               ccgaaccgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGa                                                                                                                   <  1:33823/150‑1 (MQ=255)

                               ccgaaccgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGa                                                                                                                   >  2:1165638/1‑150 (MQ=255)

                                 gaaccgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGAcc                                                                                                                 <  1:816571/150‑1 (MQ=255)

                                     cgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAAc                                                                                                             >  2:777057/1‑150 (MQ=255)

                                     cgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAAc                                                                                                             <  2:608188/150‑1 (MQ=255)

                                     cgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGCGGTTGCTTGGTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAAc                                                                                                             >  1:361672/1‑150 (MQ=255)

                                      gaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAAcc                                                                                                            <  2:1137493/150‑1 (MQ=255)

                                       aaaTGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCt                                                                                                           >  1:325347/1‑150 (MQ=255)

                                         aTGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACGCCAGCTCAGTGACCGAACCTGt                                                                                                         <  1:1157788/150‑1 (MQ=255)

                                           ggCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGttt                                                                                                       >  2:953280/1‑150 (MQ=255)

                                            gCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGtttt                                                                                                      >  1:749474/1‑150 (MQ=255)

                                             cTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTc                                                                                                     <  2:902308/150‑1 (MQ=255)

                                              tGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCa                                                                                                    >  1:432102/1‑150 (MQ=255)

                                              tGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCa                                                                                                    >  2:1450415/1‑150 (MQ=255)

                                                ccAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCa                                                                                                    <  1:104128/148‑1 (MQ=255)

                                                 cAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCg                                                                                                 >  1:65838/1‑150 (MQ=255)

                                                  aaTGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCgg                                                                                                >  2:1142429/1‑150 (MQ=255)

                                                  aaTGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCgg                                                                                                <  1:87873/150‑1 (MQ=255)

                                                   aTGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGc                                                                                               <  2:1239271/150‑1 (MQ=255)

                                                    tGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCa                                                                                              >  2:154126/1‑150 (MQ=255)

                                                        atCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCAcgcg                                                                                          >  2:278102/1‑150 (MQ=255)

                                                          cGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGtt                                                                                        >  1:1110140/1‑150 (MQ=255)

                                                           ggAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTg                                                                                       >  2:661653/1‑150 (MQ=255)

                                                            gAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCCAGCGTTgc                                                                                      >  2:1347814/1‑150 (MQ=255)

                                                             aaGAGACTATTTCCTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTgcg                                                                                     >  2:687489/1‑150 (MQ=255)

                                                              agagACTATTTCTTCTTCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTgcg                                                                                     >  1:97552/1‑149 (MQ=255)

                                                               gagaCTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCa                                                                                   >  1:520938/1‑150 (MQ=255)

                                                                agaCTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCAt                                                                                  <  2:702435/150‑1 (MQ=255)

                                                                agaCTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCAt                                                                                  <  1:201217/150‑1 (MQ=255)

                                                                  aCTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGa                                                                                >  1:1415178/1‑150 (MQ=255)

                                                                   cTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGAc                                                                               <  1:996039/150‑1 (MQ=255)

                                                                   cTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGAc                                                                               >  2:307791/1‑150 (MQ=255)

                                                                      tttCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CTGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCtt                                                                            >  2:596130/1‑150 (MQ=255)

                                                                      tttCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCtt                                                                            >  2:1338053/1‑150 (MQ=255)

                                                                      tttCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCtt                                                                            >  1:21690/1‑150 (MQ=255)

                                                                         cttcATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCttt                                                                           >  2:978594/1‑148 (MQ=255)

                                                                           tcATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTg                                                                                       >  2:837685/1‑134 (MQ=255)

                                                                            cATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTgc                                                                                      <  1:837685/134‑1 (MQ=255)

                                                                             aTCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGtt                                                                      <  1:1197681/149‑1 (MQ=255)

                                                                                   ataACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCt                                                               >  1:1228031/1‑150 (MQ=255)

                                                                                   ataACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCt                                                               >  1:590087/1‑150 (MQ=255)

                                                                                      aCCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAg                                                            >  1:892376/1‑150 (MQ=255)

                                                                                      aCCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAg                                                            >  2:128209/1‑150 (MQ=255)

                                                                                       ccAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGt                                                           <  2:132347/150‑1 (MQ=255)

                                                                                       ccAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGt                                                           >  2:1436526/1‑150 (MQ=255)

                                                                                        cAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGtt                                                          <  2:852095/150‑1 (MQ=255)

                                                                                        cAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGtt                                                          >  1:1457378/1‑150 (MQ=255)

                                                                                         aaaTGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTc                                                         <  2:1388989/150‑1 (MQ=255)

                                                                                          aaTGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTcc                                                        <  2:1077375/150‑1 (MQ=255)

                                                                                          aaTGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTcc                                                        <  2:1051160/150‑1 (MQ=255)

                                                                                              aCTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCggg                                                     <  2:970397/149‑1 (MQ=255)

                                                                                                  gTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCAc                                                >  2:1313731/1‑150 (MQ=255)

                                                                                                   ttATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCAcc                                               <  1:481493/150‑1 (MQ=255)

                                                                                                    tATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCAc                                                >  1:454223/1‑148 (MQ=255)

                                                                                                     aTTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCAc                                                >  1:1088599/1‑147 (MQ=255)

                                                                                                      ttGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGa                                            <  2:850225/150‑1 (MQ=255)

                                                                                                        ggtggtTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGAc                                           <  2:909357/149‑1 (MQ=255)

                                                                                                          tggtTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCaa                                        >  1:1409469/1‑150 (MQ=255)

                                                                                                            gtTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAAt                                       <  1:260575/149‑1 (MQ=255)

                                                                                                               gCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCt                                   >  1:1108905/1‑150 (MQ=255)

                                                                                                                 ttGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCttt                                 >  2:950563/1‑150 (MQ=255)

                                                                                                                  tGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCtttt                                >  1:1270867/1‑150 (MQ=255)

                                                                                                                   gTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCtttt                                <  1:1461121/149‑1 (MQ=255)

                                                                                                                     ttAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTgcc                             >  1:1391541/1‑150 (MQ=255)

                                                                                                                        aGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCAAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTgccgcc                          >  1:1251483/1‑150 (MQ=255)

                                                                                                                         gTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGccc                         >  2:259628/1‑150 (MQ=255)

                                                                                                                         gTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGccc                         >  2:1072360/1‑150 (MQ=255)

                                                                                                                             gTTATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCTGGATCACCAGACCAATGGCTTTTGCCGCCCCCGt                     <  2:1114970/150‑1 (MQ=255)

                                                                                                                              ttATTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCg                      >  1:727293/1‑148 (MQ=255)

                                                                                                                                aTTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAgt                  >  2:621767/1‑150 (MQ=255)

                                                                                                                                aTTCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAg                   <  2:581179/149‑1 (MQ=255)

                                                                                                                                 ttCACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAgtg                 <  1:1423773/150‑1 (MQ=255)

                                                                                                                                   cACCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAgtgt                <  2:358463/149‑1 (MQ=255)

                                                                                                                                    aCCTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGtgtg               <  1:687489/149‑1 (MQ=255)

                                                                                                                                     ccTCTTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTggg             >  2:805755/1‑150 (MQ=255)

                                                                                                                                       tctTCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTgggg            >  2:529538/1‑149 (MQ=255)

                                                                                                                                         ttCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGa           <  2:520938/148‑1 (MQ=255)

                                                                                                                                          tCGGCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAAt         >  1:407541/1‑149 (MQ=255)

                                                                                                                                            ggCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGa       <  1:1339884/149‑1 (MQ=255)

                                                                                                                                             gCAGTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGata     <  1:1240310/150‑1 (MQ=255)

                                                                                                                                                gTCACTTTTTTT‑CCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATAttt  <  1:278102/150‑1 (MQ=255)

GCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑TTTTTT‑CCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTT  >  CP009273/1484073‑1484365

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