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	97.2%	376.7 / 3.3	108	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (2/1); new base C (47/58); total (49/59)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 119 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                 >  1:2343213/1‑100 (MQ=255)

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

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc                                                                                           >  1:811438/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccc                                                                                           >  1:1670266/1‑100 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccccc                                                                                          >  2:1328813/1‑100 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cccccc                                                                                          >  1:738346/1‑100 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccn                                                                                         <  2:2945469/100‑2 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccc                                                                                         <  1:1588458/100‑2 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccctgt                                                                                      >  2:228077/1‑96 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑ccccccctgtt                                                                                     >  1:2571757/1‑95 (MQ=255)

             ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCt                                                                                    <  1:898988/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTg                                                                                   <  2:2402706/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTg                                                                                   >  2:1872548/1‑100 (MQ=255)

               cgccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTgg                                                                                  <  2:2849203/98‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCGCACATAGTCTGCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTg                                                                                   >  1:2083663/1‑99 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTgg                                                                                  >  2:1031563/1‑100 (MQ=255)

                aCCGTCACTTTACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGTGATGACTATGTTGAAATCCCCCCCTGTTCTGGa                                                                                 <  2:2031562/100‑1 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGa                                                                                 >  1:2967664/1‑100 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGa                                                                                 <  1:2010143/100‑1 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGaa                                                                                >  2:1920369/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGaa                                                                                >  1:2217750/1‑100 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGaaaan                                                                             <  2:2826615/100‑2 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAAt                                                                             <  2:1371765/100‑1 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAAt                                                                             <  2:1288127/100‑1 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATg                                                                            <  1:2140143/100‑1 (MQ=255)

                     cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCAGTGCAAAGGAAAAGACGATGTTGAAGTCCCCCCCTGTTCTGGAAAATg                                                                            >  2:2081674/1‑100 (MQ=255)

                      aCTTCACATAGCTTCAAATTCTTCCCACATAGGCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAAt                                                                             <  2:1119321/98‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGAc                                                                          <  1:2987416/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGAc                                                                          >  2:2987386/1‑100 (MQ=255)

                        ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACt                                                                         <  2:1772477/100‑1 (MQ=255)

                        ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACt                                                                         <  2:2217551/100‑1 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTa                                                                        >  1:1058184/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                     >  1:2007846/1‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:661568/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:658899/100‑1 (MQ=255)

                             aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCt                                                                    <  1:1830797/100‑1 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:336353/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:2549827/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:1990758/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:497594/1‑100 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:2416096/1‑100 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  2:333138/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:1673623/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           >  2:2793077/1‑100 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:1872569/100‑1 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          >  1:1288509/1‑100 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          >  2:1313612/1‑100 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  2:1246508/100‑1 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:2662913/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:2071351/100‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  1:3684/1‑100 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  2:290396/1‑100 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:398132/100‑1 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:1404271/100‑1 (MQ=255)

                                                 caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  1:2367040/100‑1 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  2:2736833/100‑1 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             >  2:2956245/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:269732/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  1:258173/1‑100 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  1:1996953/1‑100 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  2:1522631/100‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  1:1920390/100‑1 (MQ=255)

                                                         tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAtt                                        >  2:1787609/1‑100 (MQ=255)

                                                         tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAAGGACTAGCCTGTGTGGTGCAGTATTGCTGTTGTTGAtt                                        >  1:1683491/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:2710287/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:2482595/100‑1 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  2:2932453/1‑100 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  1:1817566/1‑100 (MQ=255)

                                                            tATCCTGCTGCCATTACAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     <  1:1614986/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  2:1460268/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:770696/100‑1 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  <  1:1031574/100‑1 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  2:2923830/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               >  1:466282/1‑100 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:1375999/100‑1 (MQ=255)

                                                                       cATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAACATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           <  2:2578069/99‑1 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         <  1:2043923/100‑1 (MQ=255)

                                                                         ttGCAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        <  2:2009305/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:790309/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  2:2007824/100‑1 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   <  2:1817546/100‑1 (MQ=255)

                                                                               aGGAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  <  2:1990736/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:2464757/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  2:1058173/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  1:1442312/1‑100 (MQ=255)

                                                                                 gagaAGACTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                <  2:2514435/100‑1 (MQ=255)

                                                                                   gaagaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  2:364922/100‑1 (MQ=255)

                                                                                   gaagaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  1:3007414/100‑1 (MQ=255)

                                                                                   gaagaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTAGGCAGGTGAt              >  2:2023434/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  2:1894846/1‑100 (MQ=255)

                                                                                      gaCTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  2:1229500/1‑100 (MQ=255)

                                                                                        cTATGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:2870650/1‑100 (MQ=255)

                                                                                          aTGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  2:811433/100‑1 (MQ=255)

                                                                                           tGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  1:2303622/1‑100 (MQ=255)

                                                                                           tGTTGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCGGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:1653074/1‑100 (MQ=255)

                                                                                             ttGGAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTGTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  2:1649805/1‑100 (MQ=255)

                                                                                             ttGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:1289601/1‑100 (MQ=255)

                                                                                             ttGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:352747/1‑100 (MQ=255)

                                                                                              tGAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  1:1891928/1‑100 (MQ=255)

                                                                                               gAAATCCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  2:2115979/100‑1 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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: `‑`