BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)_6→7	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	99.1%	376.5 / ‑3.3	106	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (53/52); total (54/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 117 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTG  >  NC_000913/4638088‑4638280

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                 <  2:918229/100‑1 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                               <  1:2766271/100‑1 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                               >  2:761068/1‑100 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                               >  2:1452474/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c                                                                                             >  1:114667/1‑100 (MQ=255)

   tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cc                                                                                            >  1:2287649/1‑100 (MQ=255)

     acTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccc                                                                                          <  1:1149564/99‑1 (MQ=255)

       ttagACTTCACCGTCACTTCACCTAGCTTCAAATTCTTCCTACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGATAAGACTATTGTGAAAT‑c‑ccccc                                                                                         >  2:1977305/4‑99 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccc                                                                                        <  2:920985/100‑2 (MQ=255)

        tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑cccccct                                                                                       >  1:1557180/1‑98 (MQ=255)

             ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCt                                                                                   <  1:491285/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTg                                                                                  <  1:2887171/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTg                                                                                  >  2:517843/1‑100 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTgg                                                                                 <  1:289008/100‑1 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGaa                                                                               >  2:2684796/1‑100 (MQ=255)

                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGaaa                                                                              <  2:50263/100‑1 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAAt                                                                            <  2:2140575/100‑1 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGa                                                                          >  1:2593753/1‑100 (MQ=255)

                        ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACt                                                                        <  2:2153633/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                    <  2:36459/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                    <  1:182514/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                    >  2:2852140/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                    <  2:94176/100‑1 (MQ=255)

                            caTAGCGTCAAATTCTTCCCACATAGTCTTCGTATCCGGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                    <  2:2479887/100‑1 (MQ=255)

                              tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTg                                                                  <  1:1017249/100‑1 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                               >  1:2147892/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CACCCCCTGTACTGGAAAATGACTAGCCTGttt                                                               >  2:1943968/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                              >  2:2175174/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                              >  2:1150604/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                              >  2:760293/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                             <  1:2315444/100‑1 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                             <  2:580113/100‑1 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                            >  1:2107933/1‑100 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                           >  2:2611996/1‑100 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                           >  2:2125452/1‑100 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                          >  1:231484/1‑100 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCan                                                        <  1:1738531/100‑2 (MQ=255)

                                            tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAtt                                                    >  1:1643113/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                   >  1:2131927/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                   <  2:270590/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                  >  2:2383249/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 <  1:2452732/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  2:345206/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTCTTTGGTGCAGTATTGc                                                  >  2:2263257/1‑99 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                >  2:2301121/1‑100 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                             <  2:1392082/100‑1 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                            >  1:2516067/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           <  1:2684816/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           <  2:2433371/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           >  1:2317189/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           <  2:2736701/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                          >  1:1755486/1‑100 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                        <  1:789914/100‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                        <  2:2361863/100‑1 (MQ=255)

                                                         tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                       <  1:72196/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                  <  1:1868672/100‑1 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                 >  2:1718369/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                >  2:2094278/1‑100 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                               >  1:640727/1‑100 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                               >  1:492970/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                              <  1:332493/100‑1 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                              <  1:983686/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                             <  2:812200/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                            <  2:1180483/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                            <  1:1373339/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                            <  2:1976206/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                            <  2:1976282/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                            <  1:726167/100‑1 (MQ=255)

                                                                     gcCATTGCAAAGGAAAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                           <  2:2670017/100‑1 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                        >  2:2844211/1‑100 (MQ=255)

                                                                         ttGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                       >  1:2207121/1‑100 (MQ=255)

                                                                         ttGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                       >  1:663671/1‑100 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                      <  1:2023757/100‑1 (MQ=255)

                                                                           gCAAAGGAGCAGACTATGTTGAAATCC‑CCCCCTGTTATGGAAAATGCCTAGCCTGTTTGGTGCAGTCTTGCTGTTGTTGATTCCGATAATGCTGATTCg                      <  1:2706858/99‑1 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                     >  1:2149519/1‑100 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                     >  2:859750/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    <  1:1717559/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    >  1:361598/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    <  2:2107920/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    <  2:361597/100‑1 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                  >  2:1116037/1‑100 (MQ=255)

                                                                               aGGAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                 >  2:151757/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                >  2:941803/1‑100 (MQ=255)

                                                                                 gagaAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg               <  1:2214659/100‑1 (MQ=255)

                                                                                 gagaAGACTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg               <  1:1452481/100‑1 (MQ=255)

                                                                                    aagaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAtt            <  2:1223224/100‑1 (MQ=255)

                                                                                     agaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           >  2:1625174/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           <  2:1158149/100‑1 (MQ=255)

                                                                                      gaCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt          >  1:182809/1‑100 (MQ=255)

                                                                                       aCTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc         <  1:1056286/100‑1 (MQ=255)

                                                                                        cTATGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg        <  2:1879717/100‑1 (MQ=255)

                                                                                          aTGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa      >  2:1676413/1‑100 (MQ=255)

                                                                                          aTGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa      >  1:2137331/1‑100 (MQ=255)

                                                                                          aTGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa      >  1:1849458/1‑100 (MQ=255)

                                                                                           tGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc     >  2:2680526/1‑100 (MQ=255)

                                                                                           tGTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc     <  1:2175187/100‑1 (MQ=255)

                                                                                            gTTGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg    <  1:1150608/100‑1 (MQ=255)

                                                                                            gTTGAAATCC‑CCCCCTGGTATGGAAAATCACTAGCCTGATTGGCGCAGTCTTCCTGCTGTTGATTGCGATAATGGTGAGTCTGCAGGTGATTGTCGAAc     >  2:2409372/1‑99 (MQ=255)

                                                                                             ttGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt   >  1:1020337/1‑100 (MQ=255)

                                                                                              tGAAATCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg  >  1:354517/1‑100 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTG  >  NC_000913/4638088‑4638280

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