BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	712,150	A→T	L218Q (CTG→CAG)	ybfF ←	acyl‑CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	712,150	0	A	T	98.3%	392.4 / ‑2.7	117	L218Q (CTG→CAG)	ybfF	acyl‑CoA esterase
Reads supporting (aligned to +/- strand): ref base A (0/2); new base T (62/53); total (62/55)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 117 total aligned reads displayed.

CGACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACC  >  NC_000913/712054‑712244

cGACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTa                                                                                             <  2:394136/100‑1 (MQ=255)

cGACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTa                                                                                             <  2:1268348/100‑1 (MQ=255)

   cGGATAGCGCGTAATACCGCATCCGGTTTTTCATTATGTATCCAGTTACCCGCGCCTGCAATCACATGCGCGCGTGCCTGTGGAAATTATGCCTGTatat                                                                                          >  1:3298458/1‑100 (MQ=255)

    ggATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatc                                                                                         >  1:2540171/1‑100 (MQ=255)

     gATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatca                                                                                        <  2:548990/100‑1 (MQ=255)

       tAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatcatc                                                                                      <  2:3001298/100‑1 (MQ=255)

       tAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatcatc                                                                                      <  2:1730905/100‑1 (MQ=255)

        aGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAAtcatca                                                                                     <  2:1860526/100‑1 (MQ=255)

         gcgcGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCAc                                                                                    >  2:3280300/1‑100 (MQ=255)

           gcgTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACgg                                                                                  <  1:583620/100‑1 (MQ=255)

           gcgTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACgg                                                                                  >  1:1529839/1‑100 (MQ=255)

            cgTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGt                                                                                 <  2:2575292/100‑1 (MQ=255)

            cgTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGt                                                                                 <  2:1242312/100‑1 (MQ=255)

                  aCCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCt                                                                           >  2:2762635/1‑100 (MQ=255)

                   ccGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTc                                                                          <  1:3011355/100‑1 (MQ=255)

                    cGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCg                                                                         >  2:2593181/1‑100 (MQ=255)

                     gCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGc                                                                        >  2:1853152/1‑100 (MQ=255)

                      cATCCGGTTTTTCATCATTGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCt                                                                       >  1:893275/1‑100 (MQ=255)

                       aTCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAACTTGAACCTGTAAATCATCACGGTACTGCTCGCTa                                                                      >  1:3096242/1‑100 (MQ=255)

                        tCCGGTTTTTCAGCATGGACCCAGTGGCCCGCGCCTGCAATCACATGCGCCTGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTa                                                                      <  1:1504973/99‑1 (MQ=255)

                        tCCGGTTTTTCAGCATGCACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTaa                                                                     <  1:2526025/100‑1 (MQ=255)

                          cGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAAc                                                                    >  2:1499170/1‑99 (MQ=255)

                          cGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACa                                                                   <  1:877447/100‑1 (MQ=255)

                          cGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACa                                                                   >  2:2012501/1‑100 (MQ=255)

                          cGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACa                                                                   <  1:2315973/100‑1 (MQ=255)

                            gTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATa                                                                 <  2:1364325/100‑1 (MQ=255)

                              ttttCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACg                                                               <  2:3039084/100‑1 (MQ=255)

                                ttCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGga                                                             >  1:2316518/1‑100 (MQ=255)

                                 tCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGgag                                                            >  2:795522/1‑100 (MQ=255)

                                 tCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGgag                                                            >  2:925603/1‑100 (MQ=255)

                                 tCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGgag                                                            <  1:1530650/100‑1 (MQ=255)

                                 tCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGgag                                                            >  1:1532101/1‑100 (MQ=255)

                                 tCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGgag                                                            >  2:427260/1‑100 (MQ=255)

                                 tCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGgag                                                            >  2:3397473/1‑100 (MQ=255)

                                   aGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGaa                                                          <  2:1456332/100‑1 (MQ=255)

                                   aGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGaa                                                          >  2:2150833/1‑100 (MQ=255)

                                     cATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGACTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAAtt                                                        >  2:333678/1‑100 (MQ=255)

                                       tGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTgc                                                      >  1:822873/1‑100 (MQ=255)

                                        ggACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTgcc                                                     <  2:1601201/100‑1 (MQ=255)

                                        ggACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTgcc                                                     <  2:2293485/100‑1 (MQ=255)

                                        ggACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTgcc                                                     <  1:2343969/100‑1 (MQ=255)

                                         gACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTgccg                                                    >  1:1172493/1‑100 (MQ=255)

                                             cAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAg                                                <  1:1055098/100‑1 (MQ=255)

                                             cAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAg                                                >  2:2371119/1‑100 (MQ=255)

                                               gTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAggg                                              <  1:365785/100‑1 (MQ=255)

                                                tGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGa                                             <  2:2606172/100‑1 (MQ=255)

                                                  aCCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATa                                           <  1:3067681/100‑1 (MQ=255)

                                                   nccGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATaa                                          >  2:901729/2‑100 (MQ=255)

                                                   cccGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATaa                                          >  1:1210051/1‑100 (MQ=255)

                                                     cgcgCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAc                                        >  2:3071474/1‑100 (MQ=255)

                                                       cgcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcag                                      <  2:2316550/100‑1 (MQ=255)

                                                         ccTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcaggg                                    <  1:1814751/100‑1 (MQ=255)

                                                          cTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggc                                   <  1:3383819/100‑1 (MQ=255)

                                                            gCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcag                                 >  1:1680224/1‑100 (MQ=255)

                                                             cAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcagg                                >  1:781852/1‑100 (MQ=255)

                                                             cAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcagg                                >  2:2969920/1‑100 (MQ=255)

                                                             cAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcagg                                >  2:1230829/1‑100 (MQ=255)

                                                              ccgatcttGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcaggg                               <  2:1130849/93‑1 (MQ=255)

                                                              aaTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcaggg                               <  1:109962/100‑1 (MQ=255)

                                                               aTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATCCGGAGAATTGCCGCCAGGGATAAAcagggcaggg                               >  1:3223904/1‑99 (MQ=255)

                                                               aTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGt                              <  2:1796249/100‑1 (MQ=255)

                                                               aTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGt                              >  1:542282/1‑100 (MQ=255)

                                                                 caccTGCGCCCTTGCCGGTGGAAATGCAGTGTGTAAATAATCATGGTACTGGTCGCTAACATACGGAGAATTGCTGCCAGGGATAAActtgcca                                  <  1:459577/94‑7 (MQ=255)

                                                                 cacaTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcaggg                               >  1:441135/1‑97 (MQ=255)

                                                                 cacaTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGa                            >  1:1907723/1‑100 (MQ=255)

                                                                 cacaTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGa                            >  2:2485255/1‑100 (MQ=255)

                                                                  acaTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACAGCTCGCTAACATACGGAGAATTGCCGACAGGGATAAACAGGGCAGGGTGa                            >  1:1131476/1‑99 (MQ=255)

                                                                    aTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATcc                         <  1:138918/100‑1 (MQ=255)

                                                                    aTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATcc                         <  2:1748109/100‑1 (MQ=255)

                                                                     tGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATccc                        >  1:2896913/1‑100 (MQ=255)

                                                                     tGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGAGATcgg                        >  1:1130834/1‑98 (MQ=255)

                                                                     ggcgcCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATccc                        <  1:2066311/99‑1 (MQ=255)

                                                                         cccGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGc                    <  2:3285691/100‑1 (MQ=255)

                                                                         cccGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGc                    >  1:520142/1‑100 (MQ=255)

                                                                         cccGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGc                    >  2:1882838/1‑100 (MQ=255)

                                                                          ccGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCACGGATAAACAGGGCAGGGTGATCCCATGcc                   >  1:2095715/1‑100 (MQ=255)

                                                                              gCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCTCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGa               <  2:1690901/100‑1 (MQ=255)

                                                                               ccTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAt              <  2:1941636/100‑1 (MQ=255)

                                                                               ccTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAt              <  1:32063/100‑1 (MQ=255)

                                                                               ccTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAt              <  2:2927789/100‑1 (MQ=255)

                                                                                cTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAtt             <  2:893287/100‑1 (MQ=255)

                                                                                 tgtgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAttt            >  1:54279/1‑100 (MQ=255)

                                                                                 tgtgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAttt            <  2:1150554/100‑1 (MQ=255)

                                                                                  gtgGAAATTGAGCCTGTAAATCATCACGGTTCTGCTCGCTAACCTACGGAGAATCGCCGCCAGCGATAACCAGGGCCGGGTGATCCCATTCCGGCAttt            >  1:2639384/1‑99 (MQ=255)

                                                                                  gtgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCAATGCCGGGAtttt           >  2:2093018/1‑100 (MQ=255)

                                                                                   tgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtc          <  1:1415369/100‑1 (MQ=255)

                                                                                   tgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtc          >  2:1670261/1‑100 (MQ=255)

                                                                                   tgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtc          <  2:2240163/100‑1 (MQ=255)

                                                                                    ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtct         <  2:1257412/100‑1 (MQ=255)

                                                                                    ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtct         >  1:301666/1‑100 (MQ=255)

                                                                                    ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtct         >  2:1143758/1‑100 (MQ=255)

                                                                                     gAAATTGAGCCTGTAAATCATCACGGTGCTGCTCGCGAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtctc        <  1:1669560/100‑1 (MQ=255)

                                                                                      aaaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTcc       <  2:1907751/100‑1 (MQ=255)

                                                                                      aaaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTcc       <  2:1654375/100‑1 (MQ=255)

                                                                                        aTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCa     <  2:1172508/100‑1 (MQ=255)

                                                                                        aTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGAAGCGTGATCCCATGCTGGGATTTTCTCCCa     >  1:875670/1‑100 (MQ=255)

                                                                                          tGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAc   >  1:2957490/1‑100 (MQ=255)

                                                                                           gAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAcc  >  2:3126790/1‑100 (MQ=255)

                                                                                           gAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAcc  >  1:2749694/1‑100 (MQ=255)

                                                                                           gAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAcc  >  1:995214/1‑100 (MQ=255)

CGACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACC  >  NC_000913/712054‑712244

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