BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,880,549	T→C	T71A (ACC→GCC)	leuE ←	leucine efflux protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,880,549	0	T	C	98.2%	354.3 / ‑3.0	109	T71A (ACC→GCC)	leuE	leucine efflux protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (61/45); total (61/48)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

Only 100 of 112 total aligned reads displayed.

TTTTTACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACC  >  NC_000913/1880452‑1880645

tttttACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCt                                                                                                <  2:2267716/100‑1 (MQ=255)

 ttttACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCtt                                                                                               <  2:2460658/100‑1 (MQ=255)

  tttACCCTTCAGGGTTGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCtt                                                                                               >  1:1721256/1‑99 (MQ=255)

  tttACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCtta                                                                                              >  1:1228319/1‑100 (MQ=255)

    tACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaat                                                                                            <  2:41120/100‑1 (MQ=255)

    tACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaat                                                                                            <  2:173259/100‑1 (MQ=255)

    tACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaat                                                                                            >  2:1877386/1‑100 (MQ=255)

     aCCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaatt                                                                                           <  1:2633367/100‑1 (MQ=255)

       ccTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaattaa                                                                                         >  1:2807487/1‑100 (MQ=255)

        cTTCAGGGTCGCGTAAAGAATTTTTCTCCCCAGATAGAGCAAATAAAACGCATCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaatta                                                                                          <  1:1523748/98‑1 (MQ=255)

        cTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGCTAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGAGTGGCCTtaattaat                                                                                        >  1:1198236/1‑100 (MQ=255)

        cTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTtaattaat                                                                                        >  1:1919785/1‑100 (MQ=255)

        cTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTtaattaat                                                                                        >  2:1981201/1‑100 (MQ=255)

         ttCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATg                                                                                       >  2:3050692/1‑100 (MQ=255)

           cAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTc                                                                                     <  2:1396512/100‑1 (MQ=255)

            aGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCg                                                                                    <  2:1908552/100‑1 (MQ=255)

            aGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCg                                                                                    <  2:3317885/100‑1 (MQ=255)

              ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGcc                                                                                  <  2:2610848/100‑1 (MQ=255)

              ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGcc                                                                                  <  1:871985/100‑1 (MQ=255)

              ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGcc                                                                                  <  2:1198872/100‑1 (MQ=255)

               gTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCa                                                                                 >  2:1965300/1‑100 (MQ=255)

                 cgcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACt                                                                               <  2:599806/100‑1 (MQ=255)

                 cgcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACt                                                                               >  2:1038739/1‑100 (MQ=255)

                 cgcgTAAAGAATTTTAATCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACt                                                                               <  1:1639843/100‑1 (MQ=255)

                  gcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACt                                                                               >  1:1225479/1‑99 (MQ=255)

                   cgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTc                                                                              >  1:3470199/1‑99 (MQ=255)

                    gTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCa                                                                            >  2:219160/1‑100 (MQ=255)

                     tAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAg                                                                           <  1:1999394/100‑1 (MQ=255)

                     tAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAg                                                                           <  1:500978/100‑1 (MQ=255)

                      aaaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGc                                                                          <  1:2281457/100‑1 (MQ=255)

                      aaaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGc                                                                          >  1:3245999/1‑100 (MQ=255)

                       aaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGcc                                                                         >  1:2575717/1‑100 (MQ=255)

                       aaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGcc                                                                         >  1:2715219/1‑100 (MQ=255)

                          aaTTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCAt                                                                      <  2:2574030/100‑1 (MQ=255)

                            ttttACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGc                                                                    <  1:2278117/100‑1 (MQ=255)

                             tttACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGcn                                                                   <  1:773828/100‑2 (MQ=255)

                             tttACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGcc                                                                   <  2:1906042/100‑1 (MQ=255)

                                 cTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGaa                                                               >  1:2796145/1‑100 (MQ=255)

                                  tCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGACTAATATCGGGGTGGCCTTAATTAATGACGCCACTCCAGACCATGCCAGaaa                                                              >  1:195983/1‑100 (MQ=255)

                                   ccccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAAc                                                             <  2:2128058/100‑1 (MQ=255)

                                   ccccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAAc                                                             >  2:2616820/1‑100 (MQ=255)

                                   ccccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAAc                                                             >  1:3188055/1‑100 (MQ=255)

                                   ccccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAAc                                                             >  2:2380781/1‑100 (MQ=255)

                                   ccccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCAAGCCCATGCCAGAAAc                                                             >  1:2646139/1‑100 (MQ=255)

                                   ccccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGGGGCCTTAATTAATGTCGCCACTCCACCCCATGCCAGAAAc                                                             >  1:2453806/1‑100 (MQ=255)

                                   ccccAGATAGAGCAAATAAAACGCACCAAGATAACGGACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTTCAGCCCATGCCAGAAAc                                                             >  1:156403/1‑100 (MQ=255)

                                    cccAGATAGAGCAAATAAAACGCATCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACa                                                            <  1:1334742/100‑1 (MQ=255)

                                    cccAGATAGAGCAAATAAAACGCACCAAGATAACGTACCATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACa                                                            >  1:2135596/1‑100 (MQ=255)

                                    cccAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACa                                                            >  2:2267505/1‑100 (MQ=255)

                                      cAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATc                                                          <  1:1490538/100‑1 (MQ=255)

                                      cAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATc                                                          <  1:2312929/100‑1 (MQ=255)

                                      cAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGAGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATc                                                          >  1:2784719/1‑100 (MQ=255)

                                       aGATAGAGCAAGTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAAc                                                             >  2:33225/1‑96 (MQ=255)

                                          tAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATa                                                      >  1:1889986/1‑100 (MQ=255)

                                            gagCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                    >  2:3321512/1‑100 (MQ=255)

                                            gagCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                    >  2:877400/1‑100 (MQ=255)

                                              gCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGc                                                  <  2:2852610/100‑1 (MQ=255)

                                              gCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGc                                                  >  2:2121996/1‑100 (MQ=255)

                                               cAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCa                                                 >  1:984965/1‑100 (MQ=255)

                                                aaaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCAt                                                <  1:1159494/100‑1 (MQ=255)

                                                aaaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCATTAATTAATGTCGCCGCTCCAGCCCATGCCAGAAACATCAATACCTCa                                                 >  2:2226374/1‑99 (MQ=255)

                                                  aTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCg                                              <  1:2935145/100‑1 (MQ=255)

                                                  aTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCg                                              >  2:3236487/1‑100 (MQ=255)

                                                  aTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCATTAATTAATGTAACCAATCAAGCCCATGCGAGGAGCATCAATACAGCATCg                                              >  2:3405930/1‑100 (MQ=255)

                                                   tAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGc                                             <  2:334179/100‑1 (MQ=255)

                                                    aaaaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGcc                                            <  2:1559533/100‑1 (MQ=255)

                                                     gatcGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCa                                           <  2:1336634/97‑1 (MQ=255)

                                                      aaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaa                                          >  2:1968980/1‑100 (MQ=255)

                                                       aCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaat                                         <  1:2184455/100‑1 (MQ=255)

                                                        cGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaata                                        <  1:877388/100‑1 (MQ=255)

                                                         gCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataa                                       >  1:626188/1‑100 (MQ=255)

                                                          cACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataaa                                      >  1:1926968/1‑100 (MQ=255)

                                                          cACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaagatc                                      >  1:1336615/1‑96 (MQ=255)

                                                            ccAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataaata                                    <  2:1551420/100‑1 (MQ=255)

                                                            ccAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataaata                                    <  1:2917273/100‑1 (MQ=255)

                                                            ccAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataaata                                    >  1:1748812/1‑100 (MQ=255)

                                                              aaGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATaca                                  >  2:2798461/1‑100 (MQ=255)

                                                                     cGTACAATGTTTAATAATATGCAGGTGGACCCAATCAATGTCGCCACTCCAGCCCATGACAGACACATCAATACCGCATCGCCAATAAATACTCCGc                              <  1:1109928/97‑1 (MQ=255)

                                                                     cGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAgg                           <  1:1807615/100‑1 (MQ=255)

                                                                      gTAGAATGTTGAATAACATCGGGGGGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAAGACCGCATCGCCAATAAAGACACCGCAgg                           <  1:3252016/99‑1 (MQ=255)

                                                                      gTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGc                          >  2:1282869/1‑100 (MQ=255)

                                                                      gTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGc                          >  1:3260445/1‑100 (MQ=255)

                                                                      gTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTAGCCAATCCAGAACATGCCAGAAACATCAATACCGCATCGCCAATAAATACACAGCAGGc                          >  1:1586196/1‑100 (MQ=255)

                                                                       tACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGcc                         >  1:1243802/1‑100 (MQ=255)

                                                                       caCAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGcc                         >  1:3297219/2‑100 (MQ=255)

                                                                        aCAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCg                        >  1:2496243/1‑100 (MQ=255)

                                                                         cAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGc                       <  2:2575751/100‑1 (MQ=255)

                                                                             gTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGa                   >  2:3424124/1‑100 (MQ=255)

                                                                             gTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGa                   >  1:2110031/1‑100 (MQ=255)

                                                                              ttGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAg                    <  1:2271002/98‑1 (MQ=255)

                                                                               tGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATa                 <  2:3260496/100‑1 (MQ=255)

                                                                               tGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATa                 >  1:350220/1‑100 (MQ=255)

                                                                                gAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATaa                >  1:2544529/1‑100 (MQ=255)

                                                                                 aataatATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAAc               <  1:2267473/100‑1 (MQ=255)

                                                                                     atatCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGcc           >  2:3061387/1‑100 (MQ=255)

                                                                                         cGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTc       <  2:2646177/100‑1 (MQ=255)

                                                                                           gggTGCCCTAAATTAATGTCCCCACTCCAGCCCATGCCAGAAACACCAATACCGCATCGCCAATAACTACACCGCAGGCCGCAAGATAACCGCCTTTCAt     <  1:3188933/100‑1 (MQ=255)

                                                                                             gtggcCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATAc   <  1:1184318/100‑1 (MQ=255)

                                                                                             gtggcCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATAc   >  2:2233962/1‑100 (MQ=255)

                                                                                              tggcCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATAcc  <  2:2287698/100‑1 (MQ=255)

TTTTTACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATACC  >  NC_000913/1880452‑1880645

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: `‑`