BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,778,478	T→C	A18A (GCT→GCC)	ydiR →	putative flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,778,478	0	T	C	97.8%	421.6 / ‑0.5	136	A18A (GCT→GCC)	ydiR	putative flavoprotein
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (71/62); total (73/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 139 total aligned reads displayed.

AAGAAAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCTGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGCCTTATGG  >  NC_000913/1778388‑1778567

aagaaAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGa                                                                                         >  1:385081/1‑93 (MQ=255)

 agaaAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGaa                                                                                        >  2:1543752/1‑93 (MQ=255)

 agaaAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGaa                                                                                        >  2:3698772/1‑93 (MQ=255)

  gaaAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGaat                                                                                       >  1:3573261/1‑92 (MQ=255)

  gaaAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAAc                                                                                       <  1:332163/93‑1 (MQ=255)

   aaaGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACt                                                                                      >  1:3648127/1‑93 (MQ=255)

     aGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGt                                                                                    <  2:1076459/93‑1 (MQ=255)

       cccTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGttt                                                                                  >  1:4294951/1‑93 (MQ=255)

       cccTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGttt                                                                                  <  1:156431/93‑1 (MQ=255)

        ccTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTg                                                                                 >  1:4741809/1‑93 (MQ=255)

        ccTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTg                                                                                 <  2:1091754/93‑1 (MQ=255)

         ctgaactgaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTgg                                                                                >  1:4275713/1‑93 (MQ=255)

         ctgaactgaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTgg                                                                                >  1:3728250/1‑93 (MQ=255)

         ctgaactgaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTgg                                                                                <  1:2733193/93‑1 (MQ=255)

          tgaactgaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTggc                                                                               >  1:1989056/1‑93 (MQ=255)

          agagctgaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTggc                                                                               <  2:1762625/89‑1 (MQ=255)

           gaactgaaGCCTACGGACAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTggcg                                                                              >  2:2195831/1‑93 (MQ=255)

              ctgaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGCCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGgcggcg                                                                           >  2:3742646/1‑93 (MQ=255)

               tgaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGgcgc                                                                          >  1:1752186/1‑93 (MQ=255)

               tgaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGgcgc                                                                          >  2:2047719/1‑93 (MQ=255)

                gaaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCt                                                                         <  2:3284406/93‑1 (MQ=255)

                 aaGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTc                                                                        >  2:4509698/1‑93 (MQ=255)

                  aGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCa                                                                       <  2:519143/93‑1 (MQ=255)

                    ccTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGc                                                                     <  2:1461986/93‑1 (MQ=255)

                      tACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCaa                                                                   <  1:3019196/93‑1 (MQ=255)

                       aCGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAAt                                                                  <  2:3396026/93‑1 (MQ=255)

                       aCGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAAt                                                                  <  1:2150865/93‑1 (MQ=255)

                        cGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATg                                                                 <  1:2487506/93‑1 (MQ=255)

                         ggAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATgg                                                                <  1:2041211/93‑1 (MQ=255)

                          gagaAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATggg                                                               <  1:3330495/93‑1 (MQ=255)

                            gaagaaACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCTGAACTGTTTGGCGGCGCTCAGCAATGGGGc                                                             >  2:1993271/1‑93 (MQ=255)

                            gaagaaACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGc                                                             >  2:253626/1‑93 (MQ=255)

                            gaagaaACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGc                                                             >  2:3161250/1‑93 (MQ=255)

                              agaaACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCa                                                           <  1:4209147/93‑1 (MQ=255)

                              agaaACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCa                                                           >  2:4020351/1‑93 (MQ=255)

                               gaaACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCaa                                                          >  2:1714847/1‑93 (MQ=255)

                               gaaACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCaa                                                          >  2:517996/1‑93 (MQ=255)

                                aaaCGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAAc                                                         >  1:3952934/1‑93 (MQ=255)

                                 aaCGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACa                                                        <  2:3572172/93‑1 (MQ=255)

                                   cGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAgg                                                      >  1:1373544/1‑93 (MQ=255)

                                    gATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGgt                                                     <  1:1107523/93‑1 (MQ=255)

                                    gATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGgt                                                     <  1:1107522/93‑1 (MQ=255)

                                     aTGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGgtg                                                    <  2:2130424/93‑1 (MQ=255)

                                     aTGAGCCAATTAAACAGCGTCAGGGTCTTTAGCGATAATGCTGTACGGTATGCCGAACTGTTTGGCCGCGCTCAGCAACGGGGCCAGCAGgtg                                                    >  2:2275035/1‑93 (MQ=255)

                                      tGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGgtgt                                                   >  2:4082124/1‑93 (MQ=255)

                                      tGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGgtgt                                                   >  2:4082123/1‑93 (MQ=255)

                                      tGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGgtgt                                                   >  1:542492/1‑93 (MQ=255)

                                           cAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGcc                                              >  1:2044844/1‑93 (MQ=255)

                                           cAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGcc                                              >  1:1212983/1‑93 (MQ=255)

                                            aaTTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCa                                             >  2:1915627/1‑93 (MQ=255)

                                                taaCAGCGTCTTGGTCTTTAACGATTATCCTTAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGt                                         <  1:956372/92‑1 (MQ=255)

                                                aaaCAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGt                                         >  1:1686141/1‑93 (MQ=255)

                                                  aCAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTAc                                       >  2:4241456/1‑93 (MQ=255)

                                                   cAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACa                                      >  2:4632226/1‑93 (MQ=255)

                                                    aGCGGCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACaa                                     >  2:1249352/1‑93 (MQ=255)

                                                      cGTCTGGGTCTTTAGCGATAATACTGAACGTTATGCTGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACaaaa                                   >  1:4293374/1‑93 (MQ=255)

                                                       gTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAAt                                  >  1:1315221/1‑93 (MQ=255)

                                                        tCTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATa                                 <  2:3460539/93‑1 (MQ=255)

                                                         cTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATAc                                >  2:4648913/1‑93 (MQ=255)

                                                         cTGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATAc                                >  1:1582502/1‑93 (MQ=255)

                                                          tGGGTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATAcc                               <  2:590278/93‑1 (MQ=255)

                                                           gggTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCg                              <  1:2120414/93‑1 (MQ=255)

                                                             gTCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGAc                            <  1:4462524/93‑1 (MQ=255)

                                                              tCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGAcc                           <  2:45243/93‑1 (MQ=255)

                                                              tCTTTAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGGCTGCCATTGTACAAAATACAGAcc                           >  2:3154430/1‑93 (MQ=255)

                                                                tttAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAg                         >  2:3238667/1‑93 (MQ=255)

                                                                tttAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAg                         <  1:305021/93‑1 (MQ=255)

                                                                 ttAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAgg                        >  2:712036/1‑93 (MQ=255)

                                                                 ttAGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAgg                        <  1:2782658/93‑1 (MQ=255)

                                                                   aGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGgcg                      >  1:1766040/1‑93 (MQ=255)

                                                                   aGCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGgcg                      >  2:357543/1‑93 (MQ=255)

                                                                    gCGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGgcgc                     <  2:4227537/93‑1 (MQ=255)

                                                                     cGATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgca                    <  1:2203808/93‑1 (MQ=255)

                                                                      gATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgcag                   >  1:3473595/1‑93 (MQ=255)

                                                                      gATAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgcag                   <  2:1089458/93‑1 (MQ=255)

                                                                       aTAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgcagg                  <  1:3442719/93‑1 (MQ=255)

                                                                        tAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgcaggc                 >  2:1746613/1‑93 (MQ=255)

                                                                        tAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgcaggc                 >  1:2430373/1‑93 (MQ=255)

                                                                        tAATCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgcagga                 >  2:2212897/1‑92 (MQ=255)

                                                                         aacccTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgcaggca                <  1:1446608/90‑1 (MQ=255)

                                                                          aTCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCgcaggcag               >  2:2767872/1‑93 (MQ=255)

                                                                           tCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATTAGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGt              <  2:1607182/93‑1 (MQ=255)

                                                                           tCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGt              >  1:1968216/1‑93 (MQ=255)

                                                                           tCCTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGt              >  1:4227152/1‑93 (MQ=255)

                                                                             cTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTa            <  1:251564/93‑1 (MQ=255)

                                                                             cTGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTa            <  2:1872538/93‑1 (MQ=255)

                                                                              tGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTAt           <  1:1453175/93‑1 (MQ=255)

                                                                              tGAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTAt           >  1:2311034/1‑93 (MQ=255)

                                                                               gAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATg          <  1:2450378/93‑1 (MQ=255)

                                                                               gAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATg          <  1:3541813/93‑1 (MQ=255)

                                                                               gAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATg          >  1:1998624/1‑93 (MQ=255)

                                                                               gAACGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATg          <  1:3226913/93‑1 (MQ=255)

                                                                                aaCGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGc         <  1:4728976/93‑1 (MQ=255)

                                                                                 aCGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGcc        >  2:3345398/1‑93 (MQ=255)

                                                                                  cGTTATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGCCt       <  1:4402137/93‑1 (MQ=255)

                                                                                    ttATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGCCTTa     >  1:1133755/1‑93 (MQ=255)

                                                                                     tATGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGCCTTAt    >  2:4482155/1‑93 (MQ=255)

                                                                                      aTGCCGAACTGTTTGGCTGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATCCCGACCAGGCGCAGGCAGTTATGCCTTATg   <  1:2744757/93‑1 (MQ=255)

                                                                                      aTGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGCCTTATg   <  2:525678/93‑1 (MQ=255)

                                                                                       tGCCGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGCCTTATgg  >  2:2820352/1‑93 (MQ=255)

AAGAAAGCCCTGAACTGAAGCCTACGGAGAAGAAACGATGAGTCAATTAAACAGCGTCTGGGTCTTTAGCGATAATCCTGAACGTTATGCTGAACTGTTTGGCGGCGCTCAGCAATGGGGCCAACAGGTGTATGCCATTGTACAAAATACCGACCAGGCGCAGGCAGTTATGCCTTATGG  >  NC_000913/1778388‑1778567

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