BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,668,886	T→C	A101A (GCA→GCG)	ccp ←	cytochrome c peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,668,886	0	T	C	99.7%	1305.8 / ‑5.0	305	A101A (GCA→GCG)	ccp	cytochrome c peroxidase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (162/142); minor base A (0/1); total (162/143)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 308 total aligned reads displayed.

GTAGCGCGGTATAACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGTGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTA  >  NC_000913/3668789‑3668983

 tAGCGCGGTATAACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCcg                                                                                               >  2:5396138/1‑101 (MQ=255)

   gcgcGGTATAACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCcgca                                                                                             >  1:5396142/1‑101 (MQ=255)

    cgcgGTATAACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCcgcac                                                                                            >  2:5396143/1‑101 (MQ=255)

       ggTATAACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCcgcacgca                                                                                         >  1:5396149/1‑101 (MQ=255)

        gTATAACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCgcacgcac                                                                                        <  1:5395780/101‑1 (MQ=255)

         tataACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGcacgcacg                                                                                       <  2:5395673/101‑1 (MQ=255)

         tataACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGcacgcacg                                                                                       >  2:5396150/1‑101 (MQ=255)

          ataACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACgg                                                                                      <  2:5395706/101‑1 (MQ=255)

            aaCGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGcc                                                                                    <  1:5395553/101‑1 (MQ=255)

            aaCGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGcc                                                                                    >  2:5396156/1‑101 (MQ=255)

            aaCGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGcc                                                                                    >  1:5396155/1‑101 (MQ=255)

             aCGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCt                                                                                   <  2:5395841/101‑1 (MQ=255)

                cgTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCaa                                                                                >  2:5396161/1‑101 (MQ=255)

                 gTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAg                                                                               >  1:5396163/1‑101 (MQ=255)

                 gTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAg                                                                               <  1:5395818/101‑1 (MQ=255)

                  ttGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAgg                                                                              >  2:5396166/1‑101 (MQ=255)

                   tggtggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGt                                                                             <  1:5395654/101‑1 (MQ=255)

                    ggtggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGtt                                                                            <  1:5395924/101‑1 (MQ=255)

                    ggtggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGtt                                                                            >  1:5396171/1‑101 (MQ=255)

                    ggtggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGtt                                                                            <  2:5395427/101‑1 (MQ=255)

                    ggtggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGtt                                                                            <  1:5395767/101‑1 (MQ=255)

                     gtggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTa                                                                           <  1:5395723/101‑1 (MQ=255)

                     gtggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTa                                                                           >  2:5396173/1‑101 (MQ=255)

                      tggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTaa                                                                          >  2:5396177/1‑101 (MQ=255)

                       ggCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTaaa                                                                         <  1:5395836/101‑1 (MQ=255)

                        gCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTaaaa                                                                        <  1:5396048/101‑1 (MQ=255)

                        gCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTaaaa                                                                        >  1:5396179/1‑101 (MQ=255)

                         cATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAg                                                                       <  1:5395738/101‑1 (MQ=255)

                          aTAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTaaa                                                                         <  1:5395789/98‑1 (MQ=255)

                           tAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGAt                                                                     >  1:5396184/1‑101 (MQ=255)

                            aGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGAtt                                                                    >  2:5396185/1‑101 (MQ=255)

                             gTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGAttt                                                                   <  1:5395989/101‑1 (MQ=255)

                               ttCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTtat                                                                 >  2:5396190/1‑101 (MQ=255)

                                tCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTtata                                                                >  1:5396191/1‑101 (MQ=255)

                                   tACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAg                                                                               >  1:5396195/1‑83 (MQ=255)

                                     cTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCaa                                                           <  2:5395619/101‑1 (MQ=255)

                                     cTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCaa                                                           >  1:5396199/1‑101 (MQ=255)

                                     cTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCaa                                                           <  1:5395642/101‑1 (MQ=255)

                                      tGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAg                                                          <  2:5395827/101‑1 (MQ=255)

                                      tGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAg                                                          >  2:5396202/1‑101 (MQ=255)

                                       gCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGc                                                         >  1:5396203/1‑101 (MQ=255)

                                        cATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCt                                                        <  2:5395951/101‑1 (MQ=255)

                                         aTCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCtt                                                       <  1:5395797/101‑1 (MQ=255)

                                         aTCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCtt                                                       >  1:5396212/1‑101 (MQ=255)

                                         aTCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCtt                                                       >  2:5396207/1‑101 (MQ=255)

                                         aTCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCtt                                                       >  2:5396211/1‑101 (MQ=255)

                                            aCCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGc                                                         <  2:5395554/96‑1 (MQ=255)

                                              ccATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGt                                                  >  1:5396217/1‑101 (MQ=255)

                                              ccATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGt                                                  >  2:5396219/1‑101 (MQ=255)

                                              ccATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGt                                                  <  1:5395889/101‑1 (MQ=255)

                                               cATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTc                                                 <  1:5395755/101‑1 (MQ=255)

                                               cATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTc                                                 <  1:5395731/101‑1 (MQ=255)

                                                aTTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCg                                                <  1:5395676/101‑1 (MQ=255)

                                                aTTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCg                                                >  2:5396223/1‑101 (MQ=255)

                                                aTTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCg                                                <  2:5395650/101‑1 (MQ=255)

                                                 ttCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGt                                               >  2:5396226/1‑101 (MQ=255)

                                                   cAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTaa                                             >  2:5396229/1‑101 (MQ=255)

                                                    aaTCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAAt                                            >  1:5396231/1‑101 (MQ=255)

                                                     aTCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATc                                           <  2:5395905/101‑1 (MQ=255)

                                                     aTCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATc                                           >  2:5396234/1‑101 (MQ=255)

                                                       cTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGc                                                         >  2:5396240/1‑85 (MQ=255)

                                                       cTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCa                                         <  1:5395708/101‑1 (MQ=255)

                                                       cTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCa                                         <  1:5395871/101‑1 (MQ=255)

                                                        ttATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCAt                                        <  1:5395688/101‑1 (MQ=255)

                                                        ttATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCAt                                        <  1:5395744/101‑1 (MQ=255)

                                                        ttATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCAt                                        >  1:5396241/1‑101 (MQ=255)

                                                         tATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATc                                       >  2:5396246/1‑101 (MQ=255)

                                                          aTTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCa                                      <  2:5395754/101‑1 (MQ=255)

                                                          aTTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCa                                      <  2:5395810/101‑1 (MQ=255)

                                                          aTTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCa                                      >  1:5396249/1‑101 (MQ=255)

                                                           ttCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCaa                                     >  2:5396250/1‑101 (MQ=255)

                                                            tCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAAc                                    <  2:5396000/101‑1 (MQ=255)

                                                            tCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAAc                                    >  1:5396253/1‑101 (MQ=255)

                                                              aaaTCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTg                                  >  1:5396257/1‑101 (MQ=255)

                                                               aaTCGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGt                                 >  2:5396259/1‑101 (MQ=255)

                                                                  cGCTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTc                              >  2:5396262/1‑101 (MQ=255)

                                                                    cTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTcgc                            >  1:5396268/1‑101 (MQ=255)

                                                                    cTTTGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTcgc                            >  2:5396267/1‑101 (MQ=255)

                                                                     tttGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTcgcg                           <  2:5395974/101‑1 (MQ=255)

                                                                       tGCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGcc                         >  2:5396273/1‑101 (MQ=255)

                                                                        gCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCt                                                        <  2:5395862/69‑1 (MQ=255)

                                                                        gCGAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCa                        >  1:5396274/1‑101 (MQ=255)

                                                                          gAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAgg                      >  2:5396277/1‑101 (MQ=255)

                                                                          gAAACGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAgg                      <  2:5395777/101‑1 (MQ=255)

                                                                            aaCGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGaa                    >  1:5396282/1‑101 (MQ=255)

                                                                             aCGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAat                   >  2:5396284/1‑101 (MQ=255)

                                                                             aCGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAat                   <  2:5395701/101‑1 (MQ=255)

                                                                              cGGGTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAata                  >  1:5396286/1‑101 (MQ=255)

                                                                                ggTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAAtata                <  1:5395969/101‑1 (MQ=255)

                                                                                ggTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTACATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAAtata                <  1:5395511/101‑1 (MQ=255)

                                                                                 gTTTATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATataa               <  2:5395729/101‑1 (MQ=255)

                                                                                  tttATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTc                              <  2:5395719/85‑1 (MQ=255)

                                                                                    tATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATataataa            >  2:5396290/1‑101 (MQ=255)

                                                                                    tATCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATataataa            <  2:5395875/101‑1 (MQ=255)

                                                                                     aTCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAg                       >  1:5396294/1‑89 (MQ=255)

                                                                                     aTCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATataataat           <  1:5395806/101‑1 (MQ=255)

                                                                                      tCAGCCAGCAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATataataata          >  2:5396296/1‑101 (MQ=255)

                                                                                       cagccagcAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAg         <  2:5395891/101‑1 (MQ=255)

                                                                                        agccagcAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAgg        <  2:5395816/101‑1 (MQ=255)

                                                                                        agccagcAGCGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAgg        >  2:5396300/1‑101 (MQ=255)

GTAGCGCGGTATAACGCGTTGGTGGCATAGTTTCATACTGCATCACCCATTCAATCTTATTCAAATCGCTTTGCGAAACGGGTTTATCAGCCAGCAGTGCCGCACGCACGGCCTCAAGGTTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTA  >  NC_000913/3668789‑3668983

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