BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,983,890	T→C	T88A (ACT→GCT)	araF ←	L‑arabinose‑binding periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,983,890	0	T	C	100.0%	374.7 / NA	120	T88A (ACT→GCT)	araF	L‑arabinose‑binding periplasmic protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (66/54); total (66/54)

Only 100 of 121 total aligned reads displayed.

TTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAA  >  NC_000913/1983802‑1983979

ttAGCAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaa                                                                                       <  2:2742541/93‑1 (MQ=255)

ttAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaa                                                                                       >  1:505314/1‑93 (MQ=255)

ttAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaa                                                                                       >  2:3544753/1‑93 (MQ=255)

ttAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaa                                                                                       >  2:2603810/1‑93 (MQ=255)

 tAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaaa                                                                                      <  1:2894970/93‑1 (MQ=255)

 tAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaaa                                                                                      <  2:379653/93‑1 (MQ=255)

  aaCAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAAt                                                                                     <  2:4559009/93‑1 (MQ=255)

   aCAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGCGCCGAGTTTGGGGTCCGGAGCACAAATa                                                                                    >  1:344054/1‑93 (MQ=255)

    cAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATaa                                                                                   <  2:4441743/93‑1 (MQ=255)

      aaCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACg                                                                                 >  1:1469306/1‑93 (MQ=255)

      aaCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACg                                                                                 >  2:267043/1‑93 (MQ=255)

       aCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGATTGGGGTCCGGAGCACAAATAACGa                                                                                >  2:1597563/1‑93 (MQ=255)

        cTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGGTTGGGGTCCGGAGCACAAATAACCaa                                                                               >  2:420994/1‑93 (MQ=255)

         tGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGaaa                                                                              >  2:4547776/1‑93 (MQ=255)

         tGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGaaa                                                                              >  1:4642058/1‑93 (MQ=255)

         tGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGaaa                                                                              <  1:2331426/93‑1 (MQ=255)

          ggTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAAc                                                                             <  1:3183620/93‑1 (MQ=255)

          ggTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAAc                                                                             >  1:2304393/1‑93 (MQ=255)

            tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           >  1:1803888/1‑93 (MQ=255)

            tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  2:4027277/93‑1 (MQ=255)

            tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           >  2:4214342/1‑93 (MQ=255)

            tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  2:1954929/93‑1 (MQ=255)

            tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           >  1:4626857/1‑93 (MQ=255)

            tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  2:134468/93‑1 (MQ=255)

              aTCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCttt                                                                         >  2:2798980/1‑93 (MQ=255)

              aTCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCttt                                                                         >  2:4064482/1‑93 (MQ=255)

              aTCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCttt                                                                         <  2:1279981/93‑1 (MQ=255)

              aTCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCttt                                                                         <  2:748235/93‑1 (MQ=255)

              aTCCACGGCAATGACTTTCATATCGTAGCCACGCGCTGTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCttt                                                                         <  2:454700/93‑1 (MQ=255)

               tCCCCGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCtttt                                                                        <  2:304167/93‑1 (MQ=255)

               tCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCtttt                                                                        >  2:2959306/1‑93 (MQ=255)

               tCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCtttt                                                                        >  2:1063515/1‑93 (MQ=255)

                ccACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTg                                                                       <  1:4661737/93‑1 (MQ=255)

                 cACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTgc                                                                      <  1:348933/93‑1 (MQ=255)

                    ggCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGcc                                                                   >  2:3532136/1‑93 (MQ=255)

                    ggCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGcc                                                                   <  1:1018492/93‑1 (MQ=255)

                     gCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCa                                                                  >  1:1067325/1‑93 (MQ=255)

                      cAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCAc                                                                 >  1:3593744/1‑93 (MQ=255)

                      cAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCAc                                                                 <  2:142730/93‑1 (MQ=255)

                      cAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCAc                                                                 <  1:2639094/93‑1 (MQ=255)

                        aTGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTg                                                               <  2:1865528/93‑1 (MQ=255)

                        aTGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTg                                                               >  2:813747/1‑93 (MQ=255)

                        aTGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTg                                                               <  1:2699198/93‑1 (MQ=255)

                           aCTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCa                                                            >  1:3258200/1‑93 (MQ=255)

                           aCTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCa                                                            >  2:2230178/1‑93 (MQ=255)

                             tttCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGc                                                          >  1:120602/1‑93 (MQ=255)

                             tttCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGc                                                          >  1:1676123/1‑93 (MQ=255)

                             tttCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGc                                                          >  2:3407017/1‑93 (MQ=255)

                               tCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCa                                                        >  2:789099/1‑93 (MQ=255)

                               tCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCa                                                        <  2:1810575/93‑1 (MQ=255)

                                cATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAg                                                       >  2:3805619/1‑93 (MQ=255)

                                cATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAg                                                       <  1:3074181/93‑1 (MQ=255)

                                  tatCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGc                                                     <  2:3271463/93‑1 (MQ=255)

                                   atCGTAGCCACGCGCTTTCGCGACGATGGTAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCt                                                    >  1:4525866/1‑93 (MQ=255)

                                   atCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCt                                                    >  1:3323896/1‑93 (MQ=255)

                                    tCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTg                                                   <  2:1441376/93‑1 (MQ=255)

                                      gTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTc                                                 >  2:3167518/1‑93 (MQ=255)

                                      gTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTc                                                 >  2:417524/1‑93 (MQ=255)

                                      gTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTc                                                 >  2:1179088/1‑93 (MQ=255)

                                      gTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTc                                                 >  1:108186/1‑93 (MQ=255)

                                         gCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGAt                                              <  1:245778/93‑1 (MQ=255)

                                          ccACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATc                                             <  2:321013/93‑1 (MQ=255)

                                          ccACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATc                                             <  1:1533386/93‑1 (MQ=255)

                                           cACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATcg                                            >  1:900189/1‑93 (MQ=255)

                                              gcgcTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGt                                         >  2:3381269/1‑93 (MQ=255)

                                              gcgcTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGt                                         <  2:3051367/93‑1 (MQ=255)

                                              gcgcTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGt                                         >  1:2049206/1‑93 (MQ=255)

                                                gcTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTc                                       >  1:1147316/1‑93 (MQ=255)

                                                 cTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGCTCGCGTTCa                                      >  1:2607716/1‑93 (MQ=255)

                                                 cTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCa                                      >  2:2391289/1‑93 (MQ=255)

                                                  tttCGCGACGATGGCAGAGCCGGGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCaa                                     >  2:1709155/1‑93 (MQ=255)

                                                  tttCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCaa                                     <  1:1623104/93‑1 (MQ=255)

                                                    tCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATg                                   >  1:2949793/1‑93 (MQ=255)

                                                      gcgACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGtt                                 >  1:4129915/1‑93 (MQ=255)

                                                        gacgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGtttt                               <  1:1187832/93‑1 (MQ=255)

                                                        gacgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGtttt                               >  1:1004534/1‑93 (MQ=255)

                                                          cgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGtttttt                             <  2:1057653/93‑1 (MQ=255)

                                                           gaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTc                            >  1:898983/1‑93 (MQ=255)

                                                           gaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTc                            >  1:2481862/1‑93 (MQ=255)

                                                            aTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAGATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTTAATGTTTTTTCg                           >  2:1543915/1‑93 (MQ=255)

                                                            aTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCg                           >  2:4658514/1‑93 (MQ=255)

                                                             tGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGc                          >  1:4344769/1‑93 (MQ=255)

                                                             tGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGc                          >  2:2158553/1‑93 (MQ=255)

                                                               gCAGAGCCGATTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCa                        <  2:1862582/93‑1 (MQ=255)

                                                               gCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCa                        <  2:1263770/93‑1 (MQ=255)

                                                                 agagCCGAGTTTGGGGTCCGGAGCNCAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATc                      <  2:704123/93‑1 (MQ=255)

                                                                  gagCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATcc                     >  1:3414120/1‑93 (MQ=255)

                                                                     ccGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggc                  >  1:1277566/1‑93 (MQ=255)

                                                                       gAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcac                <  1:1741629/93‑1 (MQ=255)

                                                                          tttGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacggc             >  1:4463737/1‑93 (MQ=255)

                                                                            tGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCaa           <  1:4258491/93‑1 (MQ=255)

                                                                             ggggTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCATGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAAt          <  1:1004532/93‑1 (MQ=255)

                                                                              gggTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATc         <  2:1317464/93‑1 (MQ=255)

                                                                                 tCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTta      <  2:2021277/93‑1 (MQ=255)

                                                                                 tCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGCTCAATGTTTTTTCGCCATCCGGCACGGCTATCTta      >  1:824117/1‑93 (MQ=255)

                                                                                  ccGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaa     >  2:1774961/1‑93 (MQ=255)

                                                                                   cGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaat    <  2:2465654/93‑1 (MQ=255)

                                                                                   cGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaat    <  1:4237372/93‑1 (MQ=255)

                                                                                     gAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaataa  <  2:3663659/93‑1 (MQ=255)

                                                                                     gAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaataa  <  2:2043785/93‑1 (MQ=255)

TTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAA  >  NC_000913/1983802‑1983979

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