BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,926,891	G→C	19.7%	S650* (TCA→TGA)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,926,891	0	G	C	19.7%	397.6 / 92.3	163	S650* (TCA→TGA)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base G (66/65); new base C (14/18); total (80/83)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 165 total aligned reads displayed.

GTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAAT  >  NC_000913/1926792‑1926990

gTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTg                                                                                                     <  2:4331863/100‑1 (MQ=255)

   gCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                                                  <  1:3310426/100‑1 (MQ=255)

    cGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGTCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                                                                 <  1:2443840/100‑1 (MQ=255)

    cGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCTACCAGATTTACCGCCATGGCCTGAGt                                                                                                  <  2:1125471/99‑1 (MQ=255)

    cGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                                                                 <  1:223342/100‑1 (MQ=255)

     gggTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGCCCTGAGTcc                                                                                                >  1:4164201/1‑100 (MQ=255)

      ggTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa                                                                                               <  1:3924796/100‑1 (MQ=255)

        tAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATg                                                                                             >  1:785706/1‑100 (MQ=255)

        gtaTGTTCCTCTCGGCATCGCGACCATAATAGCATATTCCATCGTTTTTCCTTCTTACTATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAt                                                                                              <  2:3679014/97‑1 (MQ=255)

          aTGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                                           >  2:3165026/1‑100 (MQ=255)

           tGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCg                                                                                          <  2:828023/100‑1 (MQ=255)

           tGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCg                                                                                          >  2:414461/1‑100 (MQ=255)

            gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                         <  2:3597335/100‑1 (MQ=255)

            gTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                         <  2:3614100/100‑1 (MQ=255)

             ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCggn                                                                                        <  2:1541762/100‑2 (MQ=38)

             ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGt                                                                                        <  2:1654051/100‑1 (MQ=38)

             ttCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGt                                                                                        <  2:1546523/100‑1 (MQ=38)

              tCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGtn                                                                                       <  1:1226116/100‑2 (MQ=37)

               cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTAc                                                                                      >  2:3397049/1‑100 (MQ=255)

               cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaa                                                                                      <  1:817959/100‑2 (MQ=37)

               cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaa                                                                                      <  2:3279350/100‑2 (MQ=37)

                ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                                                     <  2:461241/100‑1 (MQ=255)

                ccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaa                                                                                     <  1:1014626/100‑3 (MQ=37)

                ccTCCGTCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCTTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                                       <  1:1710916/98‑1 (MQ=255)

                  tGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatc                                                                                   <  2:1125792/100‑5 (MQ=25)

                    cgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACagag                                                                                  >  1:2518828/1‑99 (MQ=255)

                    cgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGc                                                                                 <  2:145903/100‑1 (MQ=255)

                    cgcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctg                                                                                 <  1:2087724/100‑7 (MQ=25)

                     gcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCAGGTCGGTACAGAGCa                                                                                >  2:1906250/1‑100 (MQ=255)

                     gcCAGCGCGACCAGAAAAGCATATACCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctgg                                                                                <  2:1374484/100‑8 (MQ=25)

                      ccAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCaa                                                                               >  1:532261/1‑100 (MQ=255)

                      ccAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCaa                                                                               >  1:1904497/1‑100 (MQ=255)

                       cAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctggtc                                                                              <  1:1858161/100‑10 (MQ=21)

                       cAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatctggtc                                                                              <  1:4269559/100‑10 (MQ=21)

                         gcgcGACCAGAAAAGCATATTCCATCGCTACGCCTTAGTACGAGTTAAAGCAACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaa                                                                            >  1:2382811/1‑100 (MQ=255)

                           gtgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                          <  2:4397843/98‑1 (MQ=255)

                           gcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                          <  2:3155551/100‑1 (MQ=255)

                           gcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                          <  2:4282756/100‑1 (MQ=255)

                            cgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACGGCCATGACCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                                                         >  1:3090915/1‑100 (MQ=255)

                            cgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                                                         >  1:3223100/1‑100 (MQ=255)

                             gACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGa                                                                        >  2:2265545/1‑100 (MQ=255)

                               ccAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                      <  2:3760162/100‑1 (MQ=255)

                                cAGAAAAGCATATTCCATCGCTACGACTTCGTACGATTTAAAGCGACCAGATTTACCGCCAAGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                                     <  1:2627481/100‑1 (MQ=255)

                                 aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                                    >  1:4013278/1‑100 (MQ=255)

                                      aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGtacggc                                                               >  2:1728919/1‑95 (MQ=255)

                                        cATATTCCATCGCTACGCCTTCGTACGATTTAAGACGACCAGATTTACCGCCATGGCCTGAGTCCATGCCGGTAAAGAACAATAAAAGCCGCTCATCGGt                                                             >  1:4391296/1‑100 (MQ=255)

                                          tatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAAAGCAATAAAAGATGGTCATCGGttt                                                           >  2:2894410/1‑100 (MQ=255)

                                           atTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACCGAGCAATAAAAGATGGTAATCGGtttt                                                          >  1:4188980/1‑100 (MQ=255)

                                            ttacATCGCTACGCGTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                         <  1:3549581/97‑1 (MQ=255)

                                             tCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                                        <  1:1001494/100‑1 (MQ=255)

                                               cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                      >  1:905662/1‑100 (MQ=255)

                                               cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                      >  1:2858849/1‑100 (MQ=255)

                                               cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                      <  1:3162338/100‑1 (MQ=255)

                                                 tCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTc                                                    <  1:2921580/100‑1 (MQ=255)

                                                         ccTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAtn                                            <  1:565980/100‑2 (MQ=255)

                                                           ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTa                                          >  2:4225196/1‑100 (MQ=255)

                                                              gTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg                                       <  2:2552680/100‑1 (MQ=255)

                                                               tACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGa                                      >  1:3069804/1‑100 (MQ=255)

                                                                ncGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAACGATGGTCATCGGCTTTCAGCTCGCGCAATTTAGCGAc                                     >  1:1226103/2‑100 (MQ=255)

                                                                  gATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGAccc                                   >  1:2467751/1‑100 (MQ=255)

                                                                   aTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCa                                  >  1:1418886/1‑100 (MQ=255)

                                                                    tttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAt                                 >  2:4224635/1‑100 (MQ=255)

                                                                      tAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAttt                               >  1:4367691/1‑100 (MQ=255)

                                                                      tAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAttt                               >  2:2134325/1‑100 (MQ=255)

                                                                         aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                            >  2:4077247/1‑100 (MQ=255)

                                                                         aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                            <  1:230596/100‑1 (MQ=255)

                                                                            gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                         <  2:3877157/100‑1 (MQ=255)

                                                                            gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                         <  2:3076171/100‑1 (MQ=255)

                                                                            gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                         >  2:3530943/1‑100 (MQ=255)

                                                                            gACCAGATTTACCGCCAGGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                         <  2:1205341/100‑1 (MQ=255)

                                                                             aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                         <  2:2753876/99‑1 (MQ=255)

                                                                             aCCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGt                        >  1:1518401/1‑100 (MQ=255)

                                                                                aGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGGGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTc                      >  1:2350477/1‑99 (MQ=255)

                                                                                aGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc                     >  2:3796116/1‑100 (MQ=255)

                                                                                 tctgtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                    >  1:3327067/5‑100 (MQ=25)

                                                                                 gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                    <  2:2015586/100‑1 (MQ=255)

                                                                                 gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                    <  2:3061583/100‑1 (MQ=255)

                                                                                 gATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc                     >  1:1349070/1‑99 (MQ=255)

                                                                                   tttACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCaa                  >  2:1054586/1‑100 (MQ=255)

                                                                                   tgtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCAGGTTcac                    >  2:3018956/3‑96 (MQ=25)

                                                                                      aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat               >  2:1670430/1‑100 (MQ=38)

                                                                                      aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat               >  2:3496408/1‑100 (MQ=38)

                                                                                      aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAatat               >  2:545684/1‑100 (MQ=38)

                                                                                       ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt              >  1:2332118/1‑100 (MQ=255)

                                                                                       ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt              <  1:1650476/100‑1 (MQ=255)

                                                                                       ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt              >  2:4345037/1‑100 (MQ=38)

                                                                                       ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt              >  1:1835078/1‑100 (MQ=38)

                                                                                       ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt              >  1:3968915/1‑100 (MQ=38)

                                                                                       ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATAtt              >  1:1871828/1‑100 (MQ=38)

                                                                                        cGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg             <  1:1853513/100‑1 (MQ=255)

                                                                                         gCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc            <  1:4206091/100‑1 (MQ=255)

                                                                                         gCCATGGCCTGAGTCCACGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc            <  2:3082661/100‑1 (MQ=255)

                                                                                           cATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCan          <  2:4004494/100‑2 (MQ=255)

                                                                                            aTGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACTCATTTTGCCGGTTCCCAATATTGCAcc         <  1:690010/100‑1 (MQ=255)

                                                                                             tGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCt        <  2:1309259/100‑1 (MQ=255)

                                                                                              ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGTAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTg       <  1:2415829/100‑1 (MQ=255)

                                                                                              ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTg       <  1:4224604/100‑1 (MQ=255)

                                                                                                ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag     >  1:2694967/1‑100 (MQ=255)

                                                                                                ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag     <  1:2034500/100‑1 (MQ=255)

                                                                                                  tGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGaa   <  1:292397/100‑1 (MQ=255)

GTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAAT  >  NC_000913/1926792‑1926990

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