BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,033,458	T→G	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	98.2%	396.7 / ‑4.1	109	I9I (ATA→ATC)	yedR	inner membrane protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base G (57/50); total (57/52)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

AGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAA  >  NC_000913/2033363‑2033553

gatctGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGAtgt                                                                                             <  1:1898815/96‑1 (MQ=255)

aGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGAtgt                                                                                             <  1:1116860/100‑1 (MQ=255)

 gTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGAtgtg                                                                                            >  1:2143147/1‑100 (MQ=255)

 gTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGAtgtg                                                                                            <  1:526563/100‑1 (MQ=255)

 gTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGAtgtg                                                                                            <  2:369424/100‑1 (MQ=255)

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

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

    tGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATa                                                                                         >  1:618527/1‑100 (MQ=255)

     gccgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGAtgtagatc                                                                                        >  2:1898808/1‑95 (MQ=255)

       cgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaaaa                                                                                      >  1:1065362/1‑100 (MQ=255)

       cgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaaaa                                                                                      >  1:1608936/1‑100 (MQ=255)

       cgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaaaa                                                                                      <  2:2316915/100‑1 (MQ=255)

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

            gggCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTcac                                                                                 <  2:2256879/100‑1 (MQ=255)

            gggCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTcac                                                                                 <  2:2260996/100‑1 (MQ=255)

            gggCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTcac                                                                                 <  2:784987/100‑1 (MQ=255)

             ggCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTcaca                                                                                >  1:1920188/1‑100 (MQ=255)

               cTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAtt                                                                              >  1:1784371/1‑100 (MQ=255)

                  cGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAttttt                                                                           >  2:1094123/1‑100 (MQ=255)

                   gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                          <  2:617122/100‑1 (MQ=255)

                   gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                          <  1:330454/100‑1 (MQ=255)

                    agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTcc                                                                         <  2:123887/100‑1 (MQ=255)

                    agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTcc                                                                         <  2:1733330/100‑1 (MQ=255)

                    agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                          <  2:1909471/99‑1 (MQ=255)

                      agagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCcac                                                                       <  2:1418537/100‑1 (MQ=255)

                        agTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACAt                                                                     <  1:1094125/100‑1 (MQ=255)

                           ctctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTgg                                                                  <  1:1846523/100‑1 (MQ=255)

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

                             ctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGAt                                                                >  1:1665587/1‑100 (MQ=255)

                             ctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGAt                                                                <  2:1341144/100‑1 (MQ=255)

                              tctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATg                                                               <  1:147375/100‑1 (MQ=255)

                                ttCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAg                                                             >  2:966283/1‑100 (MQ=255)

                                ctCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAg                                                             <  1:1397884/99‑1 (MQ=255)

                                 tCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGc                                                            >  1:384601/1‑100 (MQ=255)

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

                                  ccATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAgag                                                           >  2:1397881/1‑98 (MQ=255)

                                   cATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCa                                                          >  2:387011/1‑100 (MQ=255)

                                      cTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTc                                                       <  2:618526/100‑1 (MQ=255)

                                      cTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTc                                                       <  2:601988/100‑1 (MQ=255)

                                      cTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTc                                                       <  2:2284324/100‑1 (MQ=255)

                                      cTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTAGATGAGCCAGTc                                                       >  2:623885/1‑100 (MQ=255)

                                       tttGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCa                                                      <  1:1454269/100‑1 (MQ=255)

                                         tGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAg                                                    <  2:1617676/100‑1 (MQ=255)

                                         tGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAg                                                    <  2:446629/100‑1 (MQ=255)

                                         tGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAg                                                    <  1:1686719/100‑1 (MQ=255)

                                         tGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAg                                                    <  1:169753/100‑1 (MQ=255)

                                          gTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGa                                                   >  1:1427112/1‑100 (MQ=255)

                                          gTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGa                                                   >  1:402982/1‑100 (MQ=255)

                                            cGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACt                                                 <  2:402980/100‑1 (MQ=255)

                                            cGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACt                                                 <  2:1261703/100‑1 (MQ=255)

                                            cGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACt                                                 >  1:2147286/1‑100 (MQ=255)

                                             gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                >  2:405303/1‑100 (MQ=255)

                                             gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                >  1:1441275/1‑100 (MQ=255)

                                             gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                >  1:427321/1‑100 (MQ=255)

                                                ccGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAAt                                             <  1:1875349/100‑1 (MQ=255)

                                                ccGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAAt                                             <  1:1359563/100‑1 (MQ=255)

                                                 cGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGACAGTCAAGACTCAATc                                            >  1:324435/1‑100 (MQ=255)

                                                  gTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCa                                           >  2:39028/1‑100 (MQ=255)

                                                    aTTCCATTTTGAGATAACCCGGAAAACTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtt                                         >  2:281917/1‑100 (MQ=255)

                                                     ttCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttt                                        >  2:205104/1‑100 (MQ=255)

                                                      tCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttt                                       >  1:2015452/1‑100 (MQ=255)

                                                       ccATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttttt                                      >  1:1311930/1‑100 (MQ=255)

                                                       ccATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttttt                                      >  1:2029114/1‑100 (MQ=255)

                                                       ccATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttttt                                      >  1:210990/1‑100 (MQ=255)

                                                        cATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttttt                                     >  2:546978/1‑100 (MQ=255)

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

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

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

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

                                                          ttttGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGa                                   <  1:1191998/100‑1 (MQ=255)

                                                            ttGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAAc                                 <  1:1037177/100‑1 (MQ=255)

                                                             tGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACa                                <  1:1386940/100‑1 (MQ=255)

                                                             tGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACa                                >  2:1974491/1‑100 (MQ=255)

                                                             tGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACa                                >  1:2081622/1‑100 (MQ=255)

                                                             tGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACa                                <  1:1630343/100‑1 (MQ=255)

                                                               agaTAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACaaa                              >  2:1737822/1‑100 (MQ=255)

                                                               acaTAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACaaa                              >  1:1580247/3‑100 (MQ=255)

                                                                  tAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAg                           <  2:150818/100‑1 (MQ=255)

                                                                   aaCCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGa                          >  1:1073797/1‑100 (MQ=255)

                                                                      ccGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaa                       >  1:572345/1‑100 (MQ=255)

                                                                      ccGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaa                       >  1:798199/1‑100 (MQ=255)

                                                                        ggAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAAt                     <  1:1798169/100‑1 (MQ=255)

                                                                        ggAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAAt                     <  2:1065360/100‑1 (MQ=255)

                                                                            aaTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTg                 >  2:1673624/1‑100 (MQ=255)

                                                                            aaTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTg                 <  2:2015445/100‑1 (MQ=255)

                                                                             ccgatctcTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgc                <  2:1155866/93‑1 (MQ=255)

                                                                             aTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgc                <  1:205106/100‑1 (MQ=255)

                                                                               tCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcgc              >  1:1083583/1‑100 (MQ=255)

                                                                                  gACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAAc           >  1:114702/1‑100 (MQ=255)

                                                                                    cTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcagatcgg         >  1:1155868/1‑93 (MQ=255)

                                                                                    cTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACgg         >  1:2046570/1‑100 (MQ=255)

                                                                                       aGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGTACATATAGATAAATTTTGCGCAACGGCAg      >  2:109400/1‑100 (MQ=255)

                                                                                       aGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAg      >  2:856484/1‑100 (MQ=255)

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

                                                                                          aCTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGa   <  2:1580242/100‑1 (MQ=255)

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

                                                                                           ntatgatGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGaa  >  1:1132457/2‑100 (MQ=255)

                                                                                           cTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGaa  >  2:1209110/1‑100 (MQ=255)

                                                                                           cTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGaa  <  1:2092629/100‑1 (MQ=255)

                                                                                           cTATGATGTAATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGaa  >  2:1440120/1‑100 (MQ=255)

AGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAA  >  NC_000913/2033363‑2033553

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