BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,462,072	A→G	E170E (GAA→GAG)	b1399 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,462,072	0	A	G	99.6%	689.9 / NA	223	E170E (GAA→GAG)	b1399	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (107/114); total (109/114)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 229 total aligned reads displayed.

CAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAA  >  NC_000913/1461982‑1462162

cAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGgc                                                                                          >  2:2148954/1‑93 (MQ=255)

  gCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGgcgc                                                                                        >  2:4915066/1‑93 (MQ=255)

    cTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGCCACTTTTGCATGAAGCGGGGGTGGCGGATAACGTGATTTGTTTTGTGGCGCaa                                                                                      >  2:6023082/1‑93 (MQ=255)

     tGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCaaa                                                                                     <  2:1856500/93‑1 (MQ=255)

      gATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAAt                                                                                    >  1:370381/1‑93 (MQ=255)

       aTGGCATCGCCGTCGCAAAAACTGGCTGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATa                                                                                   <  1:2440080/93‑1 (MQ=255)

       aTGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATa                                                                                   >  2:353522/1‑93 (MQ=255)

        tGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATAc                                                                                  <  2:2370039/93‑1 (MQ=255)

        tGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATAc                                                                                  <  2:5695990/93‑1 (MQ=255)

         ggCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATAcc                                                                                 >  1:2232952/1‑93 (MQ=255)

          gCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCa                                                                                <  2:3686846/93‑1 (MQ=255)

            aTCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACt                                                                              <  1:516522/93‑1 (MQ=255)

            aTCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACt                                                                              >  2:3420537/1‑93 (MQ=255)

            aTCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACt                                                                              >  1:1180481/1‑93 (MQ=255)

            aTCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACt                                                                              >  2:1711958/1‑93 (MQ=255)

             tCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATAACACTg                                                                             >  1:695855/1‑93 (MQ=255)

              cGCCGTCGTAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTgg                                                                            <  1:4203721/93‑1 (MQ=255)

              cGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTgg                                                                            <  1:4282081/93‑1 (MQ=255)

                  gTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCt                                                                        >  2:4835267/1‑93 (MQ=255)

                    cGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCttt                                                                      <  2:813657/93‑1 (MQ=255)

                    cGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCttt                                                                      <  1:4506191/93‑1 (MQ=255)

                    cGCAAAAACTGGCAGATGTACAGACACTGTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCttt                                                                      <  2:3213161/93‑1 (MQ=255)

                     gCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTtc                                                                     >  2:2229054/1‑93 (MQ=255)

                      cAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTtct                                                                    >  2:1193549/1‑93 (MQ=255)

                      cAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTtct                                                                    >  1:2488022/1‑93 (MQ=255)

                        aaaaCTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTcg                                                                  <  2:2093828/93‑1 (MQ=255)

                        aaaaCTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTcg                                                                  >  1:4836737/1‑93 (MQ=255)

                         aaaCTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTcgc                                                                 <  1:4058961/93‑1 (MQ=255)

                         aaaCTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTcgc                                                                 >  1:4067702/1‑93 (MQ=255)

                           aCTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCgcgc                                                               >  2:3838749/1‑93 (MQ=255)

                            cTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCgca                                                              <  1:3191761/93‑1 (MQ=255)

                               gCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCgcagca                                                           >  1:6146498/1‑93 (MQ=255)

                                 cGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACt                                                         <  2:1619626/93‑1 (MQ=255)

                                 cGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACt                                                         <  2:2869589/93‑1 (MQ=255)

                                  gATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTg                                                        <  2:437139/93‑1 (MQ=255)

                                  gATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTg                                                        >  2:4077681/1‑93 (MQ=255)

                                  gATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTg                                                        >  2:2676545/1‑93 (MQ=255)

                                   aTGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGc                                                       >  2:32579/1‑93 (MQ=255)

                                    tGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCg                                                      >  2:6285744/1‑93 (MQ=255)

                                    tGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCg                                                      <  2:3311057/93‑1 (MQ=255)

                                     gTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGt                                                     >  2:3939838/1‑93 (MQ=255)

                                      tACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTg                                                    <  1:5518142/93‑1 (MQ=255)

                                      tACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGCGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTg                                                    >  1:2835094/1‑93 (MQ=255)

                                      gaCAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTg                                                    <  1:5807615/92‑1 (MQ=255)

                                           acacTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCaga                                               <  2:5363066/93‑1 (MQ=255)

                                           acacTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCaga                                               <  1:3515198/93‑1 (MQ=255)

                                             acTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGATGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGGCAGAGt                                             >  2:2564057/1‑93 (MQ=255)

                                             acTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGGCAGAGt                                             >  2:2417268/1‑93 (MQ=255)

                                             acTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGt                                             <  2:35798/93‑1 (MQ=255)

                                             acTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGt                                             >  2:5919913/1‑93 (MQ=255)

                                              cTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTa                                            <  1:6245105/93‑1 (MQ=255)

                                              cTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTa                                            <  1:2238761/93‑1 (MQ=255)

                                              cTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTa                                            <  2:5540824/93‑1 (MQ=255)

                                               ttttGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTag                                           <  2:2970273/93‑1 (MQ=255)

                                               ttttGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTag                                           >  1:1012488/1‑93 (MQ=255)

                                               ttttGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTag                                           >  2:1081967/1‑93 (MQ=255)

                                                tttGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTaga                                          <  2:4091249/93‑1 (MQ=255)

                                                 ttGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTagaa                                         <  2:4006608/93‑1 (MQ=255)

                                                 ttGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTagaa                                         <  1:1246032/93‑1 (MQ=255)

                                                   gCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTagaaga                                       <  2:3359512/93‑1 (MQ=255)

                                                   gCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTagaaga                                       <  1:5964259/93‑1 (MQ=255)

                                                    cATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAagag                                      >  2:5969045/1‑93 (MQ=255)

                                                     aTGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGt                                     >  1:1580625/1‑93 (MQ=255)

                                                     aTGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGt                                     <  2:4321658/93‑1 (MQ=255)

                                                     aTGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGt                                     <  2:4613303/93‑1 (MQ=255)

                                                          gCGGGTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTgg                                >  2:5414709/1‑93 (MQ=255)

                                                            gggTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCa                              >  1:4132766/1‑93 (MQ=255)

                                                            gggTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCa                              >  1:4811557/1‑93 (MQ=255)

                                                             ggTGTGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCAt                             >  1:3448054/1‑93 (MQ=255)

                                                              gtgtGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCAtt                            <  2:2659009/93‑1 (MQ=255)

                                                              gtgtGGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCAtt                            <  2:579426/93‑1 (MQ=255)

                                                               tgtgGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCAttt                           <  1:771762/93‑1 (MQ=255)

                                                               tgtgGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCAttt                           >  1:1473145/1‑93 (MQ=255)

                                                                gtgGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTa                          >  2:4980793/1‑93 (MQ=255)

                                                                 tgGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTaa                         <  1:5630025/93‑1 (MQ=255)

                                                                 tgGCGGATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTaa                         >  2:5068852/1‑93 (MQ=255)

                                                                     ggATAACGTGATTTGTTTTGAGGCGCAAATGCCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGa                     <  1:613106/93‑1 (MQ=255)

                                                                     ggATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGa                     <  1:2091785/93‑1 (MQ=255)

                                                                     ggATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGa                     >  1:2949711/1‑93 (MQ=255)

                                                                      gATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGaa                    >  1:2595571/1‑93 (MQ=255)

                                                                      gATAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGaa                    >  1:991614/1‑93 (MQ=255)

                                                                       aTAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAAc                   <  1:2009736/93‑1 (MQ=255)

                                                                       aTAACGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAAc                   <  1:2862814/93‑1 (MQ=255)

                                                                           cGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACaaaa               <  1:108998/93‑1 (MQ=255)

                                                                           cGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACaaaa               <  1:5436482/93‑1 (MQ=255)

                                                                           cGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACaaaa               >  2:1403682/1‑93 (MQ=255)

                                                                           cGTGATTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACaaaa               >  1:3811459/1‑93 (MQ=255)

                                                                               aTTTGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGcc           >  2:2759333/1‑93 (MQ=255)

                                                                                tttGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCa          >  1:3047372/1‑93 (MQ=255)

                                                                                tttGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCa          <  2:1584970/93‑1 (MQ=255)

                                                                                tttGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCa          <  2:3041120/93‑1 (MQ=255)

                                                                                tttGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCa          >  2:3733995/1‑93 (MQ=255)

                                                                                 ttGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCAt         <  2:5649564/93‑1 (MQ=255)

                                                                                 ttGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCAt         <  2:5759683/93‑1 (MQ=255)

                                                                                 ttGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCAt         <  2:1295979/93‑1 (MQ=255)

                                                                                 ttGTTTTGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCAt         >  1:4417082/1‑93 (MQ=255)

                                                                                    ttttGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTa      <  1:5261592/93‑1 (MQ=255)

                                                                                     tttGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTAc     >  1:1680890/1‑93 (MQ=255)

                                                                                       tGAGGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGa   <  2:1689361/93‑1 (MQ=255)

                                                                                        gAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGaa  >  2:6244185/1‑93 (MQ=255)

CAGCCTGATGGCATCGCCGTCGCAAAAACTGGCCGATGTACAGACACTTTTGCATGAAGCGGGTGTGGCGGATAACGTGATTTGTTTTGAAGCGCAAATACCACTGGCGCTTTCTCGCGCAGCACTGCGTGCCAGAGTAGAAGAGTGCTGGCATTTAACTGAACAAAATGCCATGTACGAA  >  NC_000913/1461982‑1462162

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