BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	865,377	C→T	G71R (GGG→AGG)	moeA ←	MoeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	865,377	0	C	T	99.6%	921.5 / ‑6.1	254	G71R (GGG→AGG)	moeA	MoeA
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (115/137); total (116/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 262 total aligned reads displayed.

CGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACATCGCTC  >  NC_000913/865287‑865465

cGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTgg                                                                                        >  2:2381403/1‑93 (MQ=255)

cGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGGTGCCTgg                                                                                        >  1:784754/1‑93 (MQ=255)

 ggTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGa                                                                                       <  2:9734775/93‑1 (MQ=255)

  gTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAg                                                                                      <  2:6122543/93‑1 (MQ=255)

     ataataCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCa                                                                                   <  2:7300603/93‑1 (MQ=255)

      taataCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCaa                                                                                  >  2:4372661/1‑93 (MQ=255)

        ataCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAata                                                                                <  2:3122048/93‑1 (MQ=255)

         taCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAatat                                                                               <  1:6197220/93‑1 (MQ=255)

         taCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAatat                                                                               <  2:10026731/93‑1 (MQ=255)

         taCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAatat                                                                               <  1:8850829/93‑1 (MQ=255)

         taCGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAatat                                                                               <  2:5668338/93‑1 (MQ=255)

             aaTGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGTCAATATCGGc                                                                           >  2:5370097/1‑93 (MQ=255)

             aaTGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATAGGc                                                                           <  2:8459256/93‑1 (MQ=255)

              aTGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCt                                                                          >  1:9119983/1‑93 (MQ=255)

                gCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTaa                                                                        <  1:6314587/93‑1 (MQ=255)

                gCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTaa                                                                        >  1:3391907/1‑93 (MQ=255)

                 cAGGTGCCCGCAGGGCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTaaa                                                                       <  1:2754186/93‑1 (MQ=255)

                 cAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTaaa                                                                       <  2:156504/93‑1 (MQ=255)

                   ggTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTCAACg                                                                     >  2:8854814/1‑93 (MQ=255)

                   ggTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACg                                                                     >  2:3157585/1‑93 (MQ=255)

                   ggTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACg                                                                     >  2:8053519/1‑93 (MQ=255)

                   ggTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACg                                                                     >  1:1013399/1‑93 (MQ=255)

                      aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                  <  2:2743783/93‑1 (MQ=255)

                      aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                  <  2:1525791/93‑1 (MQ=255)

                      aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                  <  2:7065649/93‑1 (MQ=255)

                      aCCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAc                                                                  <  1:6816617/93‑1 (MQ=255)

                       cccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAcc                                                                 >  1:9699815/1‑93 (MQ=255)

                       cccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCAcc                                                                 >  1:8323624/1‑93 (MQ=255)

                        ccGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcg                                                                <  2:3455732/93‑1 (MQ=255)

                          gCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcgcg                                                              >  1:3186592/1‑93 (MQ=255)

                          gCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcgcg                                                              >  2:2638599/1‑93 (MQ=255)

                          gCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcgcg                                                              <  1:3089791/93‑1 (MQ=255)

                          gCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACcgcg                                                              >  2:3278074/1‑93 (MQ=255)

                           cAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGt                                                             >  2:6637860/1‑93 (MQ=255)

                           cAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGt                                                             >  2:6478714/1‑93 (MQ=255)

                           cAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGt                                                             <  2:4543112/93‑1 (MQ=255)

                            aGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTa                                                            <  1:2649583/93‑1 (MQ=255)

                            aGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTa                                                            >  1:1305325/1‑93 (MQ=255)

                             ggCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAg                                                           >  1:3142632/1‑93 (MQ=255)

                              gCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGc                                                          <  2:9582171/93‑1 (MQ=255)

                              gCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGc                                                          <  1:9311554/93‑1 (MQ=255)

                               ccATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGcc                                                         <  1:992074/93‑1 (MQ=255)

                                cATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCg                                                        <  1:7497841/93‑1 (MQ=255)

                                cATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCg                                                        <  1:10026656/93‑1 (MQ=255)

                                cATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCg                                                        >  1:5433706/1‑93 (MQ=255)

                                 aTTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGt                                                       >  2:4460280/1‑93 (MQ=255)

                                  ttCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTc                                                      >  1:1664825/1‑93 (MQ=255)

                                  ttCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTc                                                      <  2:5098526/93‑1 (MQ=255)

                                   tCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTcc                                                     <  1:4574604/93‑1 (MQ=255)

                                   tCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTcc                                                     >  2:5996906/1‑93 (MQ=255)

                                    cACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCa                                                    <  2:4828453/93‑1 (MQ=255)

                                     cccATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCAt                                                   <  1:4703016/92‑1 (MQ=255)

                                     aCCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCAt                                                   <  1:6707303/93‑1 (MQ=255)

                                       cATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTg                                                 <  2:4588611/93‑1 (MQ=255)

                                         tGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCg                                               >  1:4045538/1‑93 (MQ=255)

                                              tGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGCCCATTGCGGAGtt                                          >  1:8953131/1‑93 (MQ=255)

                                               ggCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTa                                         <  2:3392081/93‑1 (MQ=255)

                                               ggCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTa                                         <  1:3021490/93‑1 (MQ=255)

                                                 cTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATc                                       <  1:2169242/93‑1 (MQ=255)

                                                 cTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATc                                       <  2:7162390/93‑1 (MQ=255)

                                                 cTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATc                                       <  2:5320105/93‑1 (MQ=255)

                                                   gACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCaa                                     >  2:2159529/1‑93 (MQ=255)

                                                   gACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCaa                                     >  2:9632271/1‑93 (MQ=255)

                                                   gACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCaa                                     <  1:4134302/93‑1 (MQ=255)

                                                    aCCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCaaa                                    <  1:351581/93‑1 (MQ=255)

                                                      cGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcc                                  <  2:1724230/93‑1 (MQ=255)

                                                       ggCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAccc                                 >  1:4898172/1‑93 (MQ=255)

                                                        gCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAAcccc                                <  1:1728294/93‑1 (MQ=255)

                                                          aaaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCgg                              >  1:2352519/1‑93 (MQ=255)

                                                           aaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGa                             >  1:1896895/1‑93 (MQ=255)

                                                           aaGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGa                             >  2:9889609/1‑93 (MQ=255)

                                                              gATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACa                          >  1:5963504/1‑93 (MQ=255)

                                                              gATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACa                          <  2:695138/93‑1 (MQ=255)

                                                               aTTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACAt                         <  1:7263438/93‑1 (MQ=255)

                                                               aTTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACAt                         <  2:123811/93‑1 (MQ=255)

                                                                tttACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATc                        >  2:6772788/1‑93 (MQ=255)

                                                                tttACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATc                        <  1:4321684/93‑1 (MQ=255)

                                                                 ttACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCa                       >  1:4370925/1‑93 (MQ=255)

                                                                 ttACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCa                       <  2:1839006/93‑1 (MQ=255)

                                                                 ttACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCa                       >  1:8213947/1‑93 (MQ=255)

                                                                   aCCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAg                     <  2:3014662/93‑1 (MQ=255)

                                                                     cGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcg                   <  2:6314962/93‑1 (MQ=255)

                                                                     cGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcg                   <  1:3647984/93‑1 (MQ=255)

                                                                       gCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAgcggc                 <  1:8157413/93‑1 (MQ=255)

                                                                         aaCGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGa               >  1:139509/1‑93 (MQ=255)

                                                                         aaCGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGa               >  2:6620790/1‑93 (MQ=255)

                                                                          aCGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGaa              >  2:7832596/1‑93 (MQ=255)

                                                                          aCGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGaa              >  1:4757832/1‑93 (MQ=255)

                                                                           cGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGaaa             <  1:1740144/93‑1 (MQ=255)

                                                                              gCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGa          <  1:4890709/93‑1 (MQ=255)

                                                                              gCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGa          >  2:3774222/1‑93 (MQ=255)

                                                                                aGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACa        >  2:9186083/1‑93 (MQ=255)

                                                                                 gCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAt       <  2:8229612/93‑1 (MQ=255)

                                                                                 gCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAt       <  1:3566669/93‑1 (MQ=255)

                                                                                  cGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtc      <  2:2099656/93‑1 (MQ=255)

                                                                                   ggTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcg     <  2:7796093/93‑1 (MQ=255)

                                                                                    gTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAGCCCCGGAACCTCAAGCGGCGAAACGACAacgc    <  2:5226629/93‑1 (MQ=255)

                                                                                    gTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcgc    >  2:3947945/1‑93 (MQ=255)

                                                                                     ttGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcgct   <  1:1544599/93‑1 (MQ=255)

                                                                                     ttGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACAtcgct   <  2:8604737/93‑1 (MQ=255)

CGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATCAAACCCCGGAACATCAAGCGGCGAAACGACATCGCTC  >  NC_000913/865287‑865465

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: `‑`