BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,048,676	A→G	S104P (TCC→CCC)	ymcD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,048,676	0	A	G	97.4%	472.8 / 1.5	155	S104P (TCC→CCC)	ymcD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (4/0); new base G (85/66); total (89/66)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

Only 100 of 163 total aligned reads displayed.

AATGCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGACGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGCCA  >  NC_000913/1048584‑1048765

  tGCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCg                                                                                         >  2:4594153/1‑93 (MQ=255)

  tGCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTCTAATTTTGACGCGTTAAAGCAGTATTAGCGTGTAGTGGACg                                                                                         >  1:3801825/1‑93 (MQ=255)

   gCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgt                                                                                        >  2:2556798/1‑93 (MQ=255)

   gCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGACgt                                                                                        >  1:4273782/1‑93 (MQ=255)

    cAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtg                                                                                       <  2:4290367/93‑1 (MQ=255)

       aTCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCgtggtg                                                                                    >  2:1473012/1‑93 (MQ=255)

        tctcTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGtggtgg                                                                                   >  2:3261999/1‑93 (MQ=255)

         ctctAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTggtggt                                                                                  <  2:265072/93‑1 (MQ=255)

            tAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGa                                                                               >  2:660160/1‑93 (MQ=255)

             aaaGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAg                                                                              <  1:2257994/93‑1 (MQ=255)

              aaGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTATTGGGCGTGGTGGTTGAgg                                                                             >  1:3534415/1‑93 (MQ=255)

              aaGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAgg                                                                             <  2:1721571/93‑1 (MQ=255)

                gCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGgtg                                                                           <  1:3756217/93‑1 (MQ=255)

                gCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGgtg                                                                           >  2:2503144/1‑93 (MQ=255)

                gCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGgtg                                                                           >  1:2268772/1‑93 (MQ=255)

                gCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGgtg                                                                           >  2:828421/1‑93 (MQ=255)

                  ccGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGtt                                                                         >  1:2673781/1‑93 (MQ=255)

                     gTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATg                                                                      >  1:3519502/1‑93 (MQ=255)

                      tACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATga                                                                     <  1:4727270/93‑1 (MQ=255)

                      tACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATga                                                                     <  1:2583044/93‑1 (MQ=255)

                        cGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATgaga                                                                   <  1:2896524/93‑1 (MQ=255)

                        cGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATgaga                                                                   >  2:3048792/1‑93 (MQ=255)

                         gAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGaa                                                                  >  1:1578246/1‑93 (MQ=255)

                         gAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGaa                                                                  <  1:3620432/93‑1 (MQ=255)

                         gAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGaa                                                                  >  1:2936828/1‑93 (MQ=255)

                          aaTATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAAc                                                                 <  2:269594/93‑1 (MQ=255)

                          aaTATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAAc                                                                 <  2:3614822/93‑1 (MQ=255)

                           atatCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAAcc                                                                >  2:3746211/1‑93 (MQ=255)

                           atatCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTAGGCGTGGTGGTTGAGGTGTTATGAGAAcc                                                                >  1:2243177/1‑93 (MQ=255)

                            tatCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGCAGTTTTGGCAACCg                                                               >  2:807596/1‑93 (MQ=255)

                             atCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGc                                                              <  1:269481/93‑1 (MQ=255)

                                  aaCCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCg                                                         >  1:2184940/1‑93 (MQ=255)

                                  aaCCGGGGATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCg                                                         >  1:1993053/1‑93 (MQ=255)

                                   aCCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGCGGTGTTATGAGAACCGCCATCGc                                                        >  2:2594362/1‑93 (MQ=255)

                                   aCCGGGTATCATGTAACTTTAATTTCGACGCGTTAAAGCAGTATAAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCTCCATCGc                                                        <  2:1604705/93‑1 (MQ=255)

                                     cGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCa                                                      >  1:1085548/1‑93 (MQ=255)

                                     cGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCa                                                      <  1:825371/93‑1 (MQ=255)

                                     cGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCa                                                      <  1:2306471/93‑1 (MQ=255)

                                      gggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCAc                                                     >  1:2416356/1‑93 (MQ=255)

                                      gggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCAc                                                     <  1:1200336/93‑1 (MQ=255)

                                      gggTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCAc                                                     <  1:463497/93‑1 (MQ=255)

                                        gTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTg                                                   >  2:2036434/1‑93 (MQ=255)

                                        gTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTg                                                   >  2:1396608/1‑93 (MQ=255)

                                         tATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTTTGAGAACCGCCATCGCCACTGt                                                  <  1:64251/93‑1 (MQ=255)

                                         tATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGt                                                  <  1:4576045/93‑1 (MQ=255)

                                         tATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGt                                                  >  2:887398/1‑93 (MQ=255)

                                          aTCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTTGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTc                                                 <  1:4190139/93‑1 (MQ=255)

                                          aTCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTc                                                 >  2:1521934/1‑93 (MQ=255)

                                          aTCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTc                                                 >  1:96052/1‑93 (MQ=255)

                                            cATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAt                                               <  2:3482642/93‑1 (MQ=255)

                                            cATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAt                                               <  2:3482643/93‑1 (MQ=255)

                                             aTTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAtt                                              >  1:4418466/1‑93 (MQ=255)

                                              tttAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAttt                                             <  2:658488/93‑1 (MQ=255)

                                              tttAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCAttt                                             >  1:4075107/1‑93 (MQ=255)

                                               ttAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTtg                                            >  2:3847925/1‑93 (MQ=255)

                                               ttAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTtg                                            <  2:273471/93‑1 (MQ=255)

                                               ttAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTtg                                            >  2:1702940/1‑93 (MQ=255)

                                                tAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTtgc                                           >  2:1222942/1‑93 (MQ=255)

                                                 aaCTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTtgct                                          <  2:2556876/93‑1 (MQ=255)

                                                  aCTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTtgctg                                         >  1:2604986/1‑93 (MQ=255)

                                                    tttAATTTCGACGCGTTAAAGCAGTATTAGCGTGTGGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGca                                       >  1:3970310/1‑93 (MQ=255)

                                                     ttcaTTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGcaa                                      >  2:4625020/4‑93 (MQ=255)

                                                      tAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGcaac                                     <  1:4286036/93‑1 (MQ=255)

                                                        aTTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGcaacaa                                   >  2:2750897/1‑93 (MQ=255)

                                                         tttCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCaacaac                                  <  1:3149446/93‑1 (MQ=255)

                                                          ttCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACt                                 <  2:3943938/93‑1 (MQ=255)

                                                          ttCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACt                                 >  2:2176839/1‑93 (MQ=255)

                                                          ttCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACt                                 <  1:3776941/93‑1 (MQ=255)

                                                           tCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTg                                <  2:1894482/93‑1 (MQ=255)

                                                            cGACGCGTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGc                               >  1:1128705/1‑93 (MQ=255)

                                                               cgcgTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCaa                            <  2:741877/93‑1 (MQ=255)

                                                                gcgTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAAc                           <  2:3731543/93‑1 (MQ=255)

                                                                gcgTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAAc                           <  2:2087552/93‑1 (MQ=255)

                                                                  gTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGt                         <  1:3301865/93‑1 (MQ=255)

                                                                  gTTAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGt                         >  2:3003145/1‑93 (MQ=255)

                                                                   ttAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTc                        >  2:3895801/1‑93 (MQ=255)

                                                                   ttAAAGCAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTc                        <  1:3142176/93‑1 (MQ=255)

                                                                      aaGCAGTATTAGCGTGTAGTGGACGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGcc                     >  2:2356376/1‑93 (MQ=255)

                                                                         cAGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCcaca                  >  1:2089114/1‑93 (MQ=255)

                                                                          aGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCacac                 >  1:2881538/1‑93 (MQ=255)

                                                                          aGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCacac                 >  1:4547786/1‑93 (MQ=255)

                                                                          aGTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCacac                 <  2:568078/93‑1 (MQ=255)

                                                                          aGTATTAGCGTGTAGTGGACGTGGGGCTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCacac                 >  1:1336556/1‑93 (MQ=255)

                                                                           gTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACAcc                >  1:3582787/1‑93 (MQ=255)

                                                                           gTATTAGCGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACAcc                >  1:157322/1‑93 (MQ=255)

                                                                              ttAGCGTGGAGGGGGCGTGGTGGTTGAGGTGTTAGGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCACCACACCAAc             >  2:1277506/1‑93 (MQ=255)

                                                                                  ggtgtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGcc         <  1:1397884/92‑1 (MQ=255)

                                                                                  cGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGcc         >  1:391483/1‑93 (MQ=255)

                                                                                  cGTGTAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGcc         <  1:1782702/93‑1 (MQ=255)

                                                                                   gtgtAGTGGGCGTGGTGGTTGAGGTGTTATGAGATCCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCg        <  1:346174/93‑1 (MQ=255)

                                                                                   gtgtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCg        <  2:2351892/93‑1 (MQ=255)

                                                                                   gtgtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCg        >  1:4059241/1‑93 (MQ=255)

                                                                                   gtgtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCg        >  1:727179/1‑93 (MQ=255)

                                                                                    tgtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCTACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGt       >  2:2442575/1‑93 (MQ=255)

                                                                                     gtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGtt      >  2:4410835/1‑93 (MQ=255)

                                                                                     gtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGtt      <  2:3692628/93‑1 (MQ=255)

                                                                                     gtAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGtt      <  2:1851783/93‑1 (MQ=255)

                                                                                      tAGTGGGCGTGTTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTg     <  2:4157038/93‑1 (MQ=255)

                                                                                      tAGTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTg     >  1:2683481/1‑93 (MQ=255)

                                                                                        gTGGGCGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGcc   >  2:444327/1‑93 (MQ=255)

AATGCAAATCTCTAAAGCCCGGTACGAATATCGCAACCGGGTATCATTTAACTTTAATTTCGACGCGTTAAAGCAGTATTAGCGTGTAGTGGACGTGGTGGTTGAGGTGTTATGAGAACCGCCATCGCCACTGTCATTTGCTGCAACAACTGCCAACGTCGCCACACCAACAGCCGTTGCCA  >  NC_000913/1048584‑1048765

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