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	98.8%	532.0 / ‑2.6	166	L204P (CTG→CCG)	yaaF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (94/70); total (95/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 183 total aligned reads displayed.

TGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATC  >  NC_000913/27811‑27995

tgcCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACt                                                                                              >  2:4015779/1‑93 (MQ=255)

    tgtgTCTTCCTCAGTGGTATTGAAATCGTTATGTGCAGTTTTGATGTCACCAATCAGGCAATATTAACGTCTGATTATCTATCTACACTGctg                                                                                          >  2:3246858/1‑91 (MQ=255)

    tgtgTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACAccgccg                                                                                          >  2:961381/1‑93 (MQ=255)

    tgtgTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACAccgccg                                                                                          >  2:648254/1‑93 (MQ=255)

    tgtgTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACAccgccg                                                                                          >  1:2945530/1‑93 (MQ=255)

    tgtgTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACAccgccg                                                                                          >  1:2593687/1‑93 (MQ=255)

    tgtgTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACAccgccg                                                                                          >  2:3246856/1‑93 (MQ=255)

     gtgtCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACcgccgc                                                                                         <  1:3679848/93‑1 (MQ=255)

        tCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGt                                                                                      >  1:258670/1‑93 (MQ=255)

        tCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGt                                                                                      <  2:1332864/93‑1 (MQ=255)

         cTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGtt                                                                                     <  2:992486/93‑1 (MQ=255)

             cGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAc                                                                                 <  1:3203150/93‑1 (MQ=255)

              gCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAAcc                                                                                <  1:4606828/93‑1 (MQ=255)

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

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

                 gTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACTAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTa                                                                             >  1:3806404/1‑93 (MQ=255)

                 gTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTa                                                                             <  1:4443554/93‑1 (MQ=255)

                   ggTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTAcc                                                                           <  1:3289981/93‑1 (MQ=255)

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

                   ggTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTAcc                                                                           >  1:3650214/1‑93 (MQ=255)

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

                    gTATTGAAATCGTCATGTGCGTTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTTCCg                                                                          <  2:1809747/93‑1 (MQ=255)

                         gAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGaaa                                                                     >  1:2361246/1‑93 (MQ=255)

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

                             tCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGc                                                                 >  2:1622215/1‑93 (MQ=255)

                               gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCtt                                                               >  1:1255778/1‑93 (MQ=255)

                               gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCtt                                                               >  2:1529615/1‑93 (MQ=255)

                               ctCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCtt                                                               <  2:633732/92‑1 (MQ=255)

                                tCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTc                                                              >  2:1005132/1‑93 (MQ=255)

                                tCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTc                                                              >  2:4376768/1‑93 (MQ=255)

                                tCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTc                                                              >  1:3662368/1‑93 (MQ=255)

                                 cATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCttga                                                             <  1:668891/93‑3 (MQ=255)

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

                                 cATGTGCGGTTTGCATGTCACCAATCAGGCAATATTAACTCCTGACTAGCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCa                                                             <  2:1503399/93‑1 (MQ=255)

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

                                  aTGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTCAAACCGTACCGGGAAAATGCTTCAc                                                            >  2:3734683/1‑93 (MQ=255)

                                  aTGTGCGGTTTCGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGCGAAAATGCTTCAc                                                            <  1:3136111/93‑1 (MQ=255)

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

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

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

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

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

                                       cGGTTTGGATGTCACCAATCAGACAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCt                                                       <  1:2463401/93‑1 (MQ=255)

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

                                           ttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGttt                                                   <  2:3437386/93‑1 (MQ=255)

                                             ggATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAg                                                 >  1:123541/1‑93 (MQ=255)

                                               aTGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGcc                                               <  1:1981567/93‑1 (MQ=255)

                                                tGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCa                                              >  2:2453165/1‑93 (MQ=255)

                                                tGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCa                                              <  1:3835174/93‑1 (MQ=255)

                                                 ttCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCAc                                             <  1:3508742/92‑1 (MQ=255)

                                                 gTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCAc                                             <  1:1024801/93‑1 (MQ=255)

                                                 gTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCAc                                             <  2:26875/93‑1 (MQ=255)

                                                 gTCACCAATCAGGCAATATTAACCCCTGACTATCTCTCTACACCGCCGCAGTTAAACTGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCAc                                             <  2:3048156/93‑1 (MQ=255)

                                                  tCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACt                                            >  2:456978/1‑93 (MQ=255)

                                                  tCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACt                                            >  1:3626507/1‑93 (MQ=255)

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

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

                                                     ccAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAcc                                         <  2:506738/93‑1 (MQ=255)

                                                     ccAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAcc                                         <  1:4563440/93‑1 (MQ=255)

                                                      cAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCg                                        >  1:895776/1‑93 (MQ=255)

                                                      cAATCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCg                                        >  1:2062383/1‑93 (MQ=255)

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

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

                                                        aTCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTa                                      <  1:4277012/93‑1 (MQ=255)

                                                         tCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAg                                     >  1:3169308/1‑93 (MQ=255)

                                                         tCAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAg                                     >  1:562155/1‑93 (MQ=255)

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

                                                          cAGGCAATATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGc                                    <  2:816119/93‑1 (MQ=255)

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

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

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

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

                                                               aaTATTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAg                               >  1:2222834/1‑93 (MQ=255)

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

                                                                  atTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTAt                            <  2:765223/93‑1 (MQ=255)

                                                                  atTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTAt                            <  2:1801447/93‑1 (MQ=255)

                                                                  atTAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTAt                            >  2:352675/1‑93 (MQ=255)

                                                                   ttAACTCCTGACTATCTCTCTACTCCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATg                           <  2:4227392/93‑1 (MQ=255)

                                                                   ttAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATg                           >  2:4157012/1‑93 (MQ=255)

                                                                    tAACTCCTGACTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGc                          <  1:4153420/93‑1 (MQ=255)

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

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

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

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

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

                                                                             aCTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCt                 >  1:731488/1‑93 (MQ=255)

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

                                                                              cTATCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtt                <  1:4289665/93‑1 (MQ=255)

                                                                                aTCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCtttc              >  2:1918864/1‑91 (MQ=255)

                                                                                aTCTCTCTACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGc              <  1:1047205/93‑1 (MQ=255)

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

                                                                                     tctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTTCggcag         >  2:2614337/1‑89 (MQ=255)

                                                                                     tctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGGGAATg         >  1:1603867/1‑93 (MQ=255)

                                                                                     tctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGaaag         >  1:1759493/1‑91 (MQ=255)

                                                                                     tctACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATg         <  2:2830340/93‑1 (MQ=255)

                                                                                      ctACACCGCCGCAGTTAGACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGc        <  1:3416690/93‑1 (MQ=255)

                                                                                       tACACCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCATGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCa       <  1:4105508/93‑1 (MQ=255)

                                                                                        acacCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTAGCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCAc      <  1:2021289/93‑1 (MQ=255)

                                                                                        acacCGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTTCGAATGCAc      >  2:3115624/1‑93 (MQ=255)

                                                                                            tGCCGCTGTTAATCCGTATCTGGAAAATGCTTCACGTCCTTTTTAGCCACTATTTTAGCGGTAGTTTGCTAAGCGGCTTTCGTATGCACGAtt  >  2:4331209/1‑92 (MQ=255)

TGCCTGTGTCTTCCGCAGTGGTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATC  >  NC_000913/27811‑27995

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