BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,332,505	A→G	Y159H (TAC→CAC)	b2229 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,332,505	0	A	G	78.9%	447.7 / 141.9	247	Y159H (TAC→CAC)	b2229	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (23/29); new base G (103/91); total (127/120)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

Only 100 of 253 total aligned reads displayed.

GTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCC  >  NC_000913/2332415‑2332596

gTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGa                                                                                           <  1:4252608/93‑1 (MQ=255)

gTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGa                                                                                           <  1:7230093/93‑1 (MQ=255)

 tATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGAt                                                                                          <  2:3865795/93‑1 (MQ=255)

  aTCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATg                                                                                         >  1:3393197/1‑93 (MQ=255)

  aTCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATg                                                                                         >  2:3577793/1‑93 (MQ=255)

  aTCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATg                                                                                         <  2:1955858/93‑1 (MQ=255)

   tCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGa                                                                                        <  1:2996220/93‑1 (MQ=255)

       tCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTTGATGACgga                                                                                    >  2:2087692/1‑91 (MQ=255)

       tCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTa                                                                                    >  2:2302336/1‑93 (MQ=255)

       tCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTa                                                                                    >  1:211913/1‑93 (MQ=255)

          ggTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACg                                                                                 >  2:6828545/1‑93 (MQ=255)

          ggTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACg                                                                                 >  1:6633111/1‑93 (MQ=255)

           gtgtCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGc                                                                                >  2:1829870/1‑93 (MQ=255)

            tgtCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGcc                                                                               >  1:2605190/1‑93 (MQ=255)

             gtCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGccc                                                                              >  1:6089163/1‑93 (MQ=255)

              tCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGcccc                                                                             >  2:2562798/1‑93 (MQ=255)

               ccTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCa                                                                            <  2:4089187/93‑1 (MQ=255)

                 ttCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATc                                                                          >  1:6333734/1‑93 (MQ=255)

                   ccATGTCATAAGTTGTTGCAGACTGACTGCGCGCAGTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCa                                                                        >  2:4526297/1‑93 (MQ=255)

                     aTGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGa                                                                      >  2:1391466/1‑93 (MQ=255)

                      tGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTGACGCCCCATCCAGAc                                                                     >  1:6694262/1‑93 (MQ=255)

                        tCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCa                                                                   >  2:4460824/1‑93 (MQ=255)

                        tCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCa                                                                   >  1:1040750/1‑93 (MQ=255)

                         cATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCAt                                                                  >  2:1736482/1‑93 (MQ=255)

                         cATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCAt                                                                  <  2:4048567/93‑1 (MQ=255)

                          aTAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCAtt                                                                 <  1:1744156/93‑1 (MQ=255)

                           tAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCCTTa                                                                >  1:324467/1‑93 (MQ=255)

                           tAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTa                                                                >  1:5892969/1‑93 (MQ=255)

                           tAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTa                                                                >  1:3778970/1‑93 (MQ=255)

                            aaGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTaa                                                               <  2:2645668/93‑1 (MQ=255)

                            aaGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTaa                                                               <  2:3207949/93‑1 (MQ=255)

                            aaGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTaa                                                               <  2:5230491/93‑1 (MQ=255)

                             aGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGATGTAACGCCCCATCCAGACCATTAAg                                                              <  1:2078172/93‑1 (MQ=255)

                              gttgttGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGt                                                             >  1:6260869/1‑93 (MQ=255)

                                tgttgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGc                                                           <  2:7189798/93‑1 (MQ=255)

                                  ttgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTg                                                         >  2:1349343/1‑93 (MQ=255)

                                  ttgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTg                                                         <  1:7168458/93‑1 (MQ=255)

                                   tgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTgg                                                        <  2:6078722/93‑1 (MQ=255)

                                   tgCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTgg                                                        <  1:5455183/93‑1 (MQ=255)

                                     cAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGc                                                      <  1:2390454/93‑1 (MQ=255)

                                      aGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCa                                                     <  1:475529/93‑1 (MQ=255)

                                      aGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCa                                                     <  2:5113473/93‑1 (MQ=255)

                                       gactgactGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCAt                                                    >  2:5298935/1‑93 (MQ=255)

                                       gactgactGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCAt                                                    <  2:3725963/93‑1 (MQ=255)

                                        actgactgCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATc                                                   <  2:5527249/93‑1 (MQ=255)

                                         ctgactgCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCg                                                  >  1:1302996/1‑93 (MQ=255)

                                          tgactgCGCGTATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGc                                                 >  1:6867426/1‑93 (MQ=255)

                                          tgactgCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGc                                                 >  2:3345156/1‑93 (MQ=255)

                                             ctgCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCt                                              <  1:1260234/93‑1 (MQ=255)

                                              tgCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTg                                             >  2:6657712/1‑93 (MQ=255)

                                               gcgcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGa                                            <  2:773581/93‑1 (MQ=255)

                                                cgcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGAt                                           <  1:4302392/93‑1 (MQ=255)

                                                 gcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATc                                          <  2:1177949/93‑1 (MQ=255)

                                                 gcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATc                                          <  2:3474728/93‑1 (MQ=255)

                                                 gcgcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATc                                          >  1:7430461/1‑93 (MQ=255)

                                                   gcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGa                                        >  2:4209672/1‑93 (MQ=255)

                                                   gcATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGa                                        >  1:3647145/1‑93 (MQ=255)

                                                    cATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaa                                       >  2:5712295/1‑93 (MQ=255)

                                                     aTTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGTCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaa                                      >  2:7258078/1‑93 (MQ=255)

                                                     aTTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaa                                      >  2:2781793/1‑93 (MQ=255)

                                                     aTTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaa                                      >  1:4662103/1‑93 (MQ=255)

                                                     aTTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaa                                      <  1:948195/93‑1 (MQ=255)

                                                     aTTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaa                                      >  1:499332/1‑93 (MQ=255)

                                                      ttCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaaa                                     <  1:5999293/93‑1 (MQ=255)

                                                      ttCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaaa                                     <  1:280390/93‑1 (MQ=255)

                                                       tCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaaaa                                    >  2:671956/1‑93 (MQ=255)

                                                         cGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaaaaaa                                  >  2:5922906/1‑93 (MQ=255)

                                                         cGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGaaaaaaa                                  <  1:3547191/93‑1 (MQ=255)

                                                          gTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAAt                                 >  2:1623474/1‑93 (MQ=255)

                                                          gTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAAt                                 >  2:3020932/1‑93 (MQ=255)

                                                          gTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAAt                                 >  2:839517/1‑93 (MQ=255)

                                                            tGTCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCa                               >  1:5748534/1‑93 (MQ=255)

                                                            tGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCa                               <  2:2427031/93‑1 (MQ=255)

                                                              tCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACTCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCata                             <  2:4938349/93‑1 (MQ=255)

                                                              tCAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCata                             <  2:1290247/93‑1 (MQ=255)

                                                              tCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCata                             <  1:507267/93‑1 (MQ=255)

                                                               cAGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCatat                            <  2:6845897/93‑1 (MQ=255)

                                                                aGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATc                           >  2:4879098/1‑93 (MQ=255)

                                                                aGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATc                           <  2:6558343/93‑1 (MQ=255)

                                                                aGTTTTCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATc                           <  2:3744814/93‑1 (MQ=255)

                                                                  ttttCGTGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGc                         >  1:3199827/1‑93 (MQ=255)

                                                                  ttttCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGc                         <  2:6133423/93‑1 (MQ=255)

                                                                     tCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAg                      >  1:5646548/1‑93 (MQ=255)

                                                                        tGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCa                   >  2:4976111/1‑93 (MQ=255)

                                                                        tGGCGCTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCa                   <  1:6794889/93‑1 (MQ=255)

                                                                            gcTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAgcgc               <  2:568237/93‑1 (MQ=255)

                                                                             cTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGcc              <  1:7057637/93‑1 (MQ=255)

                                                                             cTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGcc              <  2:4168825/93‑1 (MQ=255)

                                                                             cTTCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGcc              <  2:6225389/93‑1 (MQ=255)

                                                                              ttCCGGTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCt             <  2:4829145/93‑1 (MQ=255)

                                                                                 cGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGt          >  2:5740049/1‑93 (MQ=255)

                                                                                  ggTGTGGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGtt         >  2:7508177/1‑93 (MQ=255)

                                                                                    tgtgGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTa       <  1:5223076/93‑1 (MQ=255)

                                                                                    tgtgGTGGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTa       <  1:401233/93‑1 (MQ=255)

                                                                                     gtggtgGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTAt      >  1:1519719/1‑93 (MQ=255)

                                                                                     gtggtgGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTAt      <  1:1382312/93‑1 (MQ=255)

                                                                                     gtggtgGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTAt      >  1:7079348/1‑93 (MQ=255)

                                                                                      tggtggATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATg     <  2:5392807/93‑1 (MQ=255)

                                                                                      tggtggATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATg     <  2:2697459/93‑1 (MQ=255)

                                                                                        gTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTc   <  2:2814130/93‑1 (MQ=255)

GTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGGCGCTTCCGGTGTGGTAGATGACGTAACGCCCCATCCAGACCATTAAGTGCTGGGCATCGCCCTGATCGAAAAAAATCATATCGCCAGGCAGCGCCTGGTTTATGTCC  >  NC_000913/2332415‑2332596

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