BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,023,531	T→C	H318H (CAT→CAC)	yigC →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,023,531	0	T	C	99.5%	608.8 / ‑4.5	187	H318H (CAT→CAC)	yigC	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (116/70); total (117/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

Only 100 of 195 total aligned reads displayed.

TGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCATTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGT  >  NC_000913/4023441‑4023619

tGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTCTTTACCGTGACGCATATTACGCAGCGTGAAGATGCGATTTACCACTc                                                                                        >  1:2330380/1‑93 (MQ=255)

tGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTc                                                                                        >  2:2915021/1‑93 (MQ=255)

   cGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCAc                                                                                     <  1:3800842/93‑1 (MQ=255)

    gACCCCACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCAcc                                                                                    <  2:3781169/93‑1 (MQ=255)

    gACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCAcc                                                                                    >  2:1237438/1‑93 (MQ=255)

     aCCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCt                                                                                   >  1:3934479/1‑93 (MQ=255)

      ccACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCta                                                                                  >  2:2873488/1‑93 (MQ=255)

       cacaCCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCTACCtat                                                                                 >  1:446539/1‑93 (MQ=255)

         cacCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAc                                                                               <  1:2706203/93‑1 (MQ=255)

         cacCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAc                                                                               >  1:169298/1‑93 (MQ=255)

         cacCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAc                                                                               >  1:4612026/1‑93 (MQ=255)

            cGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCgg                                                                            <  2:310587/93‑1 (MQ=255)

              gTTACTATAATGAAGTCGATAGTTTCCCGGGATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGCTTTACCACTCCACCTATCCCaggg                                                                          >  1:4024588/1‑92 (MQ=255)

               ttACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGGGAAGATGCGATTTACCACTCCACCTATACCGGGCg                                                                         >  1:169746/1‑93 (MQ=255)

                  cTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTcc                                                                      >  2:1071543/1‑93 (MQ=255)

                     taATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGcc                                                                   >  2:3254480/1‑93 (MQ=255)

                        tGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGa                                                                <  2:2369519/93‑1 (MQ=255)

                        tGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGa                                                                >  2:3088155/1‑93 (MQ=255)

                        tGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCAACCTATACCGGGCGTCCGCCAGa                                                                >  2:2582277/1‑93 (MQ=255)

                        tGAAGTCGATAGTTTCCCGGTATTTACCATGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGa                                                                <  1:662442/93‑1 (MQ=255)

                         gAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGAt                                                               >  1:4550202/1‑93 (MQ=255)

                          aaGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATg                                                              <  2:3705735/93‑1 (MQ=255)

                           aGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGa                                                             >  2:442792/1‑93 (MQ=255)

                           aGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGa                                                             <  1:3320937/93‑1 (MQ=255)

                           aGTCGATAGTTTCCCGGTATTTACAGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGa                                                             >  1:1090643/1‑93 (MQ=255)

                            gTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAg                                                            >  2:2613376/1‑93 (MQ=255)

                             tCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGc                                                           >  1:1661280/1‑93 (MQ=255)

                              cGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGcc                                                          <  1:1311380/93‑1 (MQ=255)

                               gATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGccc                                                         >  2:2505263/1‑93 (MQ=255)

                               gATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGccc                                                         <  1:2033185/93‑1 (MQ=255)

                                 tAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcgc                                                       <  1:497313/93‑1 (MQ=255)

                                 tAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcgc                                                       >  1:1783396/1‑93 (MQ=255)

                                 tAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcgc                                                       >  1:2901679/1‑93 (MQ=255)

                                 tAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcgc                                                       >  2:2569377/1‑93 (MQ=255)

                                  aGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcgcg                                                      >  2:2002121/1‑93 (MQ=255)

                                  aGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcgcg                                                      >  1:3812463/1‑93 (MQ=255)

                                   gTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCgg                                                     <  2:72966/93‑1 (MQ=255)

                                     ttCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTg                                                   <  1:3882335/93‑1 (MQ=255)

                                     ttCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTg                                                   >  1:3593870/1‑93 (MQ=255)

                                      tCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGc                                                  <  2:3090277/93‑1 (MQ=255)

                                      tCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGAGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTgt                                                  >  2:539449/1‑92 (MQ=255)

                                        ccGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCCCTCCACCTATACCGGGCGTCCGCCAGATGAGCACGCGGTGCTg                                                <  1:3714883/93‑1 (MQ=255)

                                        ccGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTg                                                <  1:770250/93‑1 (MQ=255)

                                         cGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTgg                                               >  2:3324980/1‑93 (MQ=255)

                                         cGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGCGCTgg                                               >  2:4627936/1‑93 (MQ=255)

                                          ggTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTggg                                              >  1:2935587/1‑93 (MQ=255)

                                           gTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGCGCCTGCGGTGCTGGgt                                             >  1:552885/1‑93 (MQ=255)

                                           gTATTTACCGGGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATCCCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgt                                             <  2:1702030/93‑1 (MQ=255)

                                             aTTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgtgt                                           <  2:1098478/93‑1 (MQ=255)

                                             aTTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgtgt                                           <  1:417000/93‑1 (MQ=255)

                                              tttACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTc                                          >  1:3997645/1‑93 (MQ=255)

                                               ttACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCTAGATGAGCCCGCGGTGCTGGGTGTCg                                         <  2:2877239/93‑1 (MQ=255)

                                               ttACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCg                                         >  2:2659182/1‑93 (MQ=255)

                                                 aCCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCa                                       >  1:1663303/1‑93 (MQ=255)

                                                   cGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACt                                     <  2:3009840/93‑1 (MQ=255)

                                                   cGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACt                                     <  1:2979434/93‑1 (MQ=255)

                                                    gTGACGCATATTACCCAGCGTGAAGATGCGATTTACAACTCCACCTATACCGGGTGTCCGCCAGATGAGCCCGCGGGGCTGGGTGTCGCACta                                    >  1:3861693/1‑92 (MQ=255)

                                                     tGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGGACTGa                                   >  1:148273/1‑93 (MQ=255)

                                                     tGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGa                                   <  2:1096020/93‑1 (MQ=255)

                                                      gACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGaa                                  <  2:2086464/93‑1 (MQ=255)

                                                      gACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGaa                                  >  2:3593931/1‑93 (MQ=255)

                                                       aCGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAAc                                 <  1:850935/93‑1 (MQ=255)

                                                       aCGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAAc                                 <  1:4102265/93‑1 (MQ=255)

                                                         gCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGa                               >  2:1828855/1‑93 (MQ=255)

                                                           atatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCCCTGAACGAAg                             >  2:1556580/1‑93 (MQ=255)

                                                           atatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAg                             >  1:2781008/1‑93 (MQ=255)

                                                           atatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAg                             <  1:3414525/93‑1 (MQ=255)

                                                           atatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAg                             >  2:702879/1‑93 (MQ=255)

                                                            tatTACCCATCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgt                            >  1:1153780/1‑93 (MQ=255)

                                                            tatTACCCAGCGTGAAGATGCGATTTACCCCTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgt                            <  1:1180858/93‑1 (MQ=255)

                                                            tatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgt                            <  2:3537109/93‑1 (MQ=255)

                                                            tatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgt                            >  2:580008/1‑93 (MQ=255)

                                                            tatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTAGACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgt                            <  1:3479764/93‑1 (MQ=255)

                                                              ttACCCAGCGTGAAGATGCGATTTATCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgtgt                          >  1:1075113/1‑93 (MQ=255)

                                                              ttACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgtgt                          >  1:1075114/1‑93 (MQ=255)

                                                                tcccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGttt                        <  2:4060421/92‑1 (MQ=255)

                                                                aCCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGttt                        >  2:4689810/1‑93 (MQ=255)

                                                                aCCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGttt                        <  2:1864958/93‑1 (MQ=255)

                                                                    aGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGc                    >  1:9846/1‑93 (MQ=255)

                                                                      cGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCg                  <  1:1562511/93‑1 (MQ=255)

                                                                      cGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCg                  <  2:2392377/93‑1 (MQ=255)

                                                                      cGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCg                  >  2:3867766/1‑93 (MQ=255)

                                                                       gTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGa                 <  1:2059482/93‑1 (MQ=255)

                                                                       gTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGa                 >  2:4086532/1‑93 (MQ=255)

                                                                        tGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAt                >  2:3265178/1‑93 (MQ=255)

                                                                         gAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAtt               <  1:330252/93‑1 (MQ=255)

                                                                          aaGATGCGATTTACCACTCCACCTGTACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTc              >  2:2259818/1‑93 (MQ=255)

                                                                          aaGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTc              <  1:2601124/93‑1 (MQ=255)

                                                                           aGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGTCCGCGGGGCTGGGGGTCGCACTGAACGAAGTGTTTGTGCCGATCCt             >  1:4236251/1‑93 (MQ=255)

                                                                           aGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCt             >  2:2159659/1‑93 (MQ=255)

                                                                           aGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCt             <  1:2121603/93‑1 (MQ=255)

                                                                            gATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTg            <  1:432991/93‑1 (MQ=255)

                                                                             aTGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGc           <  1:2421095/93‑1 (MQ=255)

                                                                              tGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCa          >  2:417885/1‑93 (MQ=255)

                                                                              tGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCa          >  2:2775430/1‑93 (MQ=255)

                                                                                 gATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCaaaa       <  1:3247343/93‑1 (MQ=255)

                                                                                  aTTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCaaaaa      <  1:2571645/93‑1 (MQ=255)

                                                                                    ttACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACa    >  1:2467930/1‑93 (MQ=255)

                                                                                     tACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAg   >  2:2732995/1‑93 (MQ=255)

                                                                                      aCCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGt  >  1:1906326/1‑93 (MQ=255)

TGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCATTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGT  >  NC_000913/4023441‑4023619

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