BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,804,546	T→C	E50G (GAG→GGG)	rfaG ←	glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,804,546	0	T	C	99.6%	800.0 / ‑4.9	242	E50G (GAG→GGG)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (117/124); total (117/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 244 total aligned reads displayed.

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

ggATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCa                                                                                         >  1:1130920/1‑93 (MQ=255)

ggATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCa                                                                                         >  1:2253305/1‑93 (MQ=255)

ggATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCa                                                                                         <  2:4872529/93‑1 (MQ=255)

   tGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAAt                                                                                      <  1:6302660/93‑1 (MQ=255)

     ctctTTGAGATCATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGc                                                                                    <  1:5743353/93‑1 (MQ=255)

      tctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCt                                                                                   >  1:3485404/1‑93 (MQ=255)

      tctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCt                                                                                   <  1:290949/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  1:677599/93‑1 (MQ=255)

         ttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTc                                                                                <  2:1101816/93‑1 (MQ=255)

          tGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCg                                                                               <  2:647976/93‑1 (MQ=255)

          tGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCg                                                                               <  1:6083377/93‑1 (MQ=255)

          tGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCg                                                                               <  1:3224799/93‑1 (MQ=255)

           gagaTGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCgg                                                                              <  2:669569/93‑1 (MQ=255)

             gatgatTTTGTACCCAGGCATAATATTCTGCATTGCGTGCATGGTTGGGATGGGACTTAACTCGCACCTGAATAAGCCCAAATGCTTTCGGGc                                                                            >  2:21453/1‑93 (MQ=255)

             gatgatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGc                                                                            >  1:961509/1‑93 (MQ=255)

             gatgatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGc                                                                            >  2:6163334/1‑93 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         >  2:4816305/1‑93 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         >  1:6519786/1‑93 (MQ=255)

                 atTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATc                                                                        >  1:2325871/1‑93 (MQ=255)

                 atTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATc                                                                        <  1:3146007/93‑1 (MQ=255)

                   tttGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGc                                                                      <  1:919909/93‑1 (MQ=255)

                    ttGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGcc                                                                     <  2:3796834/93‑1 (MQ=255)

                     tGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCt                                                                    <  2:3059893/93‑1 (MQ=255)

                     tGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCt                                                                    <  2:6312402/93‑1 (MQ=255)

                      gTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCtt                                                                   >  2:2231723/1‑93 (MQ=255)

                      gTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCtt                                                                   <  1:2002444/93‑1 (MQ=255)

                       tACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTc                                                                  <  2:5630399/93‑1 (MQ=255)

                       tACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTc                                                                  <  2:6604268/93‑1 (MQ=255)

                        aCCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTcc                                                                 >  2:861031/1‑93 (MQ=255)

                        aCCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTcc                                                                 <  2:3534009/93‑1 (MQ=255)

                        aCCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTcc                                                                 >  1:5167295/1‑93 (MQ=255)

                          ccAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCa                                                               <  2:1549354/93‑1 (MQ=255)

                           cAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCAc                                                              >  1:1472568/1‑93 (MQ=255)

                            aGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACg                                                             <  2:3440025/93‑1 (MQ=255)

                            aGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACg                                                             >  1:3499921/1‑93 (MQ=255)

                             ggCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGa                                                            >  2:1470535/1‑93 (MQ=255)

                             ggCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGa                                                            >  1:3954577/1‑93 (MQ=255)

                               cATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACt                                                          <  1:5642057/93‑1 (MQ=255)

                                 taataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtgt                                                        >  2:5492618/1‑93 (MQ=255)

                                  aataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtgtg                                                       >  1:2153878/1‑93 (MQ=255)

                                     atTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtat                                                    >  2:2858460/1‑93 (MQ=255)

                                     atTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtat                                                    >  2:1710711/1‑93 (MQ=255)

                                      ttCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtata                                                   <  2:5062454/93‑1 (MQ=255)

                                      ttCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtata                                                   >  2:4490337/1‑93 (MQ=255)

                                       tCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatat                                                  >  1:1896521/1‑93 (MQ=255)

                                        cTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatata                                                 <  1:2800185/93‑1 (MQ=255)

                                          gCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAcc                                               >  2:2713396/1‑93 (MQ=255)

                                           cATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAccc                                              <  2:5107926/93‑1 (MQ=255)

                                            aTTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCg                                             >  2:3111554/1‑93 (MQ=255)

                                               gCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAAc                                          <  1:1280651/93‑1 (MQ=255)

                                                 gTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAt                                        <  2:3923509/93‑1 (MQ=255)

                                                 gTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAt                                        <  2:6363537/93‑1 (MQ=255)

                                                    cATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggt                                     >  2:2066245/1‑93 (MQ=255)

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    >  1:2523839/1‑93 (MQ=255)

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    >  2:5083268/1‑93 (MQ=255)

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    <  1:1025369/93‑1 (MQ=255)

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    <  1:3308223/93‑1 (MQ=255)

                                                      tggttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtgg                                   >  2:4604685/1‑93 (MQ=255)

                                                      tggttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtgg                                   >  2:3237781/1‑93 (MQ=255)

                                                        gttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGcc                                 >  1:448621/1‑93 (MQ=255)

                                                        gttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGcc                                 >  2:6706518/1‑93 (MQ=255)

                                                         ttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGccc                                <  2:2407620/93‑1 (MQ=255)

                                                              aTGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGc                           >  2:1951431/1‑93 (MQ=255)

                                                               tGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCg                          >  2:5760199/1‑93 (MQ=255)

                                                               tGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCg                          >  2:365823/1‑93 (MQ=255)

                                                                gggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCgg                         >  2:3614545/1‑93 (MQ=255)

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGc                        <  1:3687849/93‑1 (MQ=255)

                                                                 ggACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGc                        <  2:6030394/93‑1 (MQ=255)

                                                                  gACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCa                       >  2:3114604/1‑93 (MQ=255)

                                                                    cTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAAc                     >  1:6041190/1‑93 (MQ=255)

                                                                    cTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTAGATACCGGAACATGGTGGCCCCGTGCGGCAAc                     <  1:3094879/93‑1 (MQ=255)

                                                                     ttAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACt                    <  1:4390387/93‑1 (MQ=255)

                                                                     ttAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACt                    >  1:468001/1‑93 (MQ=255)

                                                                     ttAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACt                    >  2:2779382/1‑93 (MQ=255)

                                                                      tAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTg                   <  1:6231511/93‑1 (MQ=255)

                                                                       aaCTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGt                  <  1:2396679/93‑1 (MQ=255)

                                                                       aaCTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGt                  <  1:1428740/93‑1 (MQ=255)

                                                                       aaCTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGt                  <  1:4832638/93‑1 (MQ=255)

                                                                        aCTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGtt                 >  2:4170825/1‑93 (MQ=255)

                                                                        aCTGGAACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGtt                 <  1:2096840/93‑1 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGa               <  2:2934207/93‑1 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGa               >  2:5001498/1‑93 (MQ=255)

                                                                           ggCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGAt              <  1:3637993/93‑1 (MQ=255)

                                                                            gCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATg             >  2:2346270/1‑93 (MQ=255)

                                                                            gCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATg             <  1:1724524/93‑1 (MQ=255)

                                                                             cACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGc            >  2:824499/1‑93 (MQ=255)

                                                                             cACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGc            <  1:3079920/93‑1 (MQ=255)

                                                                             cACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGc            <  2:1648464/93‑1 (MQ=255)

                                                                              aCCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCa           <  1:5279038/93‑1 (MQ=255)

                                                                              aCCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCa           >  1:4246546/1‑93 (MQ=255)

                                                                               ccTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCaa          <  2:2037117/93‑1 (MQ=255)

                                                                               ccTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCaa          <  2:1389189/93‑1 (MQ=255)

                                                                                cTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAAt         <  2:2686106/93‑1 (MQ=255)

                                                                                  gAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgc       <  1:4095611/93‑1 (MQ=255)

                                                                                    aTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgcgc     <  1:1721570/93‑1 (MQ=255)

                                                                                    aTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgcgc     >  1:477856/1‑93 (MQ=255)

                                                                                    aTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgcgc     <  1:3387610/93‑1 (MQ=255)

                                                                                      aaGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCAt   <  2:5305686/93‑1 (MQ=255)

                                                                                      aaGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCAt   >  2:437454/1‑93 (MQ=255)

                                                                                       aGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATa  >  2:4099488/1‑93 (MQ=255)

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

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