BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,339,701	C→T	G265G (GGC→GGT)	yrbH →	putative isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,339,701	0	C	T	99.6%	851.6 / ‑5.4	250	G265G (GGC→GGT)	yrbH	putative isomerase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (136/112); total (138/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

Only 100 of 263 total aligned reads displayed.

AAATCTTGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGCGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTGAGGCACTGAA  >  NC_000913/3339611‑3339790

aaaTCTTGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGgtgt                                                                                         <  2:741114/93‑1 (MQ=255)

aaaTCTTGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGgtgt                                                                                         >  1:2486124/1‑93 (MQ=255)

 aaTCTTGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGtgtg                                                                                        <  2:2172611/93‑1 (MQ=255)

  aTCTTGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTgg                                                                                       >  1:3344097/1‑93 (MQ=255)

   tCTTGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGa                                                                                      <  1:169687/93‑1 (MQ=255)

   tCTTGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGa                                                                                      <  1:1831525/93‑1 (MQ=255)

   tCTTAGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGa                                                                                      <  2:1600179/93‑1 (MQ=255)

     ttGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATg                                                                                    >  1:2699204/1‑93 (MQ=255)

     ttGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATg                                                                                    >  1:2017034/1‑93 (MQ=255)

       ggTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGtt                                                                                  >  1:3599885/1‑93 (MQ=255)

       ggTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGtt                                                                                  <  2:4605172/93‑1 (MQ=255)

        gTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTc                                                                                 <  1:1158266/93‑1 (MQ=255)

         tATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCg                                                                                >  1:3822039/1‑93 (MQ=255)

         tATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCg                                                                                >  2:4355524/1‑93 (MQ=255)

           tGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTc                                                                              <  1:1889274/93‑1 (MQ=255)

           tGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTc                                                                              >  2:282658/1‑93 (MQ=255)

             aCTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAg                                                                            >  1:3971139/1‑93 (MQ=255)

              cTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGt                                                                           <  1:2332244/93‑1 (MQ=255)

              cTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGt                                                                           >  1:2573370/1‑93 (MQ=255)

              cTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGt                                                                           >  1:597562/1‑93 (MQ=255)

               tGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGtt                                                                          >  1:142183/1‑93 (MQ=255)

               tGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGtt                                                                          >  2:1338881/1‑93 (MQ=255)

                gTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTa                                                                         >  1:2210729/1‑93 (MQ=255)

                gTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTa                                                                         >  2:466379/1‑93 (MQ=255)

                gTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTa                                                                         <  1:4675471/93‑1 (MQ=255)

                 tCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTaa                                                                        <  2:3216409/93‑1 (MQ=255)

                  cATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAg                                                                       <  1:1818757/93‑1 (MQ=255)

                  cATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAg                                                                       >  1:3005311/1‑93 (MQ=255)

                   aTTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGt                                                                      <  2:2270273/93‑1 (MQ=255)

                   aTTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGt                                                                      >  1:1200801/1‑93 (MQ=255)

                   aTTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGt                                                                      >  1:4395469/1‑93 (MQ=255)

                    tttGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTa                                                                     >  1:2395612/1‑93 (MQ=255)

                    tttGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTa                                                                     <  1:3635410/93‑1 (MQ=255)

                    tttGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTa                                                                     >  1:1279928/1‑93 (MQ=255)

                    tttGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTa                                                                     >  2:2402460/1‑93 (MQ=255)

                      tGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTAtt                                                                   <  1:370606/93‑1 (MQ=255)

                      tGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTAtt                                                                   <  1:4215861/93‑1 (MQ=255)

                      tGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTAtt                                                                   <  2:1774551/93‑1 (MQ=255)

                       gCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTg                                                                  <  1:4355177/93‑1 (MQ=255)

                       gCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTg                                                                  >  1:1042860/1‑93 (MQ=255)

                         gATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTCCGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGcc                                                                <  1:4083651/93‑1 (MQ=255)

                         gATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGcc                                                                >  2:2174396/1‑93 (MQ=255)

                         gATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGcc                                                                >  2:4152577/1‑93 (MQ=255)

                          aTGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCg                                                               >  2:1879511/1‑93 (MQ=255)

                          aTGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCg                                                               >  1:2946536/1‑93 (MQ=255)

                           tGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGa                                                              <  2:167443/93‑1 (MQ=255)

                           tGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGa                                                              >  2:3701369/1‑93 (MQ=255)

                           tGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGa                                                              <  2:3904997/93‑1 (MQ=255)

                            gACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGAt                                                             >  2:3679368/1‑93 (MQ=255)

                            gACAATATGATGATTGAAGGCATCTTTACCGACCGAGATTTACGCCGTGTCTTCGGTATGCGCGTGGATGTTCTTCAGTTAAGTATTGCCCAt                                                             >  2:2838290/1‑93 (MQ=255)

                             aCAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGAtg                                                            >  2:44578/1‑93 (MQ=255)

                                atatGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtga                                                         >  1:1020758/1‑93 (MQ=255)

                                atatGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtga                                                         >  1:3439934/1‑93 (MQ=255)

                                atatGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGAATTGCCGATGtga                                                         >  1:4314874/1‑93 (MQ=255)

                                 tatGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtgat                                                        >  2:4322500/1‑93 (MQ=255)

                                 tatGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtgat                                                        >  2:602972/1‑93 (MQ=255)

                                   tgatgaTTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtgatga                                                      <  2:4520045/93‑1 (MQ=255)

                                   tgatgaTTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtgatga                                                      <  2:1039167/93‑1 (MQ=255)

                                   tgatgaTTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtgatga                                                      >  2:3819418/1‑93 (MQ=255)

                                   tgatgaTTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtgatga                                                      <  2:185182/93‑1 (MQ=255)

                                   agatgatTGAAGGCATCTTTACCGACGGTGATTTCCGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGtgatga                                                      <  2:3535745/92‑1 (MQ=255)

                                    gatgatTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGAc                                                     <  2:2435488/93‑1 (MQ=255)

                                    gatgatTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGAc                                                     >  2:1404549/1‑93 (MQ=255)

                                    gatgatTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGAc                                                     <  1:1508326/93‑1 (MQ=255)

                                    gatgatTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGACGATGtgatgcc                                                     >  1:349409/1‑91 (MQ=255)

                                      tgatTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGc                                                   <  2:2474930/93‑1 (MQ=255)

                                         ttGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCgg                                                >  2:1539115/1‑93 (MQ=255)

                                          tGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCggg                                               <  2:2138248/93‑1 (MQ=255)

                                           gAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCgggg                                              >  2:937279/1‑93 (MQ=255)

                                            aaGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCggggg                                             <  2:3405846/93‑1 (MQ=255)

                                            aaGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCggggg                                             >  2:301083/1‑93 (MQ=255)

                                            aaGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCggggg                                             >  1:3849512/1‑93 (MQ=255)

                                             aGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCgggggg                                            >  2:2338973/1‑93 (MQ=255)

                                              ggCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGa                                           >  2:3077865/1‑93 (MQ=255)

                                              ggCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACACCGGGGGGa                                           <  1:2312386/93‑1 (MQ=255)

                                               gCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGaa                                          >  1:3096224/1‑93 (MQ=255)

                                                 aTCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATa                                        <  1:2232143/93‑1 (MQ=255)

                                                  tCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATAc                                       >  1:2905747/1‑93 (MQ=255)

                                                  tCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATAc                                       >  2:4404044/1‑93 (MQ=255)

                                                   cTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGCGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACg                                      >  1:755002/1‑93 (MQ=255)

                                                     ttACCGACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGACGACGCCGGGGGGAATACgtg                                    <  2:2667672/93‑1 (MQ=255)

                                                          gACGGTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTgcgc                               >  2:2901907/1‑93 (MQ=255)

                                                             ggTGATTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCt                            >  2:196810/1‑93 (MQ=255)

                                                                 aTTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCa                        >  1:4121316/1‑93 (MQ=255)

                                                                 aTTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCa                        >  1:2130409/1‑93 (MQ=255)

                                                                 aTTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCa                        >  1:2864855/1‑93 (MQ=255)

                                                                 aTTTACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCa                        <  1:1919449/93‑1 (MQ=255)

                                                                  tttACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGCATACGTGTGCGCCCTGGCAt                       >  1:2284758/1‑93 (MQ=255)

                                                                  tttACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCAt                       >  1:2550086/1‑93 (MQ=255)

                                                                   ttACGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCAtt                      >  2:4215137/1‑93 (MQ=255)

                                                                      cGCCGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTg                   <  1:1677223/93‑1 (MQ=255)

                                                                         cGTGTCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGcc                <  1:3040615/93‑1 (MQ=255)

                                                                           tgtCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGt              <  2:1183644/93‑1 (MQ=255)

                                                                             tCTTCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTg            >  1:1955566/1‑93 (MQ=255)

                                                                               ttCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTGAg          <  2:455548/93‑1 (MQ=255)

                                                                               ttCGATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGACGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTGAg          <  1:1749427/93‑1 (MQ=255)

                                                                                  gATATGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTGAGGCa       >  1:4598117/1‑93 (MQ=255)

                                                                                   atatGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTGAGGCAc      >  1:645578/1‑93 (MQ=255)

                                                                                     atGGGTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTGAGGCACTg    <  2:172829/93‑1 (MQ=255)

                                                                                       gggTGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTGAGGCACTGaa  >  1:3184271/1‑93 (MQ=255)

AAATCTTGGTATGACTGTCATTTGCGATGACAATATGATGATTGAAGGCATCTTTACCGACGGTGATTTACGCCGTGTCTTCGATATGGGCGTGGATGTTCGTCAGTTAAGTATTGCCGATGTGATGACGCCGGGGGGAATACGTGTGCGCCCTGGCATTCTGGCCGTTGAGGCACTGAA  >  NC_000913/3339611‑3339790

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