BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	27,903	T→C	L204P (CTG→CCG)	yaaF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	27,903	0	T	C	99.7%	1149.8 / ‑6.2	354	L204P (CTG→CCG)	yaaF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (188/165); total (189/165)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 380 total aligned reads displayed.

CCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGA  >  NC_000913/27813‑27993

ccTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGc                                                                                          >  1:7777412/1‑93 (MQ=255)

ccTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGc                                                                                          <  1:1678730/93‑1 (MQ=255)

 cTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGcc                                                                                         >  1:419481/1‑93 (MQ=255)

   gtgtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACcgccgc                                                                                       >  2:8745490/1‑93 (MQ=255)

   gtgtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACcgccgc                                                                                       >  2:7337494/1‑93 (MQ=255)

    tgtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCa                                                                                      >  1:1001054/1‑93 (MQ=255)

        ttCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTa                                                                                  <  1:8448684/93‑1 (MQ=255)

        ttCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTa                                                                                  >  2:3665015/1‑93 (MQ=255)

         tCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTaa                                                                                 <  1:9179341/93‑1 (MQ=255)

         tCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTaa                                                                                 >  2:4777190/1‑93 (MQ=255)

         tCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTaa                                                                                 >  1:3582259/1‑93 (MQ=255)

          ccGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTaaa                                                                                >  2:3487127/1‑93 (MQ=255)

            gCATTGGAATTGATATCTTCATGTGCGGTTTGGATGGCACCAATCAGGCAATACTTACTCCTGTCGATCTCGCGACATCGCCTCAGTTAAAcc                                                                              >  2:9244841/1‑93 (MQ=255)

            gCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAcc                                                                              >  1:1281271/1‑93 (MQ=255)

            gCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATAAGGCAATATTAACTCCTGACTATCTCTCTACAGCGCCGCAGTTAAAcc                                                                              >  2:4427258/1‑93 (MQ=255)

            gCAGTAGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAcc                                                                              <  1:2079795/93‑1 (MQ=255)

             cAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCg                                                                             >  2:3059763/1‑93 (MQ=255)

              aGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGt                                                                            <  1:6514726/93‑1 (MQ=255)

                tGGTATTGAAATCGTCATGTGCGGTTTGTATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTAc                                                                          <  1:2599533/93‑1 (MQ=255)

                tGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTAc                                                                          <  1:3844290/93‑1 (MQ=255)

                tGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTAc                                                                          <  1:2400461/93‑1 (MQ=255)

                 ggTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTAcc                                                                         <  2:266434/93‑1 (MQ=255)

                   tATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCgg                                                                       >  2:6329251/1‑93 (MQ=255)

                     ttGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGa                                                                     >  2:5365055/1‑93 (MQ=255)

                     ttGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGa                                                                     <  1:7139885/93‑1 (MQ=255)

                       gAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACAGTACCGGGaaa                                                                   >  2:4674460/1‑93 (MQ=255)

                        aaaTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGaaaa                                                                  <  2:4026697/93‑1 (MQ=255)

                        aaaTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGaaaa                                                                  <  1:4483939/93‑1 (MQ=255)

                          aTCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATg                                                                >  2:3125034/1‑93 (MQ=255)

                             gTCATGTGCGGTTTTGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCtt                                                             <  2:8636793/93‑1 (MQ=255)

                              tCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTc                                                            <  2:778888/93‑1 (MQ=255)

                               cATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCa                                                           >  1:202215/1‑93 (MQ=255)

                               cATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCa                                                           >  2:8996322/1‑93 (MQ=255)

                               cATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCa                                                           >  2:4255027/1‑93 (MQ=255)

                                aTGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCAc                                                          <  2:8476566/93‑1 (MQ=255)

                                 tgtgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACg                                                         >  2:8943958/1‑93 (MQ=255)

                                 tgtgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACg                                                         >  2:9753603/1‑93 (MQ=255)

                                  gtgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGc                                                        <  1:50072/93‑1 (MQ=255)

                                   tgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGcc                                                       <  1:5629015/93‑1 (MQ=255)

                                   tgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGcc                                                       <  1:5472715/93‑1 (MQ=255)

                                   tgCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGcc                                                       >  2:9070049/1‑93 (MQ=255)

                                    gCGGTTTGGATGTCACCAATCAGGCAATATTACCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGccc                                                      <  1:4644590/93‑1 (MQ=255)

                                    gCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGccc                                                      <  2:7166188/93‑1 (MQ=255)

                                      ggTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCCCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGGTTCACGCCCTg                                                    <  2:5861498/93‑1 (MQ=255)

                                       gTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGt                                                   >  2:3905848/1‑93 (MQ=255)

                                       gTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGt                                                   >  2:3384560/1‑93 (MQ=255)

                                       gTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGt                                                   >  2:5834364/1‑93 (MQ=255)

                                        tttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGtt                                                  >  1:2520211/1‑93 (MQ=255)

                                        tttGGATGTCACCAATCAGGCAATATTAACTCCTGACCATCTCTCTACGCCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGtt                                                  <  1:8146303/93‑1 (MQ=255)

                                         ttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGttt                                                 >  1:2979608/1‑93 (MQ=255)

                                          tGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTa                                                <  1:7547021/93‑1 (MQ=255)

                                             aTGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGcc                                             >  1:6017610/1‑93 (MQ=255)

                                              tGTCACCAATAAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCa                                            >  2:1720618/1‑93 (MQ=255)

                                               gTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCAc                                           >  2:5242531/1‑93 (MQ=255)

                                                 cACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTa                                         <  1:6208983/93‑1 (MQ=255)

                                                  cccAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAc                                        <  1:1899214/92‑1 (MQ=255)

                                                  aCCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAc                                        >  2:1443419/1‑93 (MQ=255)

                                                   ccAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGATTAGCCACTAcc                                       >  1:6757851/1‑93 (MQ=255)

                                                     aaTCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGt                                     >  2:4615785/1‑93 (MQ=255)

                                                      aTCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTa                                    >  2:2274500/1‑93 (MQ=255)

                                                       tCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAg                                   <  2:961311/93‑1 (MQ=255)

                                                         aGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCg                                 >  2:2937717/1‑93 (MQ=255)

                                                          ggCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCgg                                >  2:7473248/1‑93 (MQ=255)

                                                          ggCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCgg                                >  1:7377339/1‑93 (MQ=255)

                                                           gCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGc                               >  1:5199532/1‑93 (MQ=255)

                                                            cAATATTACCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCa                              <  1:2781233/93‑1 (MQ=255)

                                                            cAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCa                              >  2:8665537/1‑93 (MQ=255)

                                                            cAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCa                              >  1:4150291/1‑93 (MQ=255)

                                                              atatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGt                            <  1:9047516/93‑1 (MQ=255)

                                                               tatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTa                           >  2:3621565/1‑93 (MQ=255)

                                                               tatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTa                           <  2:7909017/93‑1 (MQ=255)

                                                               tatTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTa                           >  2:5848111/1‑93 (MQ=255)

                                                                  tAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGc                        >  2:6859473/1‑93 (MQ=255)

                                                                   aaCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCa                       <  2:7567958/93‑1 (MQ=255)

                                                                    actactGAATATATCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaa                      <  1:6403681/93‑1 (MQ=255)

                                                                    aCTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaa                      <  2:7707520/93‑1 (MQ=255)

                                                                     cTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaaa                     <  1:5882028/93‑1 (MQ=255)

                                                                     cTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaaa                     <  2:9435500/93‑1 (MQ=255)

                                                                     cTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCaaa                     >  2:5616337/1‑93 (MQ=255)

                                                                      cccTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAg                    <  2:1871168/92‑1 (MQ=255)

                                                                         tGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCgg                 >  1:9439708/1‑93 (MQ=255)

                                                                          gACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGc                <  2:9338225/93‑1 (MQ=255)

                                                                           aCTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCt               <  2:3811234/93‑1 (MQ=255)

                                                                           aCTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATCCAAAGCGGCt               <  2:6549372/93‑1 (MQ=255)

                                                                            cTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt              >  1:3753303/1‑93 (MQ=255)

                                                                             tATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTg             <  2:4086418/93‑1 (MQ=255)

                                                                             tATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTg             <  1:3220814/93‑1 (MQ=255)

                                                                             tATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGAATGCAAAGCGGtttt             >  1:8666833/1‑92 (MQ=255)

                                                                               tctctcTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCg           <  1:1528392/93‑1 (MQ=255)

                                                                               tctctcTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCg           >  2:1026169/1‑93 (MQ=255)

                                                                               tctctcTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCg           <  2:207547/93‑1 (MQ=255)

                                                                                ctctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTTCgt          >  1:9784017/1‑92 (MQ=255)

                                                                                 tctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGaa         <  2:2455880/93‑1 (MQ=255)

                                                                                 cctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGaa         <  2:5566354/92‑1 (MQ=255)

                                                                                  ctctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAAt        <  2:1261290/93‑1 (MQ=255)

                                                                                    ctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGc      <  2:8253050/93‑1 (MQ=255)

                                                                                     tACACCGCCGCAGTTAAACCGTACCGGGATAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTTCGaacggc     >  1:4832226/1‑89 (MQ=255)

                                                                                     tACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCa     >  1:3327760/1‑93 (MQ=255)

                                                                                     tACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCa     >  1:6296336/1‑93 (MQ=255)

                                                                                        acTGCCGCAGTTAAACCGTATCGGGAAAATTCTTCACGCCCTGTTTAGCCACTACCGTAGTGGTAGAATGCAAAGCGGCTTGCGAAGGCACGa  >  1:8947910/1‑93 (MQ=255)

CCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGA  >  NC_000913/27813‑27993

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