BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	368,655	C→T	S15S (AGC→AGT)	mhpA →	putative 3‑(3‑hydroxyphenyl)propanoate/3‑hydroxycinnamate hydroxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	368,655	0	C	T	99.7%	1435.3 / ‑5.7	338	S15S (AGC→AGT)	mhpA	putative 3‑(3‑hydroxyphenyl)propanoate/3‑hydroxycinnamate hydroxylase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (161/176); minor base G (1/0); total (162/176)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

GTTGCGCTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGCGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTC  >  NC_000913/368557‑368751

 ttGCGCTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGtt                                                                                               >  2:683225/1‑101 (MQ=255)

 ttGCGCTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGtt                                                                                               <  2:682585/101‑1 (MQ=255)

  tGCGCTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTc                                                                                              <  2:682774/101‑1 (MQ=255)

  tGCGCTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTc                                                                                              <  2:682758/101‑1 (MQ=255)

   gcgcTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACCCCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCa                                                                                             <  1:682752/101‑1 (MQ=255)

   gcgcTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCa                                                                                             >  1:683228/1‑101 (MQ=255)

    cgcTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAg                                                                                            >  2:683231/1‑101 (MQ=255)

    cgcTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAg                                                                                            <  1:682656/101‑1 (MQ=255)

    cgcTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAg                                                                                            >  1:683232/1‑101 (MQ=255)

     gcTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAgg                                                                                           <  1:682371/101‑1 (MQ=255)

     gcTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAgg                                                                                           >  1:683236/1‑101 (MQ=255)

      cTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGt                                                                                          <  2:682714/101‑1 (MQ=255)

      cTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGt                                                                                          >  1:683241/1‑101 (MQ=255)

       tGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTg                                                                                         >  2:683243/1‑101 (MQ=255)

       tGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTg                                                                                         >  1:683247/1‑101 (MQ=255)

         tACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGc                                                                                       <  1:682808/101‑1 (MQ=255)

         tACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGc                                                                                       <  2:682551/101‑1 (MQ=255)

         tACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGc                                                                                       >  2:683250/1‑101 (MQ=255)

            tgtgCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGAt                                                                                    <  2:682602/101‑1 (MQ=255)

            tgtgCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGAt                                                                                    >  1:683254/1‑101 (MQ=255)

             gtgCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATc                                                                                   <  1:682745/101‑1 (MQ=255)

                cgcgcAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCt                                                                                >  1:683258/1‑101 (MQ=255)

                 gcgcAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCACTg                                                                               <  1:682685/101‑1 (MQ=255)

                  cgcAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTgg                                                                              >  1:683263/1‑101 (MQ=255)

                   gcAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGt                                                                             <  2:682678/101‑1 (MQ=255)

                   gcAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGt                                                                             >  2:683264/1‑101 (MQ=255)

                    cAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCt                                                                                <  1:682630/97‑1 (MQ=255)

                       cGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCg                                                                         <  1:682981/101‑1 (MQ=255)

                        gACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCgg                                                                        <  2:682720/101‑1 (MQ=255)

                        gACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCgg                                                                        >  1:683268/1‑101 (MQ=255)

                           aTTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGccc                                                                     >  2:683270/1‑101 (MQ=255)

                           aTTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGccc                                                                     <  1:682663/101‑1 (MQ=255)

                            ttttGTCCGAGTCGTGAGGGACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCg                                                                    >  1:683274/1‑101 (MQ=255)

                               tGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGtt                                                                 <  2:682989/101‑1 (MQ=255)

                                gTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTg                                                                <  2:682528/101‑1 (MQ=255)

                                  ccGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGtt                                                                 >  1:683279/1‑98 (MQ=255)

                                   cGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGc                                                             >  2:683282/1‑101 (MQ=255)

                                    gAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCt                                                            >  1:683283/1‑101 (MQ=255)

                                    gAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCt                                                            <  2:682964/101‑1 (MQ=255)

                                     aGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCtg                                                           <  1:682729/101‑1 (MQ=255)

                                      gTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCtga                                                          >  2:683288/1‑101 (MQ=255)

                                       tCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCtgat                                                         <  2:682733/101‑1 (MQ=255)

                                       tCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCtgat                                                         <  2:682617/101‑1 (MQ=255)

                                        cGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCtgatg                                                        >  2:683292/1‑101 (MQ=255)

                                          tGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTgatgat                                                      >  1:683299/1‑101 (MQ=255)

                                          tGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTgatgat                                                      >  2:683297/1‑101 (MQ=255)

                                           gAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGatgatg                                                     <  2:682804/101‑1 (MQ=255)

                                           gAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGatgatg                                                     <  2:683149/101‑1 (MQ=255)

                                            aGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATgg                                                    >  1:683304/1‑101 (MQ=255)

                                            aGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATgg                                                    >  2:683303/1‑101 (MQ=255)

                                             ggTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGc                                                   <  1:682817/101‑1 (MQ=255)

                                             ggTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGc                                                   <  2:682376/101‑1 (MQ=255)

                                             ggTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGc                                                   <  2:682481/101‑1 (MQ=255)

                                              gTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCg                                                  >  2:683310/1‑101 (MQ=255)

                                               tACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGa                                                 <  1:682791/101‑1 (MQ=255)

                                               tACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGa                                                 >  1:683311/1‑101 (MQ=255)

                                                aCTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCg                                                                    <  2:682741/81‑1 (MQ=255)

                                                  tGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACt                                              >  2:683319/1‑101 (MQ=255)

                                                  tGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACt                                              >  1:683315/1‑101 (MQ=255)

                                                   gAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTa                                             <  1:682953/101‑1 (MQ=255)

                                                    aaaTGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTAt                                            >  2:683322/1‑101 (MQ=255)

                                                    aaaTGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTAt                                            <  2:682687/101‑1 (MQ=255)

                                                     aaTGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTAtc                                           >  1:683323/1‑101 (MQ=255)

                                                         gCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCgg                                       <  1:682504/101‑1 (MQ=255)

                                                           aaTACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGcc                                     <  1:682595/101‑1 (MQ=255)

                                                           aaTACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGcc                                     >  2:683329/1‑101 (MQ=255)

                                                            aTACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCa                                    <  2:682788/101‑1 (MQ=255)

                                                            aTACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCa                                    <  1:682768/101‑1 (MQ=255)

                                                            aTACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCa                                    <  2:682653/101‑1 (MQ=255)

                                                             tACAACACCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAg                                   <  2:682660/101‑1 (MQ=255)

                                                              acaacaCCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGa                                  <  2:682692/101‑1 (MQ=255)

                                                               caacaCCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGAt                                 >  1:683331/1‑101 (MQ=255)

                                                                aacaCCCTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTg                                                                >  1:683334/1‑69 (MQ=255)

                                                                    cccTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCa                            >  1:683339/1‑101 (MQ=255)

                                                                      cTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCAtt                          <  2:682858/101‑1 (MQ=255)

                                                                      cTGACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCAtt                          >  2:683341/1‑101 (MQ=255)

                                                                        gACATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGa                        >  2:683346/1‑101 (MQ=255)

                                                                         aCATCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGAc                       >  1:683348/1‑101 (MQ=255)

                                                                           aTCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGt                     >  1:683353/1‑101 (MQ=255)

                                                                            tCCAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTg                    >  1:683356/1‑101 (MQ=255)

                                                                             ccAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGc                   >  2:683361/1‑101 (MQ=255)

                                                                             ccAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGc                   <  1:682455/101‑1 (MQ=255)

                                                                             ccAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGc                   >  1:683360/1‑101 (MQ=255)

                                                                              cAGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCt                  <  1:682773/101‑1 (MQ=255)

                                                                               aGCCTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCtg                 <  2:682889/101‑1 (MQ=255)

                                                                                 ccTGCTGTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCtggt               <  2:682873/101‑1 (MQ=255)

                                                                                  ctgctgTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCtggtg              <  1:682844/101‑1 (MQ=255)

                                                                                  ctgctgTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCtggtg              >  1:683366/1‑101 (MQ=255)

                                                                                   tgctgTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCtggtgg             <  1:683163/101‑1 (MQ=255)

                                                                                   tgctgTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCtggtgg             <  2:682839/101‑1 (MQ=255)

                                                                                    gctgTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCtggtg              >  2:683369/1‑99 (MQ=255)

                                                                                      tgTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGccc                                                                     <  1:682712/42‑1 (MQ=255)

                                                                                       gTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGga         <  1:682828/101‑1 (MQ=255)

                                                                                       gTTAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGga         >  2:683374/1‑101 (MQ=255)

                                                                                        ttAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGgag        >  1:683376/1‑101 (MQ=255)

                                                                                        ttAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGgag        <  1:682872/101‑1 (MQ=255)

                                                                                        ttAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGgag        >  2:683378/1‑101 (MQ=255)

                                                                                         tAACCATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGgaga       <  2:682825/101‑1 (MQ=255)

                                                                                            ccATAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAAc    >  1:683384/1‑101 (MQ=255)

                                                                                              aTAGTGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTc  <  2:683012/101‑1 (MQ=255)

GTTGCGCTGTACTGTGCGCGCAACGACATTTTGTCCGAGTCGTGAGGTACTGAAATGGCAATACAACACCCTGACATCCAGCCTGCTGTTAACCATAGCGTTCAGGTGGCGATCGCTGGTGCCGGCCCGGTTGGGCTGATGATGGCGAACTATCTCGGCCAGATGGGCATTGACGTGCTGGTGGTGGAGAAACTC  >  NC_000913/368557‑368751

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: `‑`