BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,113,475	G→A	16.9%	L61L (TTG→TTA) ‡	uspD →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,475	0	G	A	16.9%	378.2 / 59.5	158	L61L (TTG→TTA)	uspD	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base G (63/68); new base A (17/10); total (80/78)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 173 total aligned reads displayed.

CATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113377‑4113573

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

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

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

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

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

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

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

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

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

        cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAAt                                                                                           <  2:3440261/100‑1 (MQ=255)

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

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

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

         gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                          <  1:25196/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:3389663/100‑4 (MQ=38)

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

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

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

              ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                     <  2:486023/100‑6 (MQ=38)

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

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

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

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

                 gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAt                                                                                  >  1:1900133/1‑100 (MQ=255)

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

                  aCGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTAAAGAATAAGTCGGATa                                                                                 <  1:664784/100‑1 (MQ=255)

                   ngTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                >  2:1196881/2‑100 (MQ=255)

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

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

                         ttCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCt                                                                          >  2:340430/1‑100 (MQ=255)

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

                           cATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGt                                                                        <  2:3129097/100‑1 (MQ=255)

                            atatTTATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGAAAAAGAAGATATTCTTCATTTGTTGAAGAATAAGTAGAATAACAAACTGta                                                                       >  1:1235697/1‑100 (MQ=255)

                            atatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                       <  1:2700463/100‑1 (MQ=255)

                              atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                     >  1:422084/1‑100 (MQ=255)

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

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

                                    gatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGa                                                               <  1:2148728/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)

                                      tgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACg                                                             <  2:3660779/100‑1 (MQ=255)

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

                                             ctcGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                      <  2:1337958/98‑1 (MQ=255)

                                               cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAataag                                                    >  1:1337950/1‑98 (MQ=255)

                                               cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                    >  2:2750500/1‑100 (MQ=255)

                                               cGAGTTGTATCCGGGTATCAACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAATAAGATAACAAGCTGTATAAACTGACGAAAAatat                                                     >  1:474789/1‑99 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                       aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                            >  2:2231561/1‑100 (MQ=255)

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

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

                                                           gggTATCTACTTCCGTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                        <  1:2597419/100‑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                                 >  2:3743228/1‑100 (MQ=255)

                                                                     ttCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                              >  1:1739836/1‑100 (MQ=255)

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

                                                                       tccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                            <  1:3647762/99‑1 (MQ=255)

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

                                                                         cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                          >  1:964320/1‑100 (MQ=255)

                                                                              aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                     <  2:3554164/100‑1 (MQ=255)

                                                                                agaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                   <  1:863087/100‑1 (MQ=255)

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

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

                                                                                 gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  <  2:1142184/100‑1 (MQ=255)

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

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

                                                                                  aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 <  2:1742745/100‑1 (MQ=255)

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

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

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

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

                                                                                      tatTCTTCAACAATTGAAGAATAAGTCGGCTAACAAGCTGTATAAACTGACGAAAAATATTCAATGTCCGAAGACAAAACTGCGTATTGAACGCGgaga              >  1:1155135/1‑99 (MQ=255)

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

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

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

                                                                                           ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc        >  1:3073768/1‑100 (MQ=255)

                                                                                           ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc        >  2:3088062/1‑100 (MQ=255)

                                                                                             cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      >  1:3138331/1‑100 (MQ=255)

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

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

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

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

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

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

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

CATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113377‑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: `‑`