BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	709,880	T→C	intergenic (‑11/+278)	fur ← / ← fldA	Fur/flavodoxin 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	709,880	0	T	C	99.1%	1453.5 / ‑0.4	451	intergenic (‑11/+278)	fur/fldA	Fur/flavodoxin 1
Reads supporting (aligned to +/- strand): ref base T (1/3); new base C (257/190); total (258/193)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

Only 100 of 470 total aligned reads displayed.

TGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAA  >  NC_000913/709790‑709969

tGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGccc                                                                                         <  2:3147355/93‑1 (MQ=255)

   agaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGt                                                                                      <  1:1759641/93‑1 (MQ=255)

   agaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGt                                                                                      <  2:2233188/93‑1 (MQ=255)

    gaaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGtt                                                                                     >  1:3375184/1‑93 (MQ=255)

     aaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTa                                                                                    <  2:6274363/93‑1 (MQ=255)

      aCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTAc                                                                                   >  2:1491897/1‑93 (MQ=255)

       cTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACt                                                                                  >  1:347216/1‑93 (MQ=255)

        ttCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTa                                                                                 >  1:724270/1‑93 (MQ=255)

        ttCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTa                                                                                 >  1:2358986/1‑93 (MQ=255)

         tCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGACTCTGCCCTGTTACTac                                                                                >  1:2475704/1‑92 (MQ=255)

         tCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTaa                                                                                >  2:7406547/1‑93 (MQ=255)

          ccAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAg                                                                               <  2:4740405/93‑1 (MQ=255)

            aGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCg                                                                             <  2:3585797/93‑1 (MQ=255)

            aGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCg                                                                             >  2:9156209/1‑93 (MQ=255)

            aGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGAAATCTGTCCTGTTACTAAGCg                                                                             <  2:2081745/93‑1 (MQ=255)

             ggATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCgg                                                                            >  1:9312134/1‑93 (MQ=255)

              gATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGt                                                                           >  1:5950837/1‑93 (MQ=255)

               aTTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGtt                                                                          >  2:2727576/1‑93 (MQ=255)

               aTTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGtt                                                                          >  2:8344151/1‑93 (MQ=255)

                tttttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTc                                                                         >  2:4597020/1‑93 (MQ=255)

                tttttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTc                                                                         >  2:8763783/1‑93 (MQ=255)

                  tttAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCAc                                                                       <  1:4274938/93‑1 (MQ=255)

                     aaaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAGGCGGTTCAtttt                                                                    >  2:8499904/1‑92 (MQ=255)

                     aaaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTc                                                                    >  2:7050406/1‑93 (MQ=255)

                     aaaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTc                                                                    <  1:7618068/93‑1 (MQ=255)

                      aaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCa                                                                   <  2:4887738/93‑1 (MQ=255)

                      aaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTct                                                                   >  2:2425894/1‑92 (MQ=255)

                      aaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCa                                                                   <  1:2911830/93‑1 (MQ=255)

                      aaCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCa                                                                   >  1:8863192/1‑93 (MQ=255)

                       aCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCAt                                                                  >  1:9643919/1‑93 (MQ=255)

                       aCGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCAt                                                                  >  2:8221832/1‑93 (MQ=255)

                         gAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTTACTTCATTa                                                                >  1:6721807/1‑93 (MQ=255)

                         gAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTa                                                                >  1:6366742/1‑93 (MQ=255)

                         gAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTa                                                                >  2:4058959/1‑93 (MQ=255)

                         gAGGAAGCGTTACTTTCACGCCAGCTTTCCTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGACCTGTTACTAAGCGGTTCACTTCATTa                                                                >  2:1651405/1‑93 (MQ=255)

                          aGGAAGCGTTACTTTCAGGCCATCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGTATCTGCCCTGTTACTAAGCGGTTCACTTCATTag                                                               >  1:7686412/1‑93 (MQ=255)

                          aGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTag                                                               >  2:9722298/1‑93 (MQ=255)

                          aGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTag                                                               <  1:9853404/93‑1 (MQ=255)

                          aGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTag                                                               >  1:4413946/1‑93 (MQ=255)

                             aaGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTagaag                                                            >  2:782244/1‑93 (MQ=255)

                              aGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTagaaga                                                           <  1:9430786/93‑1 (MQ=255)

                              aGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTagaaga                                                           >  2:9365653/1‑93 (MQ=255)

                               gCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAgaagaa                                                          >  2:2797385/1‑93 (MQ=255)

                               gCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAgaagaa                                                          <  2:6751199/93‑1 (MQ=255)

                                 gTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGt                                                        >  2:1171458/1‑93 (MQ=255)

                                  ttACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                                       >  1:5011829/1‑93 (MQ=255)

                                  ttACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                                       <  1:5483109/93‑1 (MQ=255)

                                   tACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGa                                                      >  2:5261600/1‑93 (MQ=255)

                                      tttCAGGCCAGCTTTCTTTAGGGCGGTGTTGTTATCAGTCATGCGGTATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAg                                                   >  2:3879673/1‑93 (MQ=255)

                                       ttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGa                                                  >  2:3948851/1‑93 (MQ=255)

                                       ttCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGa                                                  >  1:2564508/1‑93 (MQ=255)

                                        tCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGaa                                                 <  2:5159470/93‑1 (MQ=255)

                                          aGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAAtt                                               <  2:2992712/93‑1 (MQ=255)

                                            gCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTg                                             >  1:1112048/1‑93 (MQ=255)

                                             ccAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGc                                            <  1:603192/93‑1 (MQ=255)

                                              cAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGct                                           >  1:8018623/1‑93 (MQ=255)

                                               aGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctc                                          <  1:1461648/93‑1 (MQ=255)

                                               aGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGCTTCACTTCATTAGAAGAAGTGACAGAATTTGctc                                          >  1:9519919/1‑93 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATGAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         <  2:5042703/93‑1 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         <  2:7758465/93‑1 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         >  2:5465534/1‑93 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         <  1:5694017/93‑1 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         <  2:1391285/93‑1 (MQ=255)

                                                gCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGctct                                         <  1:10022036/93‑1 (MQ=255)

                                                gCTTTCTTTAGGGCGGAATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGCGACAGAATTTGctct                                         >  1:9078372/1‑93 (MQ=255)

                                                  tttctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTg                                       <  1:7880636/93‑1 (MQ=255)

                                                   ttctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTga                                      >  1:8427852/1‑93 (MQ=255)

                                                    tctttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgag                                     >  1:3598408/1‑93 (MQ=255)

                                                       ttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATa                                  <  2:7497479/93‑1 (MQ=255)

                                                       ttAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATa                                  >  2:9480783/1‑93 (MQ=255)

                                                        tAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATaa                                 >  1:8261314/1‑93 (MQ=255)

                                                        tAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATaa                                 <  1:3333609/93‑1 (MQ=255)

                                                        tAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATaa                                 <  1:8905674/93‑1 (MQ=255)

                                                        tAGGGCGGTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATaa                                 <  1:5119898/93‑1 (MQ=255)

                                                               gTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATACTGCGTAt                          >  2:648103/1‑93 (MQ=255)

                                                               gTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTAt                          >  2:5611962/1‑93 (MQ=255)

                                                               gTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTAt                          >  2:6099030/1‑93 (MQ=255)

                                                               gTATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTAt                          <  1:9925737/93‑1 (MQ=255)

                                                                tATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATc                         >  1:8685189/1‑93 (MQ=255)

                                                                tATTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATc                         <  1:7477095/93‑1 (MQ=255)

                                                                 aTTGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCa                        >  1:1892449/1‑93 (MQ=255)

                                                                  ttGTTATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCAt                       <  1:4597021/93‑1 (MQ=255)

                                                                     ttATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtat                    >  1:6892995/1‑93 (MQ=255)

                                                                      tATCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtata                   <  2:1053848/93‑1 (MQ=255)

                                                                       aTCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAg                  >  2:6842008/1‑93 (MQ=255)

                                                                       aTCAGTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAg                  >  2:3384074/1‑93 (MQ=255)

                                                                           gTCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAtt              >  1:400309/1‑93 (MQ=255)

                                                                            tCATGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTg             >  2:6695512/1‑93 (MQ=255)

                                                                              aTGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGcc           <  1:3178501/93‑1 (MQ=255)

                                                                               tGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGCCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCa          <  1:4875648/93‑1 (MQ=255)

                                                                               tGCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCa          <  1:3924164/93‑1 (MQ=255)

                                                                                gCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCAc         <  1:8158021/93‑1 (MQ=255)

                                                                                gCGGAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCAc         <  2:1846143/93‑1 (MQ=255)

                                                                                  ggAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGc       <  1:7116072/93‑1 (MQ=255)

                                                                                   gAATCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGcc      >  2:4157730/1‑93 (MQ=255)

                                                                                     aTCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGCATTGCCACGCCTa    >  2:5686056/1‑93 (MQ=255)

                                                                                     aTCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTa    >  1:3671543/1‑93 (MQ=255)

                                                                                     aTCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTa    >  1:2730459/1‑93 (MQ=255)

                                                                                      tCTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaa   >  2:80023/1‑93 (MQ=255)

                                                                                       cTGCCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaaa  <  2:5938819/93‑1 (MQ=255)

TGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAA  >  NC_000913/709790‑709969

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