BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,033,458	T→G	22.8%	I9I (ATA→ATC)	yedR ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,033,458	0	T	G	22.8%	301.7 / 89.1	154	I9I (ATA→ATC)	yedR	inner membrane protein
Reads supporting (aligned to +/- strand): ref base T (61/57); new base G (16/19); total (77/77)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

CAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAATACA  >  NC_000913/2033362‑2033557

cAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTatgatg                                                                                                  <  2:509872/100‑1 (MQ=255)

cAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTatgatg                                                                                                  <  1:3570561/100‑1 (MQ=255)

 agggttccgccAGGTCTTCGAGAGATTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTAtgt                                                                                                 <  1:766626/94‑1 (MQ=255)

 aGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTAtgt                                                                                                 <  1:291461/100‑1 (MQ=255)

 aGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTAtgt                                                                                                 >  1:2927214/1‑100 (MQ=255)

 aGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTAtgt                                                                                                 >  2:2053990/1‑100 (MQ=255)

    ttgcgtccAGGGCTGCGAGAGAGTCTCTATTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGAtgtg                                                                                                <  2:3111830/92‑1 (MQ=255)

    ttGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGAt                                                                                              <  2:1491380/100‑1 (MQ=255)

     tGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATa                                                                                             >  2:2448449/1‑100 (MQ=255)

       ccgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATaaa                                                                                           >  1:220846/1‑100 (MQ=255)

       ccgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaaa                                                                                           >  1:3947356/1‑100 (MQ=255)

         gccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAg                                                                                         <  1:1009566/100‑1 (MQ=255)

         gccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAg                                                                                         <  2:872862/100‑1 (MQ=255)

          ccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAg                                                                                         >  1:3011179/1‑99 (MQ=255)

           cAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTc                                                                                       >  1:1997798/1‑100 (MQ=255)

           cAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTc                                                                                       <  2:3467131/100‑1 (MQ=255)

           cAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTc                                                                                       <  1:2158737/100‑1 (MQ=255)

            aGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTca                                                                                      >  2:1895842/1‑100 (MQ=255)

             gggCTGCGAGGGGGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGGCTTAGTACTATTATGTGATAAAAGTc                                                                                       <  2:1237701/98‑1 (MQ=255)

                   cGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACAttttt                                                                               >  2:2243725/1‑100 (MQ=255)

                   cGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACAttttt                                                                               <  1:2151824/100‑1 (MQ=255)

                     agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTcc                                                                             <  1:3680754/100‑1 (MQ=255)

                     agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTcc                                                                             <  1:3220899/100‑1 (MQ=255)

                       agagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCcac                                                                           <  2:1679886/100‑1 (MQ=255)

                        gagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTCGAGATAACCCGGAAAATTCTGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCcaca                                                                          <  2:2270580/100‑1 (MQ=255)

                          gTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCTACAtt                                                                        >  2:2781837/1‑100 (MQ=255)

                          gTCTCTCTTCCATCTTTGTCGAACCGTAGTCCATTTTGAGATAACCCGGAGAGTTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACAt                                                                         >  2:3536642/1‑99 (MQ=255)

                           tctctcTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTg                                                                       <  1:967356/100‑1 (MQ=255)

                            ctctctTCCATGTTTGTCGAACCGTATTCCATCTTGCGATAACCCGGCAAATTGAAAATTAGTACTATTATGTGATAAGAGACACATTTTTacac                                                                           >  2:3859009/1‑95 (MQ=255)

                             tctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGa                                                                     <  2:1905781/100‑1 (MQ=255)

                             tctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGa                                                                     <  2:2156909/100‑1 (MQ=255)

                             tctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGa                                                                     >  2:2287567/1‑100 (MQ=255)

                             tctctTCCACCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGa                                                                     >  1:3691799/1‑100 (MQ=255)

                                 ttCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGa                                                                  >  1:1996228/1‑99 (MQ=255)

                                 ttCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAg                                                                 <  2:1469278/100‑1 (MQ=255)

                                  tCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGc                                                                >  2:2123877/1‑100 (MQ=255)

                                   ccATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGcc                                                               <  2:221087/100‑1 (MQ=255)

                                    cATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCa                                                              >  1:1766229/1‑100 (MQ=255)

                                      tCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGt                                                            <  1:2547973/100‑1 (MQ=255)

                                        tttGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCa                                                          <  2:975018/100‑1 (MQ=255)

                                           gTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGa                                                       >  1:2570344/1‑100 (MQ=255)

                                            tCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGAc                                                      <  2:832625/100‑1 (MQ=255)

                                              gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                    <  2:1608806/100‑1 (MQ=255)

                                              gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                    >  1:1267636/1‑100 (MQ=255)

                                                aCCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTTACATTTTTCCACATTGGATGAGCCAGTCAAGACTCaa                                                  >  2:3898155/1‑100 (MQ=255)

                                                aCCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCaa                                                  <  2:899100/100‑1 (MQ=255)

                                                 ccGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAAt                                                 <  1:3725910/100‑1 (MQ=255)

                                                     aTTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtt                                             >  1:2018534/1‑100 (MQ=255)

                                                      ttCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttt                                            >  1:1606364/1‑100 (MQ=255)

                                                       tCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttt                                           >  1:2730010/1‑100 (MQ=255)

                                                       tCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACAGTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttt                                           >  1:188687/1‑100 (MQ=255)

                                                        ccATTTTGAGATAACATGGAAAATTCAGACTTAGTACTATTATGTGATAAAATGCACATTTTGCTGCATTGGATTAGGTAGTCAAGAGTCAATCAttttt                                          >  1:3911419/1‑100 (MQ=255)

                                                        acATTTTGAGATAACCCGGAAAATTCAGAATTAGTACTATTATGTGATAAAAGTCACATTTTTACACATTGGATGAGCCAGTCAAGAGTCAATTAtttct                                          >  2:3519552/2‑98 (MQ=255)

                                                         cATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttttt                                         <  1:53987/100‑1 (MQ=255)

                                                          aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATTTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                        >  2:1761871/1‑100 (MQ=255)

                                                          aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCGAGACTCAATCATTTTTTg                                        <  1:415224/100‑1 (MQ=255)

                                                          aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                        <  2:258115/100‑1 (MQ=255)

                                                          aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                        <  1:2448464/100‑1 (MQ=255)

                                                          aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                        >  2:281608/1‑100 (MQ=255)

                                                          aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                        <  2:2353962/100‑1 (MQ=255)

                                                          aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                        >  2:561346/1‑100 (MQ=255)

                                                           ttttGCGATAACCCGGAACATTCTGACTTATTACTATTATGTCATAGAAGTCACATTTTTACACATTCGATGAGCCAGTCAAGACTCAATCATTTTTTg                                        <  2:2567306/99‑1 (MQ=255)

                                                           ttttGAGATAACCCGTAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGa                                       <  1:31003/100‑1 (MQ=255)

                                                            tttGAGATAACTCGGAAACTTCAGATTTAGTACTATTATGTGATAAAAGGCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGa                                       <  2:768967/99‑1 (MQ=255)

                                                            tttGAGATAACCCGGAAAATTGAGACGTAGTACTAGGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                        <  2:3194768/98‑1 (MQ=255)

                                                             ttGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAAc                                     >  1:3842457/1‑100 (MQ=255)

                                                              tGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATATAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAGTCATTTTTTGAACa                                    <  2:2667974/100‑1 (MQ=255)

                                                              tGACATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACa                                    <  1:3678329/100‑1 (MQ=255)

                                                               gagaTAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACaa                                   >  1:1793238/1‑100 (MQ=255)

                                                                  aTAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATa                                <  1:209561/100‑1 (MQ=255)

                                                                  aTAAACCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATa                                <  2:771475/100‑1 (MQ=255)

                                                                     aCCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTTAATAAatatat                             >  2:1183525/1‑100 (MQ=255)

                                                                       ccGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaa                           >  2:2975198/1‑100 (MQ=255)

                                                                        cGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaaa                          >  1:689946/1‑100 (MQ=255)

                                                                          gAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAAt                         <  2:457622/99‑1 (MQ=255)

                                                                             aaTTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTg                     <  1:2037135/100‑1 (MQ=255)

                                                                              aTTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgc                    <  1:1058669/100‑1 (MQ=255)

                                                                              aTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTTTATATTTGTTTAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgc                    >  2:1097083/1‑100 (MQ=255)

                                                                              aTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgc                    >  1:2318932/1‑100 (MQ=255)

                                                                               ttCATACTTAGTACTAGTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGGCAAGACTCAATCATTTTTTGAACAAATAGATAAAttt                       <  2:2011041/96‑1 (MQ=255)

                                                                               ttCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcg                   <  1:2294167/100‑1 (MQ=255)

                                                                               ttCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcg                   <  2:2887702/100‑1 (MQ=255)

                                                                                  aGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTAAAGCCTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCaa                >  1:1008339/1‑100 (MQ=255)

                                                                                    aCTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACg              <  1:3569029/100‑1 (MQ=255)

                                                                                     cTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACgg             >  1:3431592/1‑100 (MQ=255)

                                                                                     cTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACgg             <  1:2610041/100‑1 (MQ=255)

                                                                                      ttAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGc            <  1:4070057/100‑1 (MQ=255)

                                                                                      ttAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGc            <  2:3806097/100‑1 (MQ=255)

                                                                                       tAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAACCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCACCGGCa           >  1:1415/1‑100 (MQ=255)

                                                                                         gTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGt         >  1:1915390/1‑100 (MQ=255)

                                                                                         gTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGt         >  1:268781/1‑100 (MQ=255)

                                                                                          tACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGTACAAATAGATAAATTTTGCGCAACGGCAGTg        >  2:1954536/1‑100 (MQ=255)

                                                                                           acaattatGTGATAAAAGCCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTg        >  2:1754897/4‑99 (MQ=255)

                                                                                           aCTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGa       <  2:1846191/100‑1 (MQ=255)

                                                                                           aCTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGa       >  1:2006988/1‑100 (MQ=255)

                                                                                            cTATAATGTGATGAAAGTCGTGTTTGTCCAAATTGAATGAGCCCGTCATGACTCAATCATTTTTTGCACAAATGGATAAATTTTgcggacgggt            <  2:1046648/94‑8 (MQ=255)

                                                                                              attatGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGTCAGTGAATa    >  1:3749176/1‑100 (MQ=255)

                                                                                              attatGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAATa    >  1:2164083/1‑100 (MQ=255)

                                                                                               ttatGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAATAc   <  2:3687822/100‑1 (MQ=255)

                                                                                                tatGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAATACa  >  1:2292112/1‑100 (MQ=255)

CAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAATACA  >  NC_000913/2033362‑2033557

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