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.6%	709.1 / ‑5.2	223	H318H (CAT→CAC)	yigC	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (136/86); minor base A (0/1); total (136/87)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 234 total aligned reads displayed.

TGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCATTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGTT  >  NC_000913/4023441‑4023620

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

tGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGGATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTc                                                                                         >  1:3566140/1‑93 (MQ=255)

 ggCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTcc                                                                                        <  1:1908089/93‑1 (MQ=255)

     aCCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGGGAAGATGCGATTTACCATTCCActt                                                                                    >  2:687895/1‑91 (MQ=255)

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

       cacaCCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCtat                                                                                  >  2:1310363/1‑93 (MQ=255)

       cacaCCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCtat                                                                                  <  1:4401926/93‑1 (MQ=255)

        acacCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCtata                                                                                 >  1:441706/1‑93 (MQ=255)

          acCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAcc                                                                               <  2:2236159/93‑1 (MQ=255)

          acCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAcc                                                                               >  1:2818188/1‑93 (MQ=255)

           ccGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAGCg                                                                              >  2:3426656/1‑93 (MQ=255)

            cGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCgg                                                                             >  2:2299681/1‑93 (MQ=255)

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

            cGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCgg                                                                             >  1:3078013/1‑93 (MQ=255)

              gTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGc                                                                           >  1:1601432/1‑93 (MQ=255)

              gTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGc                                                                           <  2:2152893/93‑1 (MQ=255)

               ttACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCg                                                                          <  2:2743542/93‑1 (MQ=255)

               ttACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCg                                                                          >  2:1791278/1‑93 (MQ=255)

                tACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGt                                                                         >  1:3133617/1‑93 (MQ=255)

                 aCTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTc                                                                        >  1:1254924/1‑93 (MQ=255)

                 aCTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTc                                                                        >  2:4258307/1‑93 (MQ=255)

                 aCTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTc                                                                        <  1:1432960/93‑1 (MQ=255)

                 aCTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTc                                                                        >  1:2765408/1‑93 (MQ=255)

                  cTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTcc                                                                       <  1:1518120/93‑1 (MQ=255)

                   tataATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGAGTCCg                                                                      >  2:1711841/1‑93 (MQ=255)

                    ataATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGc                                                                     >  2:3248671/1‑93 (MQ=255)

                         gAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGAt                                                                <  1:3612323/93‑1 (MQ=255)

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

                          aaGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATg                                                               >  1:3676589/1‑93 (MQ=255)

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

                            gTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAg                                                             <  1:3205202/93‑1 (MQ=255)

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

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

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

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

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

                                aTAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcg                                                         >  1:3599772/1‑93 (MQ=255)

                                aTAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcg                                                         <  1:2011359/93‑1 (MQ=255)

                                   gTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCgg                                                      >  2:1647393/1‑93 (MQ=255)

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

                                       cccGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCt                                                  >  1:2309011/1‑93 (MQ=255)

                                         cGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTgg                                                <  2:4465825/93‑1 (MQ=255)

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

                                           gTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgt                                              >  2:3400268/1‑93 (MQ=255)

                                           gTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgt                                              >  2:1959597/1‑93 (MQ=255)

                                           gTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgt                                              >  1:2148142/1‑93 (MQ=255)

                                            tATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgtg                                             >  1:2417331/1‑93 (MQ=255)

                                            tATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGCGCGTCCGCCAGATGAGCCCGCGGTGCTGGgtg                                             >  1:3749596/1‑93 (MQ=255)

                                              tttACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgtgtg                                           >  2:927192/1‑92 (MQ=255)

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

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

                                                     tGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGa                                    >  1:490162/1‑93 (MQ=255)

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

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

                                                       aCGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAAc                                  >  1:1582664/1‑93 (MQ=255)

                                                        cGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACg                                 <  1:909967/93‑1 (MQ=255)

                                                          cATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCCCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGaa                               >  1:2906853/1‑93 (MQ=255)

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

                                                             atTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgtg                            <  1:4636053/93‑1 (MQ=255)

                                                              ttACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgtgt                           <  1:47250/93‑1 (MQ=255)

                                                               tACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGtt                          <  1:3878395/93‑1 (MQ=255)

                                                               tACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGtt                          <  1:4258442/93‑1 (MQ=255)

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

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

                                                                 cccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtg                        >  2:4059792/1‑93 (MQ=255)

                                                                 cccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtg                        >  2:2580657/1‑93 (MQ=255)

                                                                 cccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtg                        <  1:2631246/93‑1 (MQ=255)

                                                                  ccAGCGTGAAGATGCGATTTACCACTCCTCCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgt                       <  1:2977127/93‑1 (MQ=255)

                                                                  ccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgt                       >  1:1812460/1‑93 (MQ=255)

                                                                  ccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgt                       >  2:4327968/1‑93 (MQ=255)

                                                                  ccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgt                       >  2:29227/1‑93 (MQ=255)

                                                                  ccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgt                       >  2:1210303/1‑93 (MQ=255)

                                                                  ccAGCGTGAAGATGCGATTTACAACTCCACCTATACCGGGCGTCCGCCAGATGAGCACGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgt                       >  1:746341/1‑93 (MQ=255)

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

                                                                     gCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGcc                    <  1:669233/93‑1 (MQ=255)

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

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

                                                                         gAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAtt                >  1:3216781/1‑93 (MQ=255)

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

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

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

                                                                          aaGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTc               >  1:4149973/1‑93 (MQ=255)

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

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

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

                                                                           aGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACAAAGTGTTTGTGCCGATTCt              <  2:4562033/93‑1 (MQ=255)

                                                                             aTGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGc            >  2:4367420/1‑93 (MQ=255)

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

                                                                              tGCGATTTACCACGCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCa           <  2:3063301/93‑1 (MQ=255)

                                                                                 gATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCaaaa        >  2:1432864/1‑93 (MQ=255)

                                                                                 gATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGGGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCaaaa        >  2:559825/1‑93 (MQ=255)

                                                                                  aTTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCaaaaa       >  1:4297856/1‑93 (MQ=255)

                                                                                   tttACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAAc      <  1:1829690/93‑1 (MQ=255)

                                                                                   tttACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAAc      <  1:546112/93‑1 (MQ=255)

                                                                                    ttACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACa     <  2:2636519/93‑1 (MQ=255)

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

                                                                                     tgccACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAg    <  1:4473924/91‑1 (MQ=255)

                                                                                     tACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAg    <  1:3999810/93‑1 (MQ=255)

                                                                                      aCCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGt   <  2:4146176/93‑1 (MQ=255)

                                                                                       ccACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGtt  >  1:1238201/1‑93 (MQ=255)

TGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCATTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGTT  >  NC_000913/4023441‑4023620

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