BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	263,172	A→C	pseudogene (41/63 nt)	ykfN ←	pseudogene, YdiA family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	263,172	0	A	C	98.8%	591.1 / ‑2.7	168	pseudogene (41/63 nt)	ykfN	pseudogene, YdiA family
Reads supporting (aligned to +/- strand): ref base A (0/2); new base C (88/78); total (88/80)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 168 total aligned reads displayed.

GGCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAC  >  NC_000913/263075‑263267

ggCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGc                                                                                               <  1:3626313/100‑1 (MQ=255)

 gCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCt                                                                                              >  2:1490671/1‑100 (MQ=255)

   ttCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCa                                                                                            <  2:2417427/100‑1 (MQ=255)

   ttCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCa                                                                                            <  2:1541065/100‑1 (MQ=255)

    tCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCAt                                                                                           >  1:1251048/1‑100 (MQ=255)

     cGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGcagatc                                                                                          >  1:1670727/1‑95 (MQ=255)

      gCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTACTGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCAtt                                                                                          >  1:978705/1‑99 (MQ=255)

       cTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTAt                                                                                        <  2:183641/100‑1 (MQ=255)

         cACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCAGTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGa                                                                                      >  2:2389777/1‑100 (MQ=255)

         cACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGa                                                                                      >  2:1073242/1‑100 (MQ=255)

          ncTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGAt                                                                                     >  2:352139/2‑100 (MQ=255)

           cTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATc                                                                                    <  1:108790/100‑1 (MQ=255)

            tataCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCg                                                                                   >  2:1753758/1‑100 (MQ=255)

              tataTCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGt                                                                                 <  2:3475637/97‑1 (MQ=255)

              taCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGt                                                                                 <  1:2737783/100‑1 (MQ=255)

               aCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGtt                                                                                >  2:1434477/1‑100 (MQ=255)

                cATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTa                                                                               <  2:1615519/100‑1 (MQ=255)

                cATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTa                                                                               >  2:3018107/1‑100 (MQ=255)

                cATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTa                                                                               <  2:725750/100‑1 (MQ=255)

                    cTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATAc                                                                           >  2:1868833/1‑100 (MQ=255)

                    cTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATAc                                                                           >  1:3174677/1‑100 (MQ=255)

                     ttGGCTGCAGGTTTAGTTGTACACCACTCGTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATAcn                                                                          <  2:2332115/100‑2 (MQ=255)

                     ttGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACt                                                                          >  2:7976/1‑100 (MQ=255)

                       ggCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTc                                                                        <  1:1434484/100‑1 (MQ=255)

                       ggCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTc                                                                        >  2:59040/1‑100 (MQ=255)

                       ggCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTc                                                                        >  2:3653164/1‑100 (MQ=255)

                       ggCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTc                                                                        >  2:494945/1‑100 (MQ=255)

                        gCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCa                                                                       >  1:2330201/1‑100 (MQ=255)

                        gCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCa                                                                       >  1:652741/1‑100 (MQ=255)

                         cTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCaa                                                                      >  1:2293682/1‑100 (MQ=255)

                              ggTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCt                                                                 <  2:2823232/100‑1 (MQ=255)

                              ggTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCt                                                                 >  2:3599745/1‑100 (MQ=255)

                               gTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTg                                                                <  2:2246653/100‑1 (MQ=255)

                               gTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTg                                                                >  1:870166/1‑100 (MQ=255)

                                    gTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGc                                                           >  2:3445914/1‑100 (MQ=255)

                                     ttGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCa                                                          <  1:3111547/100‑1 (MQ=255)

                                        tACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGa                                                       >  2:3231049/1‑100 (MQ=255)

                                         acacCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAt                                                      <  2:2965000/100‑1 (MQ=255)

                                           accacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttg                                                    >  1:3238695/1‑100 (MQ=255)

                                           accacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttg                                                    <  1:2875355/100‑1 (MQ=255)

                                           accacTCCTAAATTTAATGTGTTGGCAATATGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttg                                                    <  2:2327283/100‑1 (MQ=255)

                                             atcTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtt                                                  <  1:2825862/98‑1 (MQ=255)

                                             atcTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtt                                                  <  1:513665/98‑1 (MQ=255)

                                              acTCCTAAATTTAATGTGTTGGCAATGTGTTCAATATAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgttt                                                 <  1:2721756/100‑1 (MQ=255)

                                               cTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtttg                                                <  1:7591/100‑1 (MQ=255)

                                               cTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtttg                                                >  2:2559429/1‑100 (MQ=255)

                                                 ccTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTc                                              <  1:3025446/100‑1 (MQ=255)

                                                 ccTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTc                                              <  2:1685523/100‑1 (MQ=255)

                                                   tAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCg                                            <  1:248863/100‑1 (MQ=255)

                                                    aaaTTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGt                                           <  2:2435658/100‑1 (MQ=255)

                                                    aaaTTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGt                                           >  2:171699/1‑100 (MQ=255)

                                                     aaTTTAATGTGTTAGCAATGTGTTCAATAAAGCTCGAACAAAGTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTa                                          >  1:2623467/1‑100 (MQ=255)

                                                       tttAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTaaa                                        >  2:3604351/1‑100 (MQ=255)

                                                          aaTGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAg                                     <  2:3089098/100‑1 (MQ=255)

                                                          aaTGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAg                                     <  2:3044707/100‑1 (MQ=255)

                                                            tgtgTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGAt                                   >  1:1377897/1‑100 (MQ=255)

                                                            tgtgTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGAt                                   <  1:2939614/100‑1 (MQ=255)

                                                            tgtgTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGAAAAAAg                                     <  2:1999100/98‑1 (MQ=255)

                                                             gtgtTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATa                                  >  2:3425326/1‑100 (MQ=255)

                                                              tgtTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGAtcg                                 >  1:3779468/1‑98 (MQ=255)

                                                              tgtTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATaa                                 >  2:3535883/1‑100 (MQ=255)

                                                                 tGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcgc                              >  2:2702731/1‑100 (MQ=255)

                                                                 tGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcgc                              >  1:3486456/1‑100 (MQ=255)

                                                                 tGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCAAGATTGTTTGTCCGTAAAAAGATAAcgc                              <  2:3204004/100‑1 (MQ=255)

                                                                  ggCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcgcg                             >  2:2592019/1‑100 (MQ=255)

                                                                     aaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCt                          >  2:1192180/1‑100 (MQ=255)

                                                                     aaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCt                          >  2:3742375/1‑100 (MQ=255)

                                                                      aTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTg                         <  1:1490678/100‑1 (MQ=255)

                                                                      aTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTg                         >  2:1571753/1‑100 (MQ=255)

                                                                        gtgtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGcc                       >  1:54781/1‑100 (MQ=255)

                                                                         tgtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGcta                      >  1:1401135/1‑98 (MQ=255)

                                                                         tgtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCg                      <  2:1395445/100‑1 (MQ=255)

                                                                          gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTGATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGcc                       >  2:2130123/1‑98 (MQ=255)

                                                                          gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     >  2:619785/1‑100 (MQ=255)

                                                                          gtTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCgg                     <  1:419683/100‑1 (MQ=255)

                                                                           ttCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCggg                    <  2:92524/100‑1 (MQ=255)

                                                                           ttCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCggg                    <  2:2893113/100‑1 (MQ=255)

                                                                           ttCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCggg                    <  1:3742408/100‑1 (MQ=255)

                                                                           ttCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTACCggg                    <  1:2164578/100‑1 (MQ=255)

                                                                            tCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgt                   <  2:588468/100‑1 (MQ=255)

                                                                            tCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgt                   >  2:2002229/1‑100 (MQ=255)

                                                                             cAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgta                  >  1:613624/1‑100 (MQ=255)

                                                                             cAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgta                  >  1:783206/1‑100 (MQ=255)

                                                                               aTAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgtagt                >  2:673319/1‑100 (MQ=255)

                                                                                tAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgtagta               <  2:1078891/100‑1 (MQ=255)

                                                                                tAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgtagta               <  1:45120/100‑1 (MQ=255)

                                                                                tAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAGTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgtagta               <  2:3572511/100‑1 (MQ=255)

                                                                                tAAAGCTCAAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgtagta               <  2:1271799/100‑1 (MQ=255)

                                                                                  aaGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGc             <  1:1868842/100‑1 (MQ=255)

                                                                                  aaGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGc             <  1:1867506/100‑1 (MQ=255)

                                                                                  aaGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGc             <  2:2588355/100‑1 (MQ=255)

                                                                                    gCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAg           <  2:3771024/100‑1 (MQ=255)

                                                                                       cGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgcg        >  1:3506099/1‑100 (MQ=255)

                                                                                       cGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgcg        <  2:3445663/100‑1 (MQ=255)

                                                                                         aaCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCa      <  2:1817713/100‑1 (MQ=255)

                                                                                         aaCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCa      <  1:1091645/100‑1 (MQ=255)

                                                                                          aCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCAt     <  1:177507/100‑1 (MQ=255)

                                                                                            aaaTTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTa   <  1:2392493/100‑1 (MQ=255)

                                                                                             aaTTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAc  >  1:2556644/1‑100 (MQ=255)

                                                                                             aaTTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAc  <  1:3138707/100‑1 (MQ=255)

GGCTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAC  >  NC_000913/263075‑263267

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