BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,531,984	C→T	G67G (GGC→GGT)	ptsH →	PTS system protein HPr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,531,984	0	C	T	99.2%	842.9 / ‑5.5	240	G67G (GGC→GGT)	ptsH	PTS system protein HPr
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (119/118); minor base G (2/0); total (122/118)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 245 total aligned reads displayed.

CAACGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAG  >  NC_000913/2531894‑2532073

cAACGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGa                                                                                         <  2:1578279/93‑1 (MQ=255)

 aaCGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGaa                                                                                        >  2:5160386/1‑93 (MQ=255)

 aaCGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGaa                                                                                        <  1:6233068/93‑1 (MQ=255)

 aaCGACAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGaa                                                                                        <  1:4174372/93‑1 (MQ=255)

  aCGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAg                                                                                       >  1:6726233/1‑93 (MQ=255)

  aCGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAg                                                                                       >  2:2562531/1‑93 (MQ=255)

   cGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGa                                                                                      <  2:6036779/93‑1 (MQ=255)

   cGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGa                                                                                      >  1:1066402/1‑93 (MQ=255)

   cGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGa                                                                                      <  2:6502836/93‑1 (MQ=255)

    ggCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGAc                                                                                     >  1:4493375/1‑93 (MQ=255)

     gCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACg                                                                                    >  2:6605905/1‑93 (MQ=255)

       aaaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGac                                                                                  >  1:4739793/1‑92 (MQ=255)

       aaaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAg                                                                                  <  1:1380312/93‑1 (MQ=255)

       aaaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAg                                                                                  >  1:1129119/1‑93 (MQ=255)

        aaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGc                                                                                 >  2:3090197/1‑93 (MQ=255)

        aaaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGc                                                                                 >  2:4295646/1‑93 (MQ=255)

         aaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCa                                                                                >  2:468236/1‑93 (MQ=255)

         aaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCa                                                                                <  1:4859977/93‑1 (MQ=255)

         aaGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCa                                                                                <  1:4976838/93‑1 (MQ=255)

           gcgcCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGa                                                                              <  1:269546/93‑1 (MQ=255)

             gcCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGaaa                                                                            <  1:1157788/93‑1 (MQ=255)

             gcCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGaaa                                                                            <  2:126181/93‑1 (MQ=255)

              ccAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAg                                                                           <  1:4649179/93‑1 (MQ=255)

              ccAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAg                                                                           >  1:2235318/1‑93 (MQ=255)

                  cgcgAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGt                                                                       <  1:667863/93‑1 (MQ=255)

                  cgcgAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGt                                                                       <  2:2225257/93‑1 (MQ=255)

                   gcgAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGtt                                                                      <  1:5258176/93‑1 (MQ=255)

                      aaaaGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGaa                                                                   >  1:361896/1‑93 (MQ=255)

                       aaaGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAAc                                                                  >  1:809781/1‑93 (MQ=255)

                        aaGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACa                                                                 >  1:346513/1‑93 (MQ=255)

                         aGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACAt                                                                >  1:4255561/1‑93 (MQ=255)

                         aGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACAt                                                                >  1:2194716/1‑93 (MQ=255)

                          gCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATc                                                               <  1:5998179/93‑1 (MQ=255)

                             tGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTgg                                                            >  1:2164735/1‑93 (MQ=255)

                             tGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTgg                                                            <  2:5103940/93‑1 (MQ=255)

                             tGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGGTGAACATCTgg                                                            >  1:2893386/1‑93 (MQ=255)

                              gTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGt                                                           <  1:5739964/93‑1 (MQ=255)

                               tttAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGtt                                                          <  2:2107328/93‑1 (MQ=255)

                               tttAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGtt                                                          >  1:1690372/1‑93 (MQ=255)

                                   aaCTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAAc                                                      <  1:6351190/93‑1 (MQ=255)

                                   aaCTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGGGAAGACGAGCAGCAAGGGGTTAGACATCTGGTTAAAc                                                      >  2:1874054/1‑93 (MQ=255)

                                       gCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGAt                                                  >  1:874275/1‑93 (MQ=255)

                                       gCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGAt                                                  >  2:1012969/1‑93 (MQ=255)

                                       gCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGAt                                                  <  2:3671676/93‑1 (MQ=255)

                                        cAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATg                                                 <  1:3423291/93‑1 (MQ=255)

                                        cAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATg                                                 <  2:625015/93‑1 (MQ=255)

                                        cAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATg                                                 >  1:1436449/1‑93 (MQ=255)

                                        cAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATg                                                 >  1:3029497/1‑93 (MQ=255)

                                        cAGACTCTGGGACTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATg                                                 <  1:5644302/93‑1 (MQ=255)

                                         aGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATgg                                                <  1:163402/93‑1 (MQ=255)

                                         aGACTCTGGGCCTGACTCAAGGTACCGTTGGGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATgg                                                >  2:765942/1‑93 (MQ=255)

                                          gACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTACACTGATGGc                                               >  1:2416753/1‑93 (MQ=255)

                                             tctGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGa                                            >  1:5821734/1‑93 (MQ=255)

                                              ctGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGaa                                           <  1:4363523/93‑1 (MQ=255)

                                               tGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAAc                                          <  1:4801312/93‑1 (MQ=255)

                                                   ccTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGa                                      >  2:6361299/1‑93 (MQ=255)

                                                    tagtcTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAg                                     >  1:6481580/5‑93 (MQ=255)

                                                    cTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAg                                     <  1:2829273/93‑1 (MQ=255)

                                                     tGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAga                                    >  2:5223764/1‑92 (MQ=255)

                                                     tGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGt                                    >  2:3393956/1‑93 (MQ=255)

                                                     tGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGt                                    >  1:4568546/1‑93 (MQ=255)

                                                      gACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTa                                   >  1:2143791/1‑93 (MQ=255)

                                                      gACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTa                                   >  1:1803906/1‑93 (MQ=255)

                                                      gACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTa                                   >  2:3225120/1‑93 (MQ=255)

                                                        cTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGCACTCGAGTAAt                                 >  2:5656396/1‑93 (MQ=255)

                                                        cTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAAt                                 >  1:5902876/1‑93 (MQ=255)

                                                        cTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAAt                                 <  2:599960/93‑1 (MQ=255)

                                                         tCAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAAtt                                >  2:875404/1‑93 (MQ=255)

                                                          cAAGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAAttt                               >  1:6341408/1‑93 (MQ=255)

                                                           aaNGNANCGTTGTGACTATCTCNGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCNGAACTCGAGTAATTTc                              <  2:4027462/93‑1 (MQ=255)

                                                           aaGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTc                              >  1:5366255/1‑93 (MQ=255)

                                                            aGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTcc                             >  2:1975843/1‑93 (MQ=255)

                                                            aGGTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTcc                             >  1:3535752/1‑93 (MQ=255)

                                                              gTACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAAATTCCCg                           <  1:5554867/93‑1 (MQ=255)

                                                               tACCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCgg                          <  1:492155/93‑1 (MQ=255)

                                                                aCCGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCggg                         <  2:3588944/93‑1 (MQ=255)

                                                                 ccGTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGt                        >  2:143795/1‑93 (MQ=255)

                                                                   gTTGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTc                      <  1:3062124/93‑1 (MQ=255)

                                                                    ttGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCt                     >  2:965023/1‑93 (MQ=255)

                                                                    ttGTGACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCt                     >  1:6529106/1‑93 (MQ=255)

                                                                     tgtgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCtt                    <  2:5412803/93‑1 (MQ=255)

                                                                     tgtgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTAAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCtt                    <  1:3951304/93‑1 (MQ=255)

                                                                      gtgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACGTCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCttt                   <  2:2784300/93‑1 (MQ=255)

                                                                       tgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCtttt                  >  2:5521475/1‑93 (MQ=255)

                                                                       tgACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCtttt                  <  1:5166589/93‑1 (MQ=255)

                                                                        gACTATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTa                 >  2:5906156/1‑93 (MQ=255)

                                                                           tATCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTaaaa              >  2:5834463/1‑93 (MQ=255)

                                                                            aTCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTaaaaa             <  2:5000220/93‑1 (MQ=255)

                                                                            aTCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTaaaaa             >  1:6075592/1‑93 (MQ=255)

                                                                            aTCTCCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTaaaaa             <  1:46064/93‑1 (MQ=255)

                                                                              ctcCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATc           <  2:6090620/93‑1 (MQ=255)

                                                                              ctcCGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATc           <  1:2776385/93‑1 (MQ=255)

                                                                                 cGCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGt        <  2:4447612/93‑1 (MQ=255)

                                                                                  gCAGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTc       >  2:3885039/1‑93 (MQ=255)

                                                                                    aGAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTcac     >  1:5269788/1‑93 (MQ=255)

                                                                                     gAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTcaca    >  2:1415151/1‑93 (MQ=255)

                                                                                     gAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTcaca    >  2:1460413/1‑93 (MQ=255)

                                                                                     gAAGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTcaca    >  1:1937132/1‑93 (MQ=255)

                                                                                      aaGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACaa   <  1:3868292/93‑1 (MQ=255)

                                                                                       aGGTGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAg  >  2:1673389/1‑93 (MQ=255)

CAACGGCAAAAGCGCCAGCGCGAAAAGCCTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAG  >  NC_000913/2531894‑2532073

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