BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	850,742	A→G	23.3%	D98G (GAC→GGC)	ompX →	outer membrane protein X

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	850,742	0	A	G	23.3%	214.8 / 69.9	124	D98G (GAC→GGC)	ompX	outer membrane protein X
Reads supporting (aligned to +/- strand): ref base A (54/41); new base G (21/8); total (75/49)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 126 total aligned reads displayed.

CTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCT  >  NC_000913/850643‑850840

cTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGa                                                                                                    >  2:496600/1‑100 (MQ=255)

 tttCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGAc                                                                                                   >  2:1009531/1‑100 (MQ=255)

  ttCACTTACACCGAGAAAAGCCGTACTGCACGCTCTTGTGACAACAACAAAAACCCGTCCTACGGCATTACCGACGGTGCGACTACCCGCATTAAAGa                                                                                                    >  1:1033567/1‑98 (MQ=255)

    cACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTgg                                                                                                >  2:2419123/1‑100 (MQ=255)

    cACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTgg                                                                                                <  2:1310969/100‑1 (MQ=255)

    cACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGCTACCTCCAACAAAAACCAGCAGCACACCACCACTGCTGCTCCGGCTTACCGCATTAACGACTgg                                                                                                >  1:719256/1‑100 (MQ=255)

      cgtACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGc                                                                                              <  2:1482605/98‑1 (MQ=255)

      cTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGc                                                                                              <  1:2534883/100‑1 (MQ=255)

      cTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGc                                                                                              <  2:3303263/100‑1 (MQ=255)

        tACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCaa                                                                                            >  2:734569/1‑100 (MQ=255)

          cacCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGc                                                                                          >  2:458777/1‑100 (MQ=255)

            ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCa                                                                                         >  1:487891/1‑99 (MQ=255)

            ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCAt                                                                                        >  2:3241937/1‑100 (MQ=255)

            ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCAt                                                                                        >  1:1710319/1‑100 (MQ=255)

            ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGCCAACACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGc                                                                                          >  1:1116688/1‑98 (MQ=255)

              gagaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCt                                                                                      <  2:310532/100‑1 (MQ=255)

               agaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTa                                                                                     <  1:1239140/100‑1 (MQ=255)

                 aaaaGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACg                                                                                   >  1:1595292/1‑100 (MQ=255)

                   aaGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGgt                                                                                 >  1:808029/1‑100 (MQ=255)

                     gCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGgtgt                                                                               >  1:1819449/1‑100 (MQ=255)

                      ccGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTa                                                                              >  2:179857/1‑100 (MQ=255)

                          aCTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgg                                                                          >  1:1645432/1‑100 (MQ=255)

                           cTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtggg                                                                         >  2:860937/1‑100 (MQ=255)

                           cTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtggg                                                                         <  2:2317633/100‑1 (MQ=255)

                           cTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtggg                                                                         >  2:2648213/1‑100 (MQ=255)

                             gCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtg                                                                       >  1:1025739/1‑100 (MQ=255)

                             gCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtg                                                                       >  1:1408711/1‑100 (MQ=255)

                             gCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtg                                                                       >  1:2080820/1‑100 (MQ=255)

                              cAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgt                                                                      <  2:3027246/100‑1 (MQ=255)

                               aaGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtgtg                                                                     <  2:464236/100‑1 (MQ=255)

                               aaGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtg                                                                     <  1:496598/100‑1 (MQ=255)

                                 gCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtggg                                                                   <  1:1931871/100‑1 (MQ=255)

                                 gCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtggg                                                                   <  1:562354/100‑1 (MQ=255)

                                  ctctGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtgtgggt                                                                  >  2:2124558/1‑100 (MQ=255)

                                  ctctGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtgggt                                                                  <  2:2021053/100‑1 (MQ=255)

                                    ctGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTa                                                                >  1:1495323/1‑100 (MQ=255)

                                    ctGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTa                                                                <  1:1714344/100‑1 (MQ=255)

                                      ggTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATg                                                              >  2:2671237/1‑100 (MQ=255)

                                      ggTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATg                                                              >  2:1676258/1‑100 (MQ=255)

                                       gTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATgg                                                             <  2:1128493/100‑1 (MQ=255)

                                         gACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTa                                                           >  2:1612023/1‑100 (MQ=255)

                                         atcTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTa                                                           <  2:2307897/98‑1 (MQ=255)

                                           cTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa                                                         >  1:2897433/1‑100 (MQ=255)

                                           cTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa                                                         >  1:276279/1‑100 (MQ=255)

                                             acaacaAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTAGGTTATGGTaagat                                                       >  1:2307881/1‑97 (MQ=255)

                                              caacaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTc                                                      >  2:102701/1‑100 (MQ=255)

                                               aacaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTcc                                                     >  1:1363092/1‑100 (MQ=255)

                                                acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa                                                    >  2:2165610/1‑100 (MQ=255)

                                                acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGGGGGTGTGGGTTATGGTAAAt                                                        >  1:2810303/1‑96 (MQ=255)

                                                 caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg                                                   >  2:95210/1‑100 (MQ=255)

                                                 caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg                                                   <  2:888076/100‑1 (MQ=255)

                                                 caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg                                                   >  1:739415/1‑100 (MQ=255)

                                                 caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg                                                   <  2:2543890/100‑1 (MQ=255)

                                                 caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg                                                   <  2:413825/100‑1 (MQ=255)

                                                 caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg                                                   <  2:1741855/100‑1 (MQ=255)

                                                   aaaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAc                                                 >  2:1108800/1‑100 (MQ=255)

                                                    aaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAcc                                                <  2:1113296/100‑1 (MQ=255)

                                                     aaCCAGTACTACGGCATCACTGCTGGTCCGGCATACCGCATTAACGACTGGGCAAGCATCTAAGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGa                                                  >  1:322980/1‑97 (MQ=255)

                                                       tccgatctacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACt                                             <  2:1274137/94‑1 (MQ=255)

                                                         aGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGa                                           <  1:432708/100‑1 (MQ=255)

                                                          gTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGaa                                          >  2:2262207/1‑100 (MQ=255)

                                                            actacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATa                                        >  2:466218/1‑100 (MQ=255)

                                                              tacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATAcc                                      <  1:2455895/100‑1 (MQ=255)

                                                              ccgatcttCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATAcn                                      <  2:3084950/93‑2 (MQ=255)

                                                               acGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATAccc                                     <  1:2233514/100‑1 (MQ=255)

                                                               acGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACtagatcgga                                     >  1:1274133/1‑92 (MQ=255)

                                                               acGGCATCACTGCGGGGCCGGCGTACCGCATTAAGGGGTGGGCATGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATAccc                                     >  1:2700956/1‑100 (MQ=255)

                                                                cGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCg                                    >  2:2537165/1‑100 (MQ=255)

                                                                   cATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGAcc                                 <  2:2720218/100‑1 (MQ=255)

                                                                     tCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTa                               >  2:2820057/1‑100 (MQ=255)

                                                                     ncACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATAccagatcgg                               >  1:3084924/2‑93 (MQ=255)

                                                                      cACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTAc                              <  1:2425156/100‑1 (MQ=255)

                                                                        ctgctgGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACaa                            <  1:1102521/100‑1 (MQ=255)

                                                                        ctgctgGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACaa                            >  1:797523/1‑100 (MQ=255)

                                                                             ggTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                       >  1:3307407/1‑100 (MQ=255)

                                                                             ggTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                       <  1:2676689/100‑1 (MQ=255)

                                                                             ggTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                       <  2:1684396/100‑1 (MQ=255)

                                                                               tCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGac                     <  2:1065705/100‑1 (MQ=255)

                                                                               tCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGac                     <  1:2672530/100‑1 (MQ=255)

                                                                                  ggCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACAcc                  <  1:2424174/100‑1 (MQ=255)

                                                                                    cTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAg                >  1:2630/1‑100 (MQ=255)

                                                                                    cTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAg                >  2:1160185/1‑100 (MQ=255)

                                                                                     ttACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGc               <  2:1819460/100‑1 (MQ=255)

                                                                                      tACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCg              >  2:2208519/1‑100 (MQ=255)

                                                                                      tACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCg              <  2:3247400/100‑1 (MQ=255)

                                                                                       aCCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGa             >  2:1828404/1‑100 (MQ=255)

                                                                                       aCCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGa             >  1:1450259/1‑100 (MQ=255)

                                                                                       aCCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGa             >  1:515589/1‑100 (MQ=255)

                                                                                       aCCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGGAAATTCCAGACCACTGAATACCCGGCCTACAAACACGACACCAGCGa             >  1:2874623/1‑100 (MQ=255)

                                                                                        ccGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGAc            >  2:637851/1‑100 (MQ=255)

                                                                                        ccGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGAc            <  2:1508730/100‑1 (MQ=255)

                                                                                         cGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACt           <  2:1606989/100‑1 (MQ=255)

                                                                                          gCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTa          >  1:255042/1‑100 (MQ=255)

                                                                                          gCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTa          >  1:2607438/1‑100 (MQ=255)

                                                                                           cATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTAc         <  1:2262191/100‑1 (MQ=255)

                                                                                               aaCGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGtt     <  2:643116/100‑1 (MQ=255)

                                                                                               aaCGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGtt     <  2:73792/100‑1 (MQ=255)

                                                                                               aaCGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGtt     >  2:2669869/1‑100 (MQ=255)

                                                                                                aCGGCTGTGCAAGCATCTACGGTGTAGGGGCTGTGGGTTATGGTAAATGCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGttt    <  1:2747526/100‑1 (MQ=255)

                                                                                                  gACTGGGCAAGCATCTACGGTGTAGTAGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTtct  >  2:2057278/1‑100 (MQ=255)

CTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTTTCT  >  NC_000913/850643‑850840

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