BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,802,603	T→C	G495G (GGT→GGC)	b1721m →	ankyrin repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,802,603	0	T	C	99.5%	591.7 / ‑4.9	189	G495G (GGT→GGC)	b1721m	ankyrin repeat protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (88/100); total (88/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 191 total aligned reads displayed.

TTTCTGAATGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGTGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTAAATA  >  NC_000913/1802514‑1802692

  tCTGAATGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGcc                                                                                      >  1:1635371/1‑93 (MQ=255)

  tCTGAATGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGcc                                                                                      <  1:355191/93‑1 (MQ=255)

   cTGAATGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCa                                                                                     >  2:2372769/1‑93 (MQ=255)

    tGAATGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAg                                                                                    <  1:2824149/93‑1 (MQ=255)

      aaTGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGa                                                                                  >  1:1292817/1‑93 (MQ=255)

      aaTGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGa                                                                                  <  2:2775390/93‑1 (MQ=255)

      aaTGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGa                                                                                  <  2:5610079/93‑1 (MQ=255)

      aaTGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGGATTTTAGCCATGAAAATACATGATGGCGAGCCAGGa                                                                                  >  1:7774298/1‑93 (MQ=255)

          cATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATtat                                                                              <  1:7314864/93‑1 (MQ=255)

           aTTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATtata                                                                             >  2:4402283/1‑93 (MQ=255)

           aTTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATtata                                                                             <  2:4426755/93‑1 (MQ=255)

           aTTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATtata                                                                             <  1:6509156/93‑1 (MQ=255)

           aTTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATtata                                                                             <  2:6009220/93‑1 (MQ=255)

           aTTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGGGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCCGGATtata                                                                             >  1:4062924/1‑93 (MQ=255)

            ttAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATAc                                                                            <  2:4344646/93‑1 (MQ=255)

              aaCCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATAcgc                                                                          <  2:5035164/93‑1 (MQ=255)

                ccACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATAcgccg                                                                        >  2:3361323/1‑93 (MQ=255)

                 cACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATAcgccgc                                                                       <  1:2514461/93‑1 (MQ=255)

                 cACTCAAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATAcgccgc                                                                       <  2:7679370/93‑1 (MQ=255)

                  aCTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCa                                                                      >  2:7395552/1‑93 (MQ=255)

                   cTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCaa                                                                     <  2:3462292/93‑1 (MQ=255)

                   cTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCaa                                                                     <  2:6252719/93‑1 (MQ=255)

                    tCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAAt                                                                    >  2:1129211/1‑93 (MQ=255)

                    tCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAAt                                                                    <  1:3976330/93‑1 (MQ=255)

                    tCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAAt                                                                    <  1:272765/93‑1 (MQ=255)

                       aaTAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGa                                                                 <  1:2827366/93‑1 (MQ=255)

                       aaTAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGa                                                                 <  2:1235636/93‑1 (MQ=255)

                       aaTAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGa                                                                 <  2:7194395/93‑1 (MQ=255)

                         aaGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTGAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGaaa                                                               <  2:3564407/92‑1 (MQ=255)

                           gCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAaat                                                             <  1:1427277/93‑1 (MQ=255)

                             acaAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAaataa                                                           <  1:4667617/93‑1 (MQ=255)

                             acaAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAaataa                                                           <  1:3762317/93‑1 (MQ=255)

                             acaAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAaataa                                                           >  1:232109/1‑93 (MQ=255)

                              caAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAaataat                                                          <  2:3221238/93‑1 (MQ=255)

                                agagCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCa                                                        >  2:4899199/1‑93 (MQ=255)

                                   gCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCAccc                                                     <  1:713726/93‑1 (MQ=255)

                                        tAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTtgt                                                >  2:2506731/1‑93 (MQ=255)

                                         aaaCAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTtgtg                                               <  1:101456/93‑1 (MQ=255)

                                         aaaCAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTtgtg                                               >  1:7607325/1‑93 (MQ=255)

                                         aaaCAAAAAACTGGTGATGAGTAATTCAGCCATGAAATCACATGATGGCGATCCAGGACTATACGCCGCAATGGAAAATAATCACCCTCtgtg                                               >  2:7028756/1‑93 (MQ=255)

                                         aaaCAAAAAACGGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTtgtg                                               <  1:598256/93‑1 (MQ=255)

                                          aaCAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTtgtgt                                              >  1:7347681/1‑93 (MQ=255)

                                            cAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTgtgtgt                                            >  2:3702810/1‑93 (MQ=255)

                                            aaaaaaaCTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTCTACGCCGCAATGGAAAATAATCACCCTTTgtgtgt                                            <  2:7798143/92‑1 (MQ=255)

                                             aaaaaaCTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTc                                           >  1:4768193/1‑93 (MQ=255)

                                              aaaaaCTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTca                                          >  2:2928238/1‑93 (MQ=255)

                                                aaaCTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTcaca                                        <  2:7817859/93‑1 (MQ=255)

                                                  aCTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACg                                      <  2:2345230/93‑1 (MQ=255)

                                                  aCTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACg                                      >  2:6978387/1‑93 (MQ=255)

                                                   cTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACgg                                     >  2:5355888/1‑93 (MQ=255)

                                                    tGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGt                                    <  1:5884728/93‑1 (MQ=255)

                                                    tGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGt                                    >  1:1001841/1‑93 (MQ=255)

                                                    tGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGt                                    <  1:7530746/93‑1 (MQ=255)

                                                    tGGTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGt                                    >  1:436513/1‑93 (MQ=255)

                                                     ggTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGtt                                   <  2:1705707/93‑1 (MQ=255)

                                                     ggTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGtt                                   >  1:3035484/1‑93 (MQ=255)

                                                      gTGATGAGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTc                                  <  2:550159/93‑1 (MQ=255)

                                                       tgatgaGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTcc                                 <  2:3304629/93‑1 (MQ=255)

                                                        gatgaGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCct                                >  1:3673905/1‑93 (MQ=255)

                                                        gatgaGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCct                                >  1:4827512/1‑93 (MQ=255)

                                                        gatgaGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCct                                >  2:3146632/1‑93 (MQ=255)

                                                        gatgaGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCct                                >  1:5135081/1‑93 (MQ=255)

                                                         atgaGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCctc                               <  2:615530/93‑1 (MQ=255)

                                                          tgaGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCctct                              >  2:5921524/1‑93 (MQ=255)

                                                          tgaGTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCctct                              >  2:2368711/1‑93 (MQ=255)

                                                             gTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTa                           >  2:4452793/1‑93 (MQ=255)

                                                             gTATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTa                           <  2:4112974/93‑1 (MQ=255)

                                                             gTATTTAAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTa                           <  2:5036055/93‑1 (MQ=255)

                                                              tATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTaa                          >  1:619251/1‑93 (MQ=255)

                                                              tATTTTAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTaa                          >  1:6392244/1‑93 (MQ=255)

                                                                ttttAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTaaaa                        >  1:703202/1‑93 (MQ=255)

                                                                  ttAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAAtt                      <  1:5212921/93‑1 (MQ=255)

                                                                   tAGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTa                     <  2:885486/93‑1 (MQ=255)

                                                                    aGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTaa                    <  2:6365614/93‑1 (MQ=255)

                                                                    aGCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTaa                    <  1:7440298/93‑1 (MQ=255)

                                                                     gCCATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAAt                   <  2:3975781/93‑1 (MQ=255)

                                                                      ccATGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATg                  >  2:6813626/1‑93 (MQ=255)

                                                                        aTGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGc                <  1:3071492/93‑1 (MQ=255)

                                                                        aTGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATAGc                >  2:2152779/1‑93 (MQ=255)

                                                                         tGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCa               <  1:7832822/93‑1 (MQ=255)

                                                                         tGAAAATACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCa               <  2:6189723/93‑1 (MQ=255)

                                                                         tGAAAAGACATGATGGCGAGCCAGGATTATACGCCGCAAGGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCa               <  1:5031709/93‑1 (MQ=255)

                                                                           aaaaTACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATc             <  1:4402284/93‑1 (MQ=255)

                                                                           aaaaTACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATc             <  1:1994223/93‑1 (MQ=255)

                                                                           aaaaTACAGGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATc             <  2:6214172/93‑1 (MQ=255)

                                                                            aaaTACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCg            >  1:3779124/1‑93 (MQ=255)

                                                                            aaaTACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCg            <  1:7184216/93‑1 (MQ=255)

                                                                              aTACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGcc          >  2:4748373/1‑93 (MQ=255)

                                                                              aTACATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGcc          <  1:6788944/93‑1 (MQ=255)

                                                                                aCATGATGGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCtt        >  1:1807812/1‑93 (MQ=255)

                                                                                 cATGATGGCGATCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCttt       <  2:3475959/93‑1 (MQ=255)

                                                                                  atgatgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTa      <  2:7535559/93‑1 (MQ=255)

                                                                                  atgatgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTa      >  1:4926737/1‑93 (MQ=255)

                                                                                  atgatgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTa      >  1:4125467/1‑93 (MQ=255)

                                                                                    gatgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTaaa    <  1:7693118/93‑1 (MQ=255)

                                                                                     atgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTAAAt   <  2:2841101/93‑1 (MQ=255)

                                                                                     atgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTAAAt   >  1:5448139/1‑93 (MQ=255)

                                                                                      tgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTAAATa  <  2:6672289/93‑1 (MQ=255)

                                                                                      tgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTAAATa  <  1:5161566/93‑1 (MQ=255)

                                                                                      tgGCGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTAAATa  <  1:1141338/93‑1 (MQ=255)

TTTCTGAATGCATTAACCACTCCAATAGCACAAGAGCTTTTAAACAAAAAACTGGTGATGAGTATTTTAGCCATGAAAATACATGATGGTGAGCCAGGATTATACGCCGCAATGGAAAATAATCACCCTTTGTGTGTCACACGGTTCCTCTCTAAAATTAATGGCATCGCCTTTAAATA  >  NC_000913/1802514‑1802692

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