BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,587	T→C	S116G (AGT→GGT)	ydeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,587	0	T	C	99.5%	667.7 / NA	215	S116G (AGT→GGT)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (119/95); minor base G (1/0); total (120/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

GCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGAT  >  NC_000913/1616497‑1616674

gCACTGTCGACGCCGAGCTAATAAGCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCgg                                                                                       <  1:3565164/93‑1 (MQ=255)

 cACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCggg                                                                                      >  1:5433882/1‑93 (MQ=255)

   cTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGtt                                                                                    >  1:1101949/1‑93 (MQ=255)

    tgtgGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTc                                                                                   >  2:4463216/1‑93 (MQ=255)

       ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                >  2:1539646/1‑93 (MQ=255)

       ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                >  2:6254143/1‑93 (MQ=255)

       ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                <  2:1761666/93‑1 (MQ=255)

       ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                >  1:967702/1‑93 (MQ=255)

       ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                <  1:5322401/93‑1 (MQ=255)

       ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                <  1:2284671/93‑1 (MQ=255)

        gACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCt                                                                               <  1:7528028/93‑1 (MQ=255)

        gACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCt                                                                               >  2:3440746/1‑93 (MQ=255)

         aCGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCtt                                                                              <  1:6386617/93‑1 (MQ=255)

           gCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTcc                                                                            <  2:207606/93‑1 (MQ=255)

           gCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTcc                                                                            <  2:3775129/93‑1 (MQ=255)

              gAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGc                                                                         <  2:759095/93‑1 (MQ=255)

              gAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGc                                                                         >  2:6122165/1‑93 (MQ=255)

                gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctc                                                                       >  2:2511725/1‑93 (MQ=255)

                gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctc                                                                       <  1:1362839/93‑1 (MQ=255)

                gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctc                                                                       <  1:85872/93‑1 (MQ=255)

                gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctc                                                                       >  2:5780994/1‑93 (MQ=255)

                gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTACAGctc                                                                       >  2:1388421/1‑93 (MQ=255)

                gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGATCATCTTCCAGctc                                                                       >  1:1094532/1‑93 (MQ=255)

                 cTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctct                                                                      >  2:6909461/1‑93 (MQ=255)

                  taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCtt                                                                     >  2:6783375/1‑93 (MQ=255)

                  taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCtt                                                                     <  2:3230820/93‑1 (MQ=255)

                  taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCtt                                                                     >  1:1137886/1‑93 (MQ=255)

                   aataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTc                                                                    >  1:4252289/1‑93 (MQ=255)

                   aataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTc                                                                    <  1:7160992/93‑1 (MQ=255)

                    ataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTcc                                                                   >  1:460481/1‑93 (MQ=255)

                     taatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCg                                                                  >  2:7612523/1‑93 (MQ=255)

                     taatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCg                                                                  >  1:340099/1‑93 (MQ=255)

                      aatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGa                                                                 >  2:2735946/1‑93 (MQ=255)

                       atcatcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGAc                                                                >  1:5462831/1‑93 (MQ=255)

                        tcatcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACt                                                               >  2:5579827/1‑93 (MQ=255)

                        tcatcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGATCTTCCGACt                                                               >  1:4143105/1‑93 (MQ=255)

                         catcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACtt                                                              >  2:5165302/1‑93 (MQ=255)

                          atcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTg                                                             >  1:6264195/1‑93 (MQ=255)

                          atcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTg                                                             >  2:3069861/1‑93 (MQ=255)

                           tcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGc                                                            >  1:3892016/1‑93 (MQ=255)

                           tcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGc                                                            <  1:3664129/93‑1 (MQ=255)

                           tcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGc                                                            >  2:2689227/1‑93 (MQ=255)

                             gCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCtt                                                          <  2:3324609/93‑1 (MQ=255)

                              cTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCttt                                                         <  2:1690319/93‑1 (MQ=255)

                                aTGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTc                                                       <  2:5186617/93‑1 (MQ=255)

                                aTGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTc                                                       >  1:4790447/1‑93 (MQ=255)

                                   gTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGcc                                                    >  2:5746617/1‑93 (MQ=255)

                                         ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                              >  1:6795905/1‑93 (MQ=255)

                                         ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                              <  2:7146127/93‑1 (MQ=255)

                                         ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                              >  1:5897247/1‑93 (MQ=255)

                                           aCCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGa                                            <  1:1443383/93‑1 (MQ=255)

                                            ccGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaa                                           >  2:1241168/1‑93 (MQ=255)

                                            ccGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaa                                           <  2:5854408/93‑1 (MQ=255)

                                            ccGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaa                                           <  2:3529064/93‑1 (MQ=255)

                                            ccGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaa                                           >  1:413780/1‑93 (MQ=255)

                                              ggCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatc                                         >  1:7377226/1‑93 (MQ=255)

                                              ggCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatc                                         <  2:2414599/93‑1 (MQ=255)

                                               gCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatca                                        >  2:2134103/1‑93 (MQ=255)

                                                 aGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaat                                      >  2:4188678/1‑93 (MQ=255)

                                                 aGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaat                                      <  2:7382661/93‑1 (MQ=255)

                                                      ttGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCtt                                 >  1:1829184/1‑93 (MQ=255)

                                                      ttGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCtt                                 >  2:3622924/1‑93 (MQ=255)

                                                      ttGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCtt                                 >  2:3742423/1‑93 (MQ=255)

                                                      ttGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCAGACTTGATTTTCGCCTCCGGTGAATCAATCGCtt                                 >  2:2660044/1‑93 (MQ=255)

                                                       tGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCttt                                <  2:6304770/93‑1 (MQ=255)

                                                       tGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCttt                                >  2:6597037/1‑93 (MQ=255)

                                                        ggCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTc                               <  1:7762090/93‑1 (MQ=255)

                                                        ggCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTc                               <  2:7525050/93‑1 (MQ=255)

                                                         gCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCt                              >  1:4298668/1‑93 (MQ=255)

                                                         gCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCt                              <  2:2577052/93‑1 (MQ=255)

                                                         gCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGGTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGctctct                              >  2:4235912/1‑93 (MQ=255)

                                                            ttGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGtt                           >  1:247839/1‑93 (MQ=255)

                                                            ttGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGtt                           >  2:5814723/1‑93 (MQ=255)

                                                              gCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGc                         >  2:552738/1‑93 (MQ=255)

                                                               ccAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCg                        <  1:1575401/93‑1 (MQ=255)

                                                                cAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCgg                       >  1:1993003/1‑93 (MQ=255)

                                                                cAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCgg                       >  1:6302297/1‑93 (MQ=255)

                                                                  gCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGaa                     >  2:7694603/1‑93 (MQ=255)

                                                                    ggCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAag                   <  1:4968668/93‑1 (MQ=255)

                                                                     gCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAaga                  <  1:2115950/93‑1 (MQ=255)

                                                                       acaAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGACTCAATCGCTTTCTGTTGCGGAAAGAGc                >  2:5685345/1‑93 (MQ=255)

                                                                       acaAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGc                >  1:498262/1‑93 (MQ=255)

                                                                        caAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCACTCGCTTTCTGTTGCGGCAAGAGCa               >  2:2721319/1‑93 (MQ=255)

                                                                        caAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGct               >  2:5540653/1‑92 (MQ=255)

                                                                        caAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCa               <  2:1146991/93‑1 (MQ=255)

                                                                          aaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATc             <  1:2586224/93‑1 (MQ=255)

                                                                          aaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATc             <  1:3140600/93‑1 (MQ=255)

                                                                          aaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATc             >  1:7293681/1‑93 (MQ=255)

                                                                           aGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATcc            <  1:6612845/93‑1 (MQ=255)

                                                                           aGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATcc            >  1:5992443/1‑93 (MQ=255)

                                                                            ggCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTTCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCg           <  1:3790909/93‑1 (MQ=255)

                                                                            ggCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCg           >  2:5234110/1‑93 (MQ=255)

                                                                            ggCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCg           >  1:126041/1‑93 (MQ=255)

                                                                               aaTATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaa        >  1:4773928/1‑93 (MQ=255)

                                                                               aaTATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaa        >  2:50693/1‑93 (MQ=255)

                                                                                 tatTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAAc      >  2:4132743/1‑93 (MQ=255)

                                                                                   ttGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCg    <  2:2058344/93‑1 (MQ=255)

                                                                                   ttGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCg    <  1:3097565/93‑1 (MQ=255)

                                                                                   ttGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCg    >  2:3982976/1‑93 (MQ=255)

                                                                                     ggCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGAt  <  1:931869/93‑1 (MQ=255)

GCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGAT  >  NC_000913/1616497‑1616674

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