BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,567,608	Δ1 bp	39.6%	coding (500/1227 nt)	eutH ←	ethanolamine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,567,607	0	A	.	39.6%	325.4 / 246.1	159	coding (501/1227 nt)	eutH	ethanolamine transporter
Reads supporting (aligned to +/- strand): ref base A (49/46); new base . (38/25); total (87/72)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

TCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGAAAACACAATC  >  NC_000913/2567510‑2567706

 cATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATc                                                                                                  >  2:1766691/1‑100 (MQ=255)

  aTGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCg                                                                                                 <  1:1771924/100‑1 (MQ=255)

   tGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACTTAGCAACCAGACCACCAGCGATACAACCAATc                                                                                                  <  1:1651911/98‑1 (MQ=255)

   tGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCgg                                                                                                <  2:3885523/100‑1 (MQ=255)

   tGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGa                                                                                               >  2:3964049/1‑100 (MQ=255)

     ttCATCAGGATCAGGGCGAAAGTGTATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAAt                                                                                             <  1:4114971/100‑1 (MQ=255)

      tcatccGGATCAGGGCGAAAGTGTATTCCACCGGCTGGCCGTTGATCAGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACGTCC‑ATCGGaan                                                                                             <  1:510276/99‑2 (MQ=255)

       catcaGGCTCAGGGCGAAAGTGAATTCCCCCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACTAGCCCACCAGCGATACAACCAATCGGAAt                                                                                             <  2:497095/99‑1 (MQ=255)

       catcaGGATCAGGGCGAAAGTGAATTCCACCGGCTTGCCGTTGATCTGCACACCGGAGTACAGAGCAGCCAGACCACCAGCGATACAACCAATCGGAATg                                                                                            <  1:2083061/100‑1 (MQ=255)

       catcaGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGaa                                                                                              >  1:3507204/1‑98 (MQ=255)

       catcaGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAA‑CTATCGGAATgg                                                                                           <  1:3009489/100‑1 (MQ=255)

        atcaggatcaggGCGAAAGTGAATGCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATgg                                                                                           <  1:553688/100‑1 (MQ=255)

            ggatcaggGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTcac                                                                                       <  1:163460/100‑1 (MQ=255)

             gatcaggGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACaa                                                                                     >  2:2694128/1‑100 (MQ=255)

                caggGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATg                                                                                   <  2:2288188/100‑1 (MQ=255)

                  gggCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGcc                                                                                 >  2:137814/1‑100 (MQ=255)

                  gggCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGcc                                                                                 >  1:1034933/1‑100 (MQ=255)

                    gCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATAccc                                                                                >  1:2245856/1‑99 (MQ=255)

                      gAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGcc                                                                             <  1:1731841/100‑1 (MQ=255)

                      gAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGcc                                                                             <  2:2511104/100‑1 (MQ=255)

                         aGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCa                                                                         >  1:1910350/1‑100 (MQ=255)

                         aGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCa                                                                         >  2:74677/1‑100 (MQ=255)

                         aGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCa                                                                         >  1:1910291/1‑100 (MQ=255)

                          gTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCa                                                                         <  2:27503/100‑1 (MQ=255)

                          gTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCAc                                                                        >  1:2793357/1‑100 (MQ=255)

                             aaTTCCACCGGCTGGCCTTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGc                                                                      <  1:2821754/100‑1 (MQ=255)

                              aTTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCAc                                                                        >  1:2837104/1‑97 (MQ=255)

                              aTTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGcc                                                                     >  1:2066385/1‑100 (MQ=255)

                              aTTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCg                                                                    <  2:3375450/100‑1 (MQ=255)

                                      ggCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCa                                                             <  2:1784848/100‑1 (MQ=255)

                                      ggCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAg                                                            >  2:322048/1‑100 (MQ=255)

                                      ggCTGGCCGTTGATCTGCACACCGGAGTACACAGCAACCAGCCCACCAGCGATACAACCAATCGAAATGGGCACAATGCCCGCCAGCACGCCGAgagcga                                                             >  1:3444965/1‑98 (MQ=255)

                                        cTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGAt                                                          >  2:1079293/1‑100 (MQ=255)

                                           gCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAac                                                       >  2:2791120/1‑100 (MQ=255)

                                           gCCGTTGATCTACACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATaa                                                        >  2:605574/1‑100 (MQ=255)

                                            ccGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAac                                                       <  1:1959372/100‑1 (MQ=255)

                                            ccGTTGATCTGCACACCGGAGTACATAGCAACCAGACAACCAGCGATACAACTAATCGGAATGGTCACAATACCCTCAAGCACGCCGAGAGCCAGATaa                                                        >  2:467010/1‑99 (MQ=255)

                                             cGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAacg                                                      >  1:599870/1‑100 (MQ=255)

                                             cGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAacga                                                     >  1:3821062/1‑100 (MQ=255)

                                                tGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACgg                                                  <  1:2304191/100‑1 (MQ=255)

                                                 gATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACgg                                                  >  1:783789/1‑100 (MQ=255)

                                                 gATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCCCCCGCACGCCGAGCGCCAGATAACGACGGt                                                 >  1:3696162/1‑100 (MQ=255)

                                                   tctcCACACCGGAGTACATAGCAACCAGACCACCAGAGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCa                                               <  2:3004824/100‑1 (MQ=255)

                                                   tCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTc                                                <  2:3800082/100‑1 (MQ=255)

                                                   tCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTc                                                <  1:1045545/100‑1 (MQ=255)

                                                   tCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGt                                                 >  1:1653334/1‑98 (MQ=255)

                                                   tCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCa                                               >  2:924608/1‑100 (MQ=255)

                                                    cTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACCACC‑GTCGGAATGGTCGCAATGCCCGCCAGCCCGCCGAGCGCCAGATAACGACGGTc                                                >  1:2119267/1‑98 (MQ=255)

                                                    cTGCACACCGGAATACATAACAACAAAACCACCAGCGATACAACC‑ATCGGAATAGACACAATACCCGCCAACACACCGAACCCCAGATAacaa                                                     >  1:1144563/1‑91 (MQ=255)

                                                     tGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAg                                              >  2:1079129/1‑100 (MQ=255)

                                                     tGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGACAGa                                             >  1:2610267/1‑100 (MQ=255)

                                                      gCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGaa                                            >  2:3282628/1‑100 (MQ=255)

                                                       cacaCCGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAg                                           <  1:2409695/100‑1 (MQ=255)

                                                         cacCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAgg                                          >  2:1990747/1‑100 (MQ=255)

                                                         cacCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAgg                                          >  2:1406303/1‑100 (MQ=255)

                                                           ccGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGtt                                        <  1:2436290/100‑1 (MQ=255)

                                                           ccGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGtt                                        >  2:2866591/1‑100 (MQ=255)

                                                           ccGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGtt                                        >  1:3111868/1‑100 (MQ=255)

                                                           ccGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGCCACAATGCCCGCAAGCACGCCGAGCGCCAGATAACGACGGTCAGAAgggt                                        >  1:1157927/1‑98 (MQ=255)

                                                           ccGGAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGTAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTc                                       <  1:869951/100‑1 (MQ=255)

                                                            cGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTc                                       <  2:2195297/100‑1 (MQ=255)

                                                             ggAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGa                                     >  1:328018/1‑100 (MQ=255)

                                                             ggAGTACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGa                                     >  2:486757/1‑100 (MQ=255)

                                                                gTACATAGCAACCAGACCACCAGCGATACAACCAACCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGAt                                    <  1:316351/99‑1 (MQ=255)

                                                                 tACATAGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAAt                                 <  2:3083166/100‑1 (MQ=255)

                                                                  aCATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATACCGACGGTCAGAAGGTTCGATAAt                                 >  1:3014983/1‑100 (MQ=255)

                                                                  aCATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAAt                                 >  2:1813133/1‑100 (MQ=255)

                                                                  aCATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAAt                                 >  2:4140994/1‑100 (MQ=255)

                                                                    aTAGCAACCAGACCACCAGCGCTACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGc                               >  2:2590097/1‑100 (MQ=255)

                                                                    aTAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGc                               <  2:1184919/100‑1 (MQ=255)

                                                                    aTAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGc                               >  1:3721433/1‑100 (MQ=255)

                                                                    aTAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGc                               <  2:296528/100‑1 (MQ=255)

                                                                      aGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCg                             >  1:1398445/1‑99 (MQ=255)

                                                                      aGCAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGa                            >  2:3776465/1‑100 (MQ=255)

                                                                       gCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGa                            >  2:3222184/1‑100 (MQ=255)

                                                                        cAACCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAgc                          <  1:2891832/100‑1 (MQ=255)

                                                                         aaCCAGACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAgcg                         <  1:4097220/100‑1 (MQ=255)

                                                                          aCCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAgcg                         >  1:2858271/1‑100 (MQ=255)

                                                                          aCCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAgcg                         <  1:774397/100‑1 (MQ=255)

                                                                          aCCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATACCGAgc                          >  2:1833392/1‑99 (MQ=255)

                                                                            cAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAgcgc                        <  1:2532052/99‑1 (MQ=255)

                                                                              gACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCAcc                    <  1:547398/100‑1 (MQ=255)

                                                                              gACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCAcc                    <  1:3264414/100‑1 (MQ=255)

                                                                              gACCACCAGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCAcc                    <  2:1172286/100‑1 (MQ=255)

                                                                                ccaccaGCGATACAACC‑ATCGGAATGGTCACAATGCCAGCCAGCACGCCGAGCGCCAGATAACGACGGTAAGAAGGTTCGATAATGCCGAGCGCCACCgg                  >  2:3563344/1‑100 (MQ=255)

                                                                                  accaGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGa                 >  1:3460330/1‑100 (MQ=255)

                                                                                    naGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGACAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAAt               >  1:389618/2‑100 (MQ=255)

                                                                                    caGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATg              <  2:2128492/100‑1 (MQ=255)

                                                                                    caGCGATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATg              >  1:295919/1‑100 (MQ=255)

                                                                                      gCGATACAACC‑ATCGGAATGGTAACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGCCGGTCAGAAGGTTCGATAATGCCGCGCGCAACCGGAATGGa            >  1:2808044/1‑100 (MQ=255)

                                                                                        gATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCGGAAGGTTCGATAATGCCGAGCGCCACCGGagatcgg           >  2:2713678/1‑94 (MQ=255)

                                                                                        gATACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGCTCAGAAGGTCCGCTACTGCCGCCCGCCACCGGAATGGaa           >  1:360435/1‑99 (MQ=255)

                                                                                         ctACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGaaa          <  1:2590089/99‑1 (MQ=255)

                                                                                          tACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGaaaa         <  2:3632445/100‑1 (MQ=255)

                                                                                          tACAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGAAAac        <  2:3540773/100‑1 (MQ=255)

                                                                                           aCAACC‑ATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGAAAaca       <  2:1408654/100‑1 (MQ=255)

                                                                                            cAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGAAAaca       >  2:1473783/1‑100 (MQ=255)

                                                                                            cAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGAAAaca       <  2:3508837/100‑1 (MQ=255)

                                                                                                cAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGAAAACACAAt   >  2:2804875/1‑100 (MQ=255)

                                                                                                cAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGAAAACACAAt   <  1:2709885/100‑1 (MQ=255)

TCATGTTCATCAGGATCAGGGCGAAAGTGAATTCCACCGGCTGGCCGTTGATCTGCACACCGGAGTACATAGCAACCAGACCACCAGCGATACAACCAATCGGAATGGTCACAATGCCCGCCAGCACGCCGAGCGCCAGATAACGACGGTCAGAAGGTTCGATAATGCCGAGCGCCACCGGAATGGAAAACACAATC  >  NC_000913/2567510‑2567706

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