BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,113,474	T→A	16.4%	L61* (TTG→TAG) ‡	uspD →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,474	0	T	A	16.4%	398.1 / 67.0	158	L61* (TTG→TAG)	uspD	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base T (65/67); new base A (15/11); total (80/78)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

Only 100 of 178 total aligned reads displayed.

GACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113375‑4113573

gACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                                                                                                     >  2:2009049/1‑100 (MQ=255)

gACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                                                                                                     >  1:3460443/1‑100 (MQ=255)

  cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                   <  2:1760863/100‑1 (MQ=255)

  cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                   <  1:1339167/100‑1 (MQ=255)

   aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                  <  1:3743944/100‑1 (MQ=255)

    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 <  1:567954/100‑1 (MQ=255)

    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 <  1:2925405/100‑1 (MQ=255)

    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 >  1:1152445/1‑100 (MQ=255)

    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 >  2:83493/1‑100 (MQ=255)

    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAACCGAATTGTATCCGGGTATCTACTTCCCTTCAACAGAAGATATTCTTCACtcgttt                                                                                                 >  1:3087256/1‑99 (MQ=255)

    tAATGACGCTCACCTGACGTTAATGCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 <  2:1705162/100‑1 (MQ=255)

      aTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaa                                                                                               <  2:485755/100‑1 (MQ=255)

       tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                              >  1:1278559/1‑100 (MQ=255)

       tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                              <  1:2561400/100‑1 (MQ=255)

       tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATACTCTTCAATTGTTgaag                                                                                              >  1:2227819/1‑100 (MQ=255)

          cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGGTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                           >  2:1657047/1‑100 (MQ=255)

           gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                          <  2:2297651/100‑1 (MQ=255)

           gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                          <  1:25196/100‑1 (MQ=255)

           gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                          <  2:2749319/100‑1 (MQ=255)

           gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATa                                                                                          <  2:1378032/100‑1 (MQ=255)

             tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                        <  1:3622164/100‑1 (MQ=255)

              cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGt                                                                                       <  2:2110385/100‑1 (MQ=255)

              cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatct                                                                                       <  2:2227927/100‑4 (MQ=38)

                ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTTCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                     >  1:449582/1‑100 (MQ=255)

                ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                     <  1:1188918/100‑6 (MQ=38)

                 cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                    <  2:2435547/100‑1 (MQ=255)

                 cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                    >  2:205100/1‑100 (MQ=255)

                 cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                    <  2:2529564/100‑7 (MQ=38)

                    aCGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATa                                                                                 <  1:3062113/100‑1 (MQ=255)

                     cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCACTGCAACAGAAGATGTTCTTCAATTATTGAAGAACAAGCCgacta                                                                                 >  2:3424804/1‑95 (MQ=255)

                      gTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                <  1:1904334/99‑1 (MQ=255)

                         aaTTCATATTGATGATGGCTTAAGCACATTGTATCCGGATATCTAATTTCCTGCAACAGAAGATATTCTTCAATTATTGAAGAATAAATCGGATaaaa                                                                              >  1:1529242/1‑96 (MQ=255)

                          aTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                           >  2:3775583/1‑100 (MQ=255)

                           ttCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCAGAtcggaagag                                                                          >  2:1250561/1‑92 (MQ=255)

                             cATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGt                                                                        >  1:1742736/1‑100 (MQ=255)

                                atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                     <  2:2778049/100‑1 (MQ=255)

                                  tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaan                                                                   <  2:1145817/100‑2 (MQ=255)

                                   gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                  <  1:2915576/100‑1 (MQ=255)

                                   gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAACAAGTCGGATAACCAGCTGTATAAAc                                                                  <  2:2508452/100‑1 (MQ=255)

                                    atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                 <  2:3103449/100‑1 (MQ=255)

                                     tgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                                >  1:1046904/1‑100 (MQ=255)

                                       atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                              <  2:1088893/100‑1 (MQ=255)

                                       atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                              <  2:279952/100‑1 (MQ=255)

                                       atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                              <  1:1082223/100‑1 (MQ=255)

                                       atgGCTTAAGCGAGTTGTATCCGGGTATCGACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                              <  1:2656388/100‑1 (MQ=255)

                                            ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                         >  1:3691730/1‑100 (MQ=255)

                                            ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                         <  1:3666842/100‑1 (MQ=255)

                                            ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                         >  2:3096651/1‑100 (MQ=255)

                                            ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                         >  2:1889270/1‑100 (MQ=255)

                                            ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                         <  2:1268553/100‑1 (MQ=255)

                                               aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                      >  2:2795457/1‑100 (MQ=255)

                                                  gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                   <  1:1681874/100‑1 (MQ=255)

                                                    gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                 >  1:1160562/1‑100 (MQ=255)

                                                      tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAAGAAGCTGTATAAAATGACGAAAAATAGTAAATg                                               >  1:3349585/1‑100 (MQ=255)

                                                      tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAAt                                                >  1:3385660/1‑99 (MQ=255)

                                                       gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                              >  2:2597432/1‑100 (MQ=255)

                                                       gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                              >  1:2553785/1‑100 (MQ=255)

                                                       gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                              >  2:2117926/1‑100 (MQ=255)

                                                       gTATCAGGGTATCTACGACCCTGCAACAGAAGATATTCTTCAAAAGAAGAAGAATAAGTCGGATATCAAGCTGAAAAAAAAAAAAAAAAATATTCAATgg                                              <  1:1424113/100‑1 (MQ=255)

                                                        tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                             <  2:2149207/100‑1 (MQ=255)

                                                        tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                             >  2:2134131/1‑100 (MQ=255)

                                                          tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                           <  2:3574744/100‑1 (MQ=255)

                                                          tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                           <  2:1727630/100‑1 (MQ=255)

                                                           ccGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGa                                          >  1:3490715/1‑100 (MQ=255)

                                                             gggTATCTACTTCCGTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                        <  1:2597419/100‑1 (MQ=255)

                                                              tgTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                       <  2:3313186/99‑1 (MQ=255)

                                                               gTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGAc                                      <  1:2245459/100‑1 (MQ=255)

                                                                tATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACa                                     >  2:3619410/1‑100 (MQ=255)

                                                                 aTCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaa                                    <  2:2338563/100‑1 (MQ=255)

                                                                 aTCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAAAAGTCAGATAACAACCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGAc                                      <  1:2354385/98‑1 (MQ=255)

                                                                 aTCTACTGCCCTGCAACAGAAGATATTCTTCAATTGTTAAAGAATAAGTCGGATAACAAGCTGTATAAAAAAAAGAAAAATATTCAATGGCCGAAGa                                       <  1:353291/97‑1 (MQ=255)

                                                                    tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                 <  1:946246/100‑1 (MQ=255)

                                                                      cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                               <  2:3132691/100‑1 (MQ=255)

                                                                        tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                             <  2:3398705/100‑1 (MQ=255)

                                                                         cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                            >  1:789474/1‑100 (MQ=255)

                                                                          ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                           >  2:3647785/1‑100 (MQ=255)

                                                                          ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                           >  2:187278/1‑100 (MQ=255)

                                                                                  agaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                   >  2:121898/1‑100 (MQ=255)

                                                                                   gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  1:926452/1‑100 (MQ=255)

                                                                                   gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  2:2946284/1‑100 (MQ=255)

                                                                                   gaagaTATTATTGAATTGTTGAAGAATAAGTTGGATAACAAGATGTATAAAATGAGGAAAAATATTCAATGGACGAAGACAAAACTGCGTATTGAAcgag                  >  1:3174746/1‑98 (MQ=255)

                                                                                   cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  2:864815/6‑100 (MQ=38)

                                                                                   cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                     >  2:1792739/6‑97 (MQ=38)

                                                                                     acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  2:759689/4‑100 (MQ=38)

                                                                                     acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  2:3734087/4‑100 (MQ=38)

                                                                                      cttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag               >  1:1962379/3‑100 (MQ=255)

                                                                                       ttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgaga              >  2:3004822/2‑100 (MQ=255)

                                                                                        tatTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa             >  2:406665/1‑100 (MQ=255)

                                                                                        tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa             >  1:3265515/1‑100 (MQ=255)

                                                                                        tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATAGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa             >  1:2742955/1‑100 (MQ=255)

                                                                                        tatTCTTCAACAATTGAAGAACAAGTCGGATAACAAGCCGTATAAACTGACGAAAAACATTCAATGGCCGAAGACAAAACTGCGTAATGAACGCGcaga              >  1:1155309/1‑99 (MQ=255)

                                                                                           tcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg          >  2:3416772/1‑100 (MQ=255)

                                                                                           tcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg          <  2:2286238/100‑1 (MQ=255)

                                                                                            cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGc         >  1:1121212/1‑100 (MQ=255)

                                                                                              tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg       <  2:3472080/100‑1 (MQ=255)

                                                                                               cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      <  2:1946101/100‑1 (MQ=255)

                                                                                               cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      <  2:436321/100‑1 (MQ=255)

                                                                                                aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAATGACAAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa     <  1:2470819/100‑1 (MQ=255)

                                                                                                   tgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAac  >  1:1208404/1‑100 (MQ=255)

                                                                                                   tgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAac  <  1:187274/100‑1 (MQ=255)

GACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113375‑4113573

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