BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,638,193:1	(C)_6→7	100%	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	93.0%	510.0 / 28.9	158	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (5/6); new base C (91/56); total (96/62)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGAT  >  NC_000913/4638088‑4638285

   tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc                                                                                                  <  2:226762/100‑1 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑ccc                                                                                              >  1:4006018/1‑100 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccc                                                                                             <  1:684987/100‑2 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccc                                                                                             <  1:1223835/100‑2 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccc                                                                                             <  2:390976/100‑2 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccc                                                                                             <  2:32484/100‑2 (MQ=255)

         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccctg                                                                                           <  1:3326219/100‑4 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccctgt                                                                                          >  2:2056206/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccctgt                                                                                          >  1:4188955/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTACCATTGCAAAGCAGAAGACTATATTGAAA‑T‑ccc‑ccc                                                                                              >  1:1597976/1‑96 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccctgtt                                                                                         <  1:2309359/100‑6 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCTAATTCTTCCCACATAGTCTTCGTATCCTGCTTCTATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                         <  1:105569/100‑1 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTTCAAATGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                         >  1:490379/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                         <  2:2395258/100‑1 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                         >  2:495658/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                         >  2:1635490/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                         >  2:848287/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                         >  2:951145/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTg                                                                                       <  1:3918982/100‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTgg                                                                                      >  1:758784/1‑100 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGa                                                                                     >  2:267192/1‑100 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAAt                                                                                 >  2:1921302/1‑100 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATg                                                                                <  1:4250775/100‑1 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGCAACAGCCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATg                                                                                >  1:2166727/1‑99 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGAc                                                                              >  1:1025845/1‑100 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGAc                                                                              <  2:2202232/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGAc                                                                              <  2:4380173/100‑1 (MQ=255)

                       cTTCACATAGCTTAAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGAc                                                                              >  1:2730635/1‑100 (MQ=255)

                        ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGTAGAAGACTATTTTGAAA‑TCCCC‑CCCTGTGCTGGAAAATGatt                                                                             <  2:2974028/100‑3 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAg                                                                           <  2:26867/100‑1 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACAGAGTCTTCGTATCCTGCTGCCATTGCAAAGGTGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAg                                                                           <  1:4315381/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTTCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGcc                                                                         >  2:1178651/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGttt                                                                    >  2:3792291/1‑100 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGttt                                                                    >  2:4256658/1‑100 (MQ=255)

                                  ttCAATTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAACGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGTAAAATGACTAGCCTGttt                                                                    >  1:4090653/1‑99 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                                  >  2:360146/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                                  <  1:3993405/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                               >  2:314333/1‑100 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                             <  1:2092895/100‑1 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                        >  2:1792703/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                        >  2:1468857/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                       >  1:3518311/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                       <  2:3041348/99‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                      >  1:2963151/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                     >  1:2181067/1‑100 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATAGCAAAGGAGAAGACTATGTTGATA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                       <  1:320052/99‑1 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                     <  2:3485404/100‑1 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                                  >  2:61243/1‑100 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCCTCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                                  <  1:4202007/100‑1 (MQ=255)

                                                     ntCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTACTgttgtt                                                >  1:900288/2‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                                >  1:1713702/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                                >  1:3057369/1‑100 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                               >  2:1064627/1‑100 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                              <  1:1186407/100‑1 (MQ=255)

                                                        ttCGTTTCCTGCTGCCATTGCAAAAGAGAAGTCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                             <  1:4125154/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                          >  1:3214537/1‑100 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                         >  1:2470048/1‑100 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                       <  1:3714232/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                       >  2:1649729/1‑100 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                      >  2:4156230/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                     >  2:217803/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                     >  1:1767377/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                      <  1:3356306/99‑1 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                    >  1:1751456/1‑100 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                    <  1:846660/100‑1 (MQ=255)

                                                                 tgctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                    <  1:3753317/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTTTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                                  <  1:2348540/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                                 <  1:1544548/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                                 <  1:1275755/100‑1 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                                >  1:2750330/1‑100 (MQ=255)

                                                                     gcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                                >  2:684841/1‑100 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                               >  2:2867212/1‑100 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                              <  1:3107060/100‑1 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                              >  1:1825807/1‑100 (MQ=255)

                                                                         ttGCAAAGGAGAAGATTATGTTGAAATT‑CTC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                            <  2:1279506/100‑1 (MQ=255)

                                                                         ttGCAAAGGAGAAGACTATGTTGAAAGC‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                            <  1:2709831/100‑1 (MQ=255)

                                                                           gCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                          >  1:3756787/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                         <  2:298394/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                         >  1:3811154/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                         >  1:3204859/1‑100 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                       <  1:177323/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAGA‑TCCCC‑CCCTGTTCGGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTGATTGATTCCGATAATGCTGATTCGGCAg                       >  1:883730/1‑98 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                     <  2:3884696/100‑1 (MQ=255)

                                                                                 gagaAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                    <  2:3456858/100‑1 (MQ=255)

                                                                                   gaagaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt                  >  2:3560299/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg                <  2:2044075/100‑1 (MQ=255)

                                                                                      gaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt               >  2:329939/1‑100 (MQ=255)

                                                                                      gaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt               >  1:2777828/1‑100 (MQ=255)

                                                                                       aCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc              >  2:719097/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg             >  1:3452366/1‑100 (MQ=255)

                                                                                        cTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg             <  1:719096/100‑1 (MQ=255)

                                                                                         tATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa            <  2:1203271/100‑1 (MQ=255)

                                                                                           tGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc          >  2:2321355/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc          >  2:246413/1‑100 (MQ=255)

                                                                                           tGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc          <  2:3040107/100‑1 (MQ=255)

                                                                                            gTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg         >  2:362523/1‑100 (MQ=255)

                                                                                              tGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg       >  2:1266542/1‑100 (MQ=255)

                                                                                              tGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg       >  1:2486549/1‑100 (MQ=255)

                                                                                              tGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg       >  2:1881062/1‑100 (MQ=255)

                                                                                              tGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg       >  1:385571/1‑100 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGAT  >  NC_000913/4638088‑4638285

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