BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	447,207	C→T	A167T (GCT→ACT)	cyoE ←	heme O synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	447,207	0	C	T	99.6%	1003.5 / NA	238	A167T (GCT→ACT)	cyoE	heme O synthase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (123/114); minor base G (1/0); total (124/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

TACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCG  >  NC_000913/447109‑447304

 aCCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTcag                                                                                                >  1:662003/1‑101 (MQ=255)

  ccACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTcagc                                                                                               <  2:661615/101‑1 (MQ=255)

   cACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTcagca                                                                                              <  2:661440/101‑1 (MQ=255)

    aCTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTcagcag                                                                                             <  1:661724/101‑1 (MQ=255)

    aCTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTcagcag                                                                                             <  2:661755/101‑1 (MQ=255)

      tGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGa                                                                                           >  2:662008/1‑101 (MQ=255)

      tGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGa                                                                                           <  2:661443/101‑1 (MQ=255)

       ggCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGAt                                                                                          <  2:661500/101‑1 (MQ=255)

        gCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATc                                                                                         >  2:662010/1‑101 (MQ=255)

        gCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATc                                                                                         <  1:661631/101‑1 (MQ=255)

        gCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATc                                                                                         >  2:662012/1‑101 (MQ=255)

         cAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCg                                                                                        <  1:661552/101‑1 (MQ=255)

         cAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCg                                                                                        >  1:662017/1‑101 (MQ=255)

         cAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCg                                                                                        >  1:662014/1‑101 (MQ=255)

         cAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCg                                                                                        <  2:661651/101‑1 (MQ=255)

           aTACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCt                                                                                      >  2:662019/1‑101 (MQ=255)

           aTACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCt                                                                                      <  1:661536/101‑1 (MQ=255)

           aTACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCt                                                                                      <  2:661384/101‑1 (MQ=255)

            tACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTg                                                                                     <  2:661805/101‑1 (MQ=255)

            tACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTg                                                                                     <  2:661696/101‑1 (MQ=255)

             aCCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTgc                                                                                    >  1:662022/1‑101 (MQ=255)

              ccGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTgcg                                                                                   <  2:661597/101‑1 (MQ=255)

              ccGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTgcg                                                                                   <  2:661548/101‑1 (MQ=255)

               cGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTgcgc                                                                                  >  2:662024/1‑101 (MQ=255)

                 gAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGcgccg                                                                                >  1:662027/1‑101 (MQ=255)

                  aaTGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGcgccgc                                                                               >  2:662029/1‑101 (MQ=255)

                    tGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTa                                                                             <  2:661563/101‑1 (MQ=255)

                    tGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTa                                                                             >  1:662031/1‑101 (MQ=255)

                      tttGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATc                                                                           <  1:661571/101‑1 (MQ=255)

                      tttGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATc                                                                           >  1:662033/1‑101 (MQ=255)

                      tttGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATc                                                                           <  1:661649/101‑1 (MQ=255)

                         gccgccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGaa                                                                        >  2:662037/1‑101 (MQ=255)

                          ccgccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAAc                                                                       >  1:662039/1‑101 (MQ=255)

                           cgccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAAcc                                                                      >  2:662041/1‑100 (MQ=255)

                            gccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTc                                                                     >  1:662042/1‑101 (MQ=255)

                            gccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTc                                                                     <  1:661692/101‑1 (MQ=255)

                             ccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTc                                                                     <  1:661497/100‑1 (MQ=255)

                              cTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCAc                                                                   >  2:662046/1‑101 (MQ=255)

                              cTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCAc                                                                   <  1:661672/101‑1 (MQ=255)

                               tGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCAcc                                                                  <  2:661621/101‑1 (MQ=255)

                                ggTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCg                                                                 >  1:662048/1‑101 (MQ=255)

                                ggTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCg                                                                 <  1:661610/101‑1 (MQ=255)

                                 gTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCgg                                                                <  2:661505/101‑1 (MQ=255)

                                 gTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCgg                                                                <  1:661362/101‑1 (MQ=255)

                                   aaTCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGtt                                                              <  1:661514/101‑1 (MQ=255)

                                   aaTCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGtt                                                              <  1:661399/101‑1 (MQ=255)

                                   aaTCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGtt                                                              <  2:661745/101‑1 (MQ=255)

                                    aTCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCAc                                                                   >  2:662054/1‑95 (MQ=255)

                                    aTCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATACTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTa                                                             >  1:662051/1‑101 (MQ=255)

                                     tCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTAc                                                            >  1:662055/1‑101 (MQ=255)

                                       cTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCg                                                          >  1:662059/1‑101 (MQ=255)

                                        ttAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGc                                                         >  2:662061/1‑101 (MQ=255)

                                          aaaGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGcac                                                       <  1:661699/101‑1 (MQ=255)

                                          aaaGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGcac                                                       >  1:662064/1‑101 (MQ=255)

                                           aaGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGcaca                                                      <  1:661303/101‑1 (MQ=255)

                                             gCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGt                                                    >  2:662067/1‑101 (MQ=255)

                                              cGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTa                                                   >  1:662069/1‑101 (MQ=255)

                                              cGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTa                                                   <  2:661873/101‑1 (MQ=255)

                                               ggAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAg                                                  <  1:661598/101‑1 (MQ=255)

                                               ggAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAg                                                  >  2:662070/1‑101 (MQ=255)

                                                   aaTGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGa                                              >  2:662075/1‑101 (MQ=255)

                                                     tGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATc                                            >  1:662077/1‑101 (MQ=255)

                                                     tGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATc                                            <  1:661864/101‑1 (MQ=255)

                                                       gCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCAc                                          <  2:661710/101‑1 (MQ=255)

                                                         gATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGcgccg                                                                                >  2:662079/1‑61 (MQ=255)

                                                           tGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACcggc                                      <  1:661488/101‑1 (MQ=255)

                                                           tGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACcggc                                      >  2:662081/1‑101 (MQ=255)

                                                              cATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCggcggc                                   >  2:662086/1‑101 (MQ=255)

                                                              cATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCggcggc                                   >  1:662084/1‑101 (MQ=255)

                                                               aTAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGgcggcg                                  >  2:662088/1‑101 (MQ=255)

                                                               aTAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGgcggcg                                  <  2:661475/101‑1 (MQ=255)

                                                                 aGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCa                                >  2:662091/1‑101 (MQ=255)

                                                                 aGGAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCa                                >  1:662092/1‑101 (MQ=255)

                                                                   gAGTGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACcggc                                      <  2:661773/93‑1 (MQ=255)

                                                                      tGAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGcc                           <  1:661764/101‑1 (MQ=255)

                                                                       gAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCg                          >  1:662099/1‑101 (MQ=255)

                                                                       gAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCg                          >  2:662097/1‑101 (MQ=255)

                                                                       gAGGCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCg                          >  2:662100/1‑101 (MQ=255)

                                                                         ggCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGga                        <  2:661633/101‑1 (MQ=255)

                                                                          gCATCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGgag                       >  1:662103/1‑101 (MQ=255)

                                                                            aTCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCgg                         >  2:662106/1‑97 (MQ=255)

                                                                             tCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGc                                                         >  2:662109/1‑64 (MQ=255)

                                                                             tCTGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGc                    <  1:661607/101‑1 (MQ=255)

                                                                               tGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGa                  >  1:662110/1‑101 (MQ=255)

                                                                               tGCCACAGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGa                  <  2:661587/101‑1 (MQ=255)

                                                                                  cacaGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAAcc               <  1:661904/101‑1 (MQ=255)

                                                                                  cacaGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAAcc               >  1:662114/1‑101 (MQ=255)

                                                                                  cacaGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAAcc               >  2:662112/1‑101 (MQ=255)

                                                                                     aGGCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAAt            >  2:662117/1‑101 (MQ=255)

                                                                                      ggCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATc           >  1:662119/1‑101 (MQ=255)

                                                                                       gCTGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCa          <  1:661807/101‑1 (MQ=255)

                                                                                         tGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAAc        <  1:661743/101‑1 (MQ=255)

                                                                                         tGAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAAc        <  2:661662/101‑1 (MQ=255)

                                                                                          gAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCa          <  1:661668/98‑1 (MQ=255)

                                                                                          gAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACg       >  1:662123/1‑101 (MQ=255)

                                                                                          gAAGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACg       >  2:662124/1‑101 (MQ=255)

                                                                                           aaGATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGt      <  1:662006/101‑1 (MQ=255)

                                                                                             gATAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGc    >  1:662127/1‑101 (MQ=255)

                                                                                               tAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCGCCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCg  <  2:661540/101‑1 (MQ=255)

                                                                                               tAGTCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCg  <  2:661577/101‑1 (MQ=255)

TACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCG  >  NC_000913/447109‑447304

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