BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	945,097	T→C	S2P (TCC→CCC)	ycaD →	putative transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	945,097	0	T	C	99.2%	397.1 / ‑4.7	127	S2P (TCC→CCC)	ycaD	putative transport
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (68/58); total (69/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

Only 100 of 127 total aligned reads displayed.

GGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTT  >  NC_000913/945008‑945186

ggATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCa                                                                                        <  1:644004/93‑1 (MQ=255)

 gATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCAc                                                                                       <  1:1752049/93‑1 (MQ=255)

   tatCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGt                                                                                     >  2:3879100/1‑93 (MQ=255)

    atCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATTCCCACGta                                                                                    <  2:105057/93‑1 (MQ=255)

    atCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGta                                                                                    >  2:1665312/1‑93 (MQ=255)

     tCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGtat                                                                                   >  1:1451288/1‑93 (MQ=255)

      cgcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGtata                                                                                  <  1:2104318/93‑1 (MQ=255)

       gcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATAc                                                                                 >  1:3668097/1‑93 (MQ=255)

       gcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATAc                                                                                 >  1:1446811/1‑93 (MQ=255)

       gcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATAc                                                                                 >  1:4061524/1‑93 (MQ=255)

       gcgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATAc                                                                                 <  2:2346231/93‑1 (MQ=255)

        cgTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATAcc                                                                                <  1:2067806/93‑1 (MQ=255)

         gTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATAccc                                                                               >  2:414319/1‑93 (MQ=255)

         gTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATAccc                                                                               >  1:128663/1‑93 (MQ=255)

           aTTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAg                                                                             >  2:2559521/1‑93 (MQ=255)

              ttCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCt                                                                          <  1:2210576/93‑1 (MQ=255)

              ttCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCt                                                                          <  2:429206/93‑1 (MQ=255)

                   ccATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCAt                                                                     >  2:143950/1‑93 (MQ=255)

                     aTTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGc                                                                   <  2:2474518/93‑1 (MQ=255)

                     aTTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGc                                                                   >  1:398059/1‑93 (MQ=255)

                      ttGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCt                                                                  <  1:2326848/93‑1 (MQ=255)

                      ttGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCt                                                                  >  2:2166829/1‑93 (MQ=255)

                        gACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCttt                                                                <  2:3128956/93‑1 (MQ=255)

                        gACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCttt                                                                >  1:1650690/1‑93 (MQ=255)

                         aCAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtttt                                                               >  1:2546846/1‑93 (MQ=255)

                         aCAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtttt                                                               >  2:3604026/1‑93 (MQ=255)

                         aCAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtttt                                                               <  1:965592/93‑1 (MQ=255)

                         aCAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtttt                                                               >  2:1677526/1‑93 (MQ=255)

                         aCAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCtttt                                                               <  2:630229/93‑1 (MQ=255)

                           aaTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGc                                                             >  2:4193076/1‑93 (MQ=255)

                           aaTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGc                                                             >  2:790002/1‑93 (MQ=255)

                            aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCt                                                            >  1:4143704/1‑93 (MQ=255)

                            aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCt                                                            >  1:2941134/1‑93 (MQ=255)

                            aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCt                                                            <  2:2835085/93‑1 (MQ=255)

                            aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCt                                                            >  1:2528891/1‑93 (MQ=255)

                            aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCt                                                            <  2:384949/93‑1 (MQ=255)

                            aTGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCt                                                            <  2:4653108/93‑1 (MQ=255)

                              gTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGt                                                          >  2:3183860/1‑93 (MQ=255)

                               tttttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATTCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTc                                                         <  2:93622/93‑1 (MQ=255)

                                ttttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCt                                                        <  2:2743709/93‑1 (MQ=255)

                                ttttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCt                                                        <  2:38140/93‑1 (MQ=255)

                                ttttGGCGGTGGAATGATGCGCATGAAATTGGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCt                                                        <  2:1510641/93‑1 (MQ=255)

                                 tttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTg                                                       <  1:319234/93‑1 (MQ=255)

                                 tttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTg                                                       <  1:769419/93‑1 (MQ=255)

                                  ttGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTgg                                                      <  1:3715427/93‑1 (MQ=255)

                                   tGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGc                                                     >  2:4484474/1‑93 (MQ=255)

                                     gCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCt                                                   <  2:1063476/93‑1 (MQ=255)

                                       ggTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGc                                                 <  1:408793/93‑1 (MQ=255)

                                          ggCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCttt                                              <  2:604134/93‑1 (MQ=255)

                                          ggCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCttt                                              >  2:3023992/1‑93 (MQ=255)

                                           gCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCtttt                                             <  1:4628692/93‑1 (MQ=255)

                                           gCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCtttt                                             <  2:4235869/93‑1 (MQ=255)

                                             atgatgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgt                                           >  1:4209856/1‑93 (MQ=255)

                                               gatgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTttgttg                                         >  1:4062495/1‑93 (MQ=255)

                                                atgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGa                                        >  2:3734235/1‑93 (MQ=255)

                                                atgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGa                                        >  2:1748494/1‑93 (MQ=255)

                                                atgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGa                                        >  1:599417/1‑93 (MQ=255)

                                                 tgCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGAc                                       >  2:3413842/1‑93 (MQ=255)

                                                   cgcATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGActc                                     <  1:384825/93‑1 (MQ=255)

                                                   cgcATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGActc                                     <  1:66352/93‑1 (MQ=255)

                                                    gcATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGActct                                    >  1:4269383/1‑93 (MQ=255)

                                                     cATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTg                                   <  2:2059345/93‑1 (MQ=255)

                                                     cATGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTg                                   <  1:2723482/93‑1 (MQ=255)

                                                      aTGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTgg                                  >  2:3568360/1‑93 (MQ=255)

                                                      aTGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTgg                                  <  2:1795395/93‑1 (MQ=255)

                                                       tGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGc                                 >  2:2365765/1‑93 (MQ=255)

                                                       tGAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGc                                 <  2:3878787/93‑1 (MQ=255)

                                                        gAAATTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCg                                >  2:1433668/1‑93 (MQ=255)

                                                         aaaTTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGa                               <  1:3842571/93‑1 (MQ=255)

                                                         aaaTTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGa                               >  1:2807335/1‑93 (MQ=255)

                                                         aaaTTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGa                               >  1:2347106/1‑93 (MQ=255)

                                                          aaTTTGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGAt                              >  2:442992/1‑93 (MQ=255)

                                                            tttGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTg                            >  1:1976539/1‑93 (MQ=255)

                                                             ttGAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGc                           <  1:3387479/93‑1 (MQ=255)

                                                               gAACTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCgg                         <  1:2722244/93‑1 (MQ=255)

                                                                aaCTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgt                        >  1:1857265/1‑93 (MQ=255)

                                                                  cTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgtgt                      >  1:2376645/1‑93 (MQ=255)

                                                                  cTTCCTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGgtgt                      <  2:2039863/93‑1 (MQ=255)

                                                                     ccTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaa                   >  1:4613904/1‑93 (MQ=255)

                                                                      cTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTTGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaaa                  <  1:4637310/93‑1 (MQ=255)

                                                                      cTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaaa                  <  1:228837/93‑1 (MQ=255)

                                                                      cTCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTaaa                  >  2:3664737/1‑93 (MQ=255)

                                                                       tCACGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAAt                 <  1:316622/93‑1 (MQ=255)

                                                                         agggTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGCTTGCGGTGTTAAATac               >  2:520142/3‑93 (MQ=255)

                                                                         aCGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATac               <  2:3433081/93‑1 (MQ=255)

                                                                          cGGTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATaca              >  1:3021672/1‑93 (MQ=255)

                                                                            gTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACt            <  2:3271146/93‑1 (MQ=255)

                                                                            gTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACt            >  2:2588771/1‑93 (MQ=255)

                                                                            gTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACt            <  1:763735/93‑1 (MQ=255)

                                                                            gTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACt            >  1:1618811/1‑93 (MQ=255)

                                                                            gTTTTAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACt            >  2:2487576/1‑93 (MQ=255)

                                                                             ttttAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTc           >  1:112690/1‑93 (MQ=255)

                                                                              tttAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCg          >  1:1088584/1‑93 (MQ=255)

                                                                              tttAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCCGGCGATTGCGGTGTCAAATACACTCg          >  1:1588045/1‑93 (MQ=255)

                                                                               ttAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGt         >  1:2838205/1‑93 (MQ=255)

                                                                               ttAATTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGt         >  1:119275/1‑93 (MQ=255)

                                                                                  aTTCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGcc      >  2:816250/1‑93 (MQ=255)

                                                                                   ttCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCg     >  2:4377712/1‑93 (MQ=255)

                                                                                   ttCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCg     >  2:3064434/1‑93 (MQ=255)

                                                                                   ttCATGCCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCg     >  2:1440761/1‑93 (MQ=255)

GGATATCGCGTATTTTCACCCATTGACAATGTTTTTGGCGGTGGCATGATGCGCATGAAATTTGAACTTCCTCACGGTTTTAATTCATGTCCACGTATACCCAGCCTGTCATGCTTTTGCTGTCTGGCCTGCTTTTGTTGACTCTGGCGATTGCGGTGTTAAATACACTCGTGCCGCTT  >  NC_000913/945008‑945186

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