BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,562,351	G→A	G118S (GGC→AGC)	b2449 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,562,351	0	G	A	99.5%	669.4 / NA	189	G118S (GGC→AGC)	b2449	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (88/100); minor base C (0/1); 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 = 1.00e+00

Only 100 of 191 total aligned reads displayed.

ATCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGG  >  NC_000913/2562261‑2562440

aTCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGc                                                                                         >  2:2050033/1‑93 (MQ=255)

aTCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGc                                                                                         >  1:1545385/1‑93 (MQ=255)

aTCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGc                                                                                         >  1:1760731/1‑93 (MQ=255)

aTCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGc                                                                                         <  1:2861548/93‑1 (MQ=255)

 tCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCt                                                                                        <  1:2486827/93‑1 (MQ=255)

  tcTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCtt                                                                                       <  1:4069781/92‑1 (MQ=255)

  tcTGATAAAATAACATTTCAGCAATGAAACTGACTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCtt                                                                                       <  1:4147852/92‑1 (MQ=255)

  ccTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCtt                                                                                       <  1:457833/93‑1 (MQ=255)

  ccTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCtt                                                                                       <  2:3135422/93‑1 (MQ=255)

   cTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCttt                                                                                      <  1:2739781/93‑1 (MQ=255)

    tGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTg                                                                                     <  2:3430491/93‑1 (MQ=255)

     ggTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGt                                                                                    >  2:1122667/1‑93 (MQ=255)

      gTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGtt                                                                                   <  1:3431273/93‑1 (MQ=255)

       tAAAATGCCATTTCAGCAAAGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGttt                                                                                  <  2:4219948/93‑1 (MQ=255)

        aaaaTACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGtttt                                                                                 <  2:1916929/93‑1 (MQ=255)

        aaaaTACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGtttt                                                                                 >  2:1182316/1‑93 (MQ=255)

          aaTACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTgg                                                                               <  2:1270890/93‑1 (MQ=255)

          aaTACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTgg                                                                               >  2:3382255/1‑93 (MQ=255)

           aTACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggt                                                                              <  2:4188879/93‑1 (MQ=255)

            tACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggta                                                                             <  1:3126003/93‑1 (MQ=255)

            tACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggta                                                                             <  2:2416361/93‑1 (MQ=255)

             aCCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggtag                                                                            >  2:273393/1‑93 (MQ=255)

                aTTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTggtaggta                                                                         >  2:2476042/1‑93 (MQ=255)

                 tttCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGgtaggtag                                                                        <  2:2625576/93‑1 (MQ=255)

                 tttCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGgtaggtag                                                                        >  1:2810857/1‑93 (MQ=255)

                    cAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGt                                                                     <  2:284595/93‑1 (MQ=255)

                     aGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTa                                                                    <  2:3296491/93‑1 (MQ=255)

                      gCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAt                                                                   >  2:1295945/1‑93 (MQ=255)

                      gCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGGAt                                                                   >  2:4290200/1‑93 (MQ=255)

                       cAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAtt                                                                  <  1:1515359/93‑1 (MQ=255)

                        aaTGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTAttt                                                                 >  1:654956/1‑93 (MQ=255)

                         aTGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTg                                                                >  2:448354/1‑93 (MQ=255)

                           gAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGaa                                                              >  2:360410/1‑93 (MQ=255)

                            aaaCTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGaaa                                                             >  2:2162116/1‑93 (MQ=255)

                             aaCTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGaaaa                                                            <  2:1592330/93‑1 (MQ=255)

                             aaCTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGaaaa                                                            >  2:4183246/1‑93 (MQ=255)

                                tGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACg                                                         >  1:4241696/1‑93 (MQ=255)

                                tGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACg                                                         >  2:3114281/1‑93 (MQ=255)

                                tGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGCATTTGAAAAACg                                                         >  2:382364/1‑93 (MQ=255)

                                 gAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAACAACGt                                                        >  1:352108/1‑93 (MQ=255)

                                 gAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGt                                                        <  2:1008427/93‑1 (MQ=255)

                                  aGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTGTGGTAGGTAGAGTATTTGAAAAACGTa                                                       <  1:3476657/93‑1 (MQ=255)

                                    ttCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAat                                                     <  2:501172/93‑1 (MQ=255)

                                     tCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAatt                                                    <  2:2092491/93‑1 (MQ=255)

                                       aaaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattga                                                  <  2:1377501/93‑1 (MQ=255)

                                        aaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgat                                                 >  2:324327/1‑93 (MQ=255)

                                        aaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgat                                                 >  2:4310970/1‑93 (MQ=255)

                                        aaCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgat                                                 <  2:3901614/93‑1 (MQ=255)

                                         aCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgatt                                                <  1:4698970/93‑1 (MQ=255)

                                         aCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgatt                                                >  1:4640708/1‑93 (MQ=255)

                                         aCTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgatt                                                <  1:2489438/93‑1 (MQ=255)

                                          cTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgattg                                               >  2:170372/1‑93 (MQ=255)

                                          cTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgattg                                               >  1:964529/1‑93 (MQ=255)

                                          cTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAattgattg                                               <  1:4130467/93‑1 (MQ=255)

                                           ttAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGc                                              <  2:3038566/93‑1 (MQ=255)

                                           ttAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGc                                              <  1:518309/93‑1 (MQ=255)

                                            tAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCg                                             >  2:747310/1‑93 (MQ=255)

                                             aaaGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGCAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGa                                            >  1:2833855/1‑93 (MQ=255)

                                              aaGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGaa                                           <  1:3688761/93‑1 (MQ=255)

                                               aGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAAt                                          <  2:926606/93‑1 (MQ=255)

                                               aGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAAt                                          >  1:4286825/1‑93 (MQ=255)

                                                gCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATa                                         <  2:857305/93‑1 (MQ=255)

                                                gCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATa                                         <  2:4136121/93‑1 (MQ=255)

                                                  ggCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAga                                       >  1:2913603/1‑92 (MQ=255)

                                                  ggCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGc                                       <  1:778959/93‑1 (MQ=255)

                                                  ggCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGc                                       <  1:627245/93‑1 (MQ=255)

                                                  ggCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGc                                       <  1:221638/93‑1 (MQ=255)

                                                   gCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGca                                      >  1:4100668/1‑93 (MQ=255)

                                                         gTCAGAGAATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACt                                >  2:1734334/1‑93 (MQ=255)

                                                             gagaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaac                            <  1:3363953/93‑1 (MQ=255)

                                                             gagaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaac                            >  1:4209955/1‑93 (MQ=255)

                                                             gagaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaac                            >  2:4186605/1‑93 (MQ=255)

                                                             gagaATGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGCATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAaac                            >  1:1236976/1‑93 (MQ=255)

                                                                 aTGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAAcc                        >  2:2871239/1‑93 (MQ=255)

                                                                 aTGCGCCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAAcc                        <  1:3720688/93‑1 (MQ=255)

                                                                   gcgcCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCaa                      <  2:260298/93‑1 (MQ=255)

                                                                    cgcCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAAt                     >  1:3150571/1‑93 (MQ=255)

                                                                    cgcCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAAt                     <  2:3912527/93‑1 (MQ=255)

                                                                    cgcCAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAAt                     <  2:2925735/93‑1 (MQ=255)

                                                                      ccAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATac                   >  1:4461523/1‑93 (MQ=255)

                                                                      ccAACGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATac                   >  1:2321863/1‑93 (MQ=255)

                                                                        aaCGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATacac                 >  2:3251976/1‑93 (MQ=255)

                                                                        aaCGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATacac                 >  1:3824722/1‑93 (MQ=255)

                                                                        aaCGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATacac                 <  2:4064148/93‑1 (MQ=255)

                                                                         aCGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAGCGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACAcc                <  1:1072822/93‑1 (MQ=255)

                                                                         aCGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACAcc                >  1:1589100/1‑93 (MQ=255)

                                                                         aCGTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACAcc                <  1:763272/93‑1 (MQ=255)

                                                                          cGTAGCAGAAAACAACAGCTTTGTTTTGGTAGATAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACAccc               <  1:3318010/93‑1 (MQ=255)

                                                                           gTAGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCg              <  2:629640/93‑1 (MQ=255)

                                                                             aGCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGc            >  1:99843/1‑93 (MQ=255)

                                                                              gCAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGcc           >  1:3642718/1‑93 (MQ=255)

                                                                               cAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGccc          <  2:2694678/93‑1 (MQ=255)

                                                                               cAGAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGccc          >  2:4354168/1‑93 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        <  1:2994554/93‑1 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        <  1:2174813/93‑1 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        >  1:2132376/1‑93 (MQ=255)

                                                                                 gAAAACAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCt        <  1:1925537/93‑1 (MQ=255)

                                                                                  aaaaCAACAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTc       >  2:1464837/1‑93 (MQ=255)

                                                                                    aacaacAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGc     >  1:3662589/1‑93 (MQ=255)

                                                                                       aacAGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCgg  >  1:558797/1‑93 (MQ=255)

ATCCTGGTAAAATACCATTTCAGCAATGAAACTGAGTTCAAACTTAAAGCGGCAGCGGTCAGAGAATGCGCCAACGTAGCAGAAAACAACGGCTTTGTTTTGGTAGGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGG  >  NC_000913/2562261‑2562440

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