BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,536,994	T→A	S204C (AGC→TGC)	cysM ←	cysteine synthase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,536,994	0	T	A	91.4%	737.0 / 49.5	223	S204C (AGC→TGC)	cysM	cysteine synthase B
Reads supporting (aligned to +/- strand): ref base T (6/13); new base A (95/108); total (102/121)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 235 total aligned reads displayed.

GATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCC  >  NC_000913/2536904‑2537083

  tgaATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGccc                                                                                       >  1:4851378/1‑93 (MQ=255)

  tgaATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGccc                                                                                       >  1:6745447/1‑93 (MQ=255)

  tgaATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGccc                                                                                       <  1:5213918/93‑1 (MQ=255)

    aaTATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCctc                                                                                     <  2:4096036/93‑1 (MQ=255)

    aaTATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCctc                                                                                     <  1:5905679/93‑1 (MQ=255)

     atatCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCctct                                                                                    <  1:7136111/93‑1 (MQ=255)

     atatCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCctct                                                                                    <  1:7710973/93‑1 (MQ=255)

       atCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTc                                                                                  <  2:3227369/93‑1 (MQ=255)

       atCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTc                                                                                  <  1:8535369/93‑1 (MQ=255)

        tCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTcc                                                                                 <  2:660291/93‑1 (MQ=255)

        tCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTcc                                                                                 >  1:1665369/1‑93 (MQ=255)

         ccAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCg                                                                                >  2:8167102/1‑93 (MQ=255)

          cAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCgg                                                                               <  1:3547319/93‑1 (MQ=255)

              aCCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGc                                                                           <  2:5884103/93‑1 (MQ=255)

              aCCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGc                                                                           >  2:3685440/1‑93 (MQ=255)

               ccTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCa                                                                          >  1:4655254/1‑93 (MQ=255)

               ccTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCa                                                                          <  2:6463591/93‑1 (MQ=255)

               ccTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCa                                                                          <  1:1186131/93‑1 (MQ=255)

               ccTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCa                                                                          <  2:6771423/93‑1 (MQ=255)

               ccTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCa                                                                          <  1:1557782/93‑1 (MQ=255)

                cTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAg                                                                         >  1:3282573/1‑93 (MQ=255)

                 tCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAgg                                                                        >  1:3395010/1‑93 (MQ=255)

                 tCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAgg                                                                        >  2:640884/1‑93 (MQ=255)

                  cATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGc                                                                       <  2:3359670/93‑1 (MQ=255)

                   aTCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGcc                                                                      <  2:3480608/93‑1 (MQ=255)

                   aTCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGcc                                                                      >  1:2490586/1‑93 (MQ=255)

                      caccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGAc                                                                   <  2:715926/93‑1 (MQ=255)

                         caGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAAt                                                                <  2:8064771/93‑1 (MQ=255)

                           gagaAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATgg                                                              >  1:5366954/1‑93 (MQ=255)

                           gagaAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATgg                                                              <  1:1336470/93‑1 (MQ=255)

                           gagaAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATgg                                                              >  1:6042394/1‑93 (MQ=255)

                           gagaAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATgg                                                              >  1:9709793/1‑93 (MQ=255)

                             gaAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTc                                                            >  2:9236637/1‑93 (MQ=255)

                               aGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCAc                                                          >  2:4328078/1‑93 (MQ=255)

                               aGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCAc                                                          <  1:6369731/93‑1 (MQ=255)

                                gCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCAcc                                                         <  2:1690081/93‑1 (MQ=255)

                                gCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCAcc                                                         <  1:7731361/93‑1 (MQ=255)

                                 cGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCg                                                        <  1:1012752/93‑1 (MQ=255)

                                   ttGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGGCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGt                                                      >  1:6684808/1‑93 (MQ=255)

                                   ttGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGt                                                      >  1:4169422/1‑93 (MQ=255)

                                   ttGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGt                                                      >  1:4061614/1‑93 (MQ=255)

                                    tGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGtt                                                     <  1:7293805/93‑1 (MQ=255)

                                    tGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGtt                                                     <  1:7523373/93‑1 (MQ=255)

                                     gAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGttt                                                    <  1:7933805/93‑1 (MQ=255)

                                      aaaaTCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTg                                                   <  1:5066229/93‑1 (MQ=255)

                                        aaTCCCCGGCAGATATTCCGTAGGCCAGCGGCGCATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGa                                                 <  2:8955518/93‑1 (MQ=255)

                                         aTCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGAt                                                <  1:9373410/93‑1 (MQ=255)

                                          tCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGAtt                                               <  1:6911661/93‑1 (MQ=255)

                                             ccGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGtt                                            <  1:1405286/93‑1 (MQ=255)

                                              cGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTc                                           <  2:7435837/93‑1 (MQ=255)

                                              cGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTc                                           <  1:9318835/93‑1 (MQ=255)

                                              cGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTc                                           <  1:3788967/93‑1 (MQ=255)

                                               ggCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTcg                                          >  1:3369179/1‑93 (MQ=255)

                                               ggCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTcg                                          <  2:9548407/93‑1 (MQ=255)

                                                 cAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTcgcg                                        >  2:4844018/1‑93 (MQ=255)

                                                   gATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCa                                      >  1:4264998/1‑93 (MQ=255)

                                                   gATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCa                                      <  1:2735797/93‑1 (MQ=255)

                                                   gATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCa                                      >  1:3817425/1‑93 (MQ=255)

                                                    atatTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCAt                                     <  1:5619602/93‑1 (MQ=255)

                                                     tatTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATa                                    <  2:9341299/93‑1 (MQ=255)

                                                      ggtCCGTAGGGGAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGGCGACAATGGTCCCCGGTTTGGATTGTTCGCGCATaa                                   <  2:5275004/91‑1 (MQ=255)

                                                      atTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATaa                                   >  2:8361352/1‑93 (MQ=255)

                                                       ttCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATaaa                                  >  1:7381148/1‑93 (MQ=255)

                                                        tCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAg                                 <  2:1646574/93‑1 (MQ=255)

                                                          cGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCg                               <  1:8696909/93‑1 (MQ=255)

                                                          cGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCg                               >  1:9905052/1‑93 (MQ=255)

                                                          cGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCg                               >  2:357718/1‑93 (MQ=255)

                                                          cGTAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCg                               >  2:3740174/1‑93 (MQ=255)

                                                            tAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTg                             <  2:6882457/93‑1 (MQ=255)

                                                            tAGGCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTg                             <  2:2539075/93‑1 (MQ=255)

                                                               gCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTgaga                          <  2:7614833/93‑1 (MQ=255)

                                                               gCCAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTgaga                          <  1:9259052/93‑1 (MQ=255)

                                                                ccAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGAc                         <  2:9472125/93‑1 (MQ=255)

                                                                ccAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGAc                         >  1:4324827/1‑93 (MQ=255)

                                                                 cAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACg                        <  2:7793711/93‑1 (MQ=255)

                                                                 cAGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACg                        >  2:743029/1‑93 (MQ=255)

                                                                  aGCGGCGAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGc                       >  2:9162138/1‑93 (MQ=255)

                                                                     ggcgAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCgg                    >  1:8751976/1‑93 (MQ=255)

                                                                     ggcgAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCgg                    <  2:1210252/93‑1 (MQ=255)

                                                                     ggcgAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACACCgg                    >  1:5989393/1‑93 (MQ=255)

                                                                      gcgAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCggg                   >  2:3112406/1‑92 (MQ=255)

                                                                      gcgAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGt                   >  2:243612/1‑93 (MQ=255)

                                                                       cgAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTg                  <  2:4869130/93‑1 (MQ=255)

                                                                       cgAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTg                  >  1:2418308/1‑93 (MQ=255)

                                                                       cgAATGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTg                  >  1:1977340/1‑93 (MQ=255)

                                                                          aTGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATa               >  2:7661938/1‑93 (MQ=255)

                                                                          aTGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATa               <  1:6822358/93‑1 (MQ=255)

                                                                           tGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAg              <  2:2049345/93‑1 (MQ=255)

                                                                           tGCCGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAg              >  2:8711501/1‑93 (MQ=255)

                                                                             ccGGGAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTg            >  1:8284252/1‑93 (MQ=255)

                                                                               gggAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGcc          <  2:4919475/93‑1 (MQ=255)

                                                                               ggaaaTGCTGCAGCCCTCTTCCGGTGGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGcc          <  2:131477/90‑1 (MQ=255)

                                                                                 gAATGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCgg        >  1:4537840/1‑93 (MQ=255)

                                                                                  aaTGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGgt       <  2:9663018/93‑1 (MQ=255)

                                                                                   aTGCTGCAGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGgtc      >  1:8124033/1‑93 (MQ=255)

                                                                                    agctgcagcCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGgtcg     <  1:9983684/92‑1 (MQ=255)

                                                                                     gctgcagcCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGgtcgt    >  1:9598191/1‑93 (MQ=255)

                                                                                      ctgcagcCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTTCCGgtcgtc   >  1:2864051/1‑93 (MQ=255)

                                                                                       tgcagcCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTcc  >  1:3055170/1‑93 (MQ=255)

                                                                                       tgcagcCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTcc  <  2:5840903/93‑1 (MQ=255)

GATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCC  >  NC_000913/2536904‑2537083

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