BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,038,965	T→C	V154A (GTC→GCC)	b1972 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,038,965	0	T	C	99.7%	1002.6 / ‑5.3	311	V154A (GTC→GCC)	b1972	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (151/159); total (151/160)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 319 total aligned reads displayed.

GCTGGGTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTG  >  NC_000913/2038875‑2039054

gCTGGGTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCt                                                                                         <  2:7100650/93‑1 (MQ=255)

gCTGGGTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCt                                                                                         <  1:1622731/93‑1 (MQ=255)

 cTGGGTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTa                                                                                        <  1:8095770/93‑1 (MQ=255)

   gggTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTAtt                                                                                      >  1:5356564/1‑92 (MQ=255)

   gggTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATc                                                                                      >  1:6185853/1‑93 (MQ=255)

   gggTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATc                                                                                      <  2:4442217/93‑1 (MQ=255)

    ggTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCt                                                                                     >  2:4381288/1‑93 (MQ=255)

       aaTTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGt                                                                                  >  1:8915611/1‑93 (MQ=255)

        aTTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTc                                                                                 <  1:1913839/93‑1 (MQ=255)

        aTTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTc                                                                                 <  2:5868862/93‑1 (MQ=255)

         ttttGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTcg                                                                                >  1:3890305/1‑93 (MQ=255)

          tttGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTcgc                                                                               <  1:3192116/93‑1 (MQ=255)

           ttgcttgcTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTcgcg                                                                              <  1:9219708/93‑1 (MQ=255)

           tcacttgcttTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTcgcg                                                                              <  1:3434334/90‑1 (MQ=255)

            tgcttgctTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGa                                                                             <  2:7833023/93‑1 (MQ=255)

            tgcttgctTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGa                                                                             <  2:2053607/93‑1 (MQ=255)

            tgcttgctTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGa                                                                             <  1:8001608/93‑1 (MQ=255)

             gcttgcttTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGAt                                                                            <  2:9564846/93‑1 (MQ=255)

              cttgcttTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGAtt                                                                           >  1:4532671/1‑93 (MQ=255)

              cttgcttTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGAtt                                                                           <  2:4989888/93‑1 (MQ=255)

                tacttTAGCGTTCACTTCACCCCAGGCGATGCAGCCAAAACTCTGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCt                                                                         <  1:1476870/91‑1 (MQ=38)

                 gcttTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTg                                                                        >  1:9047082/1‑93 (MQ=255)

                 gcttTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTGCGCCTATCTTGTCGCGATTCTg                                                                        <  1:6827510/93‑1 (MQ=255)

                  cttTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTgg                                                                       <  1:6078561/93‑1 (MQ=255)

                   tttAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGc                                                                      >  2:7789335/1‑93 (MQ=255)

                   tttAGCGTTCACTTCAACCCAGGCGATGCAGAGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGc                                                                      >  2:4870841/1‑93 (MQ=255)

                    ttAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGcc                                                                     >  2:3897638/1‑93 (MQ=255)

                     tAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGccc                                                                    >  2:2755986/1‑93 (MQ=255)

                     tAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGccc                                                                    >  1:1184863/1‑93 (MQ=255)

                     tAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGccc                                                                    <  1:5230379/93‑1 (MQ=255)

                       gCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCa                                                                  <  2:337782/93‑1 (MQ=255)

                        cGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCaa                                                                 <  2:9153897/93‑1 (MQ=255)

                          ttCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATa                                                               >  2:7280810/1‑93 (MQ=255)

                          ttCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATa                                                               >  2:1830121/1‑93 (MQ=255)

                          ttCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATa                                                               >  2:3515931/1‑93 (MQ=255)

                          ttCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATa                                                               <  1:6255746/93‑1 (MQ=255)

                           tCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATAc                                                              >  1:8612185/1‑93 (MQ=255)

                             aCTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGGATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACAt                                                            <  1:462463/93‑1 (MQ=255)

                              cTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACAtt                                                           <  2:6009247/93‑1 (MQ=255)

                              cTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACAtt                                                           <  2:9658724/93‑1 (MQ=255)

                               ttCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACAGTa                                                          >  2:1673192/1‑93 (MQ=255)

                               cgcAACCTAGGCGATGCAGCCAAACCTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTa                                                          <  2:3178069/91‑1 (MQ=38)

                                tCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTAt                                                         >  2:6660764/1‑93 (MQ=255)

                                   aCCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCtg                                                      <  1:2594238/93‑1 (MQ=255)

                                   aCCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCtg                                                      >  2:5056338/1‑93 (MQ=255)

                                    cccAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCtgt                                                     >  1:4156230/1‑93 (MQ=255)

                                      cAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTgg                                                   <  2:4904866/93‑1 (MQ=255)

                                      cAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTgg                                                   >  1:7551513/1‑93 (MQ=255)

                                        ggCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTc                                                 <  1:2961738/93‑1 (MQ=255)

                                        ggCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTc                                                 <  2:5616336/93‑1 (MQ=255)

                                        ggCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGACTCTGGCCCCAATACATTATCTGTGGTc                                                 >  1:2916367/1‑93 (MQ=255)

                                         gCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTcc                                                <  1:9462612/93‑1 (MQ=255)

                                            aTGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTg                                             >  1:9956451/1‑93 (MQ=255)

                                            aTGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTg                                             >  2:8744193/1‑93 (MQ=255)

                                            aTGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTg                                             <  2:3586447/93‑1 (MQ=255)

                                              gCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTTTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCTTGaa                                           <  2:1943527/93‑1 (MQ=255)

                                              gCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGaa                                           <  1:2447621/93‑1 (MQ=255)

                                              gCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGaa                                           <  1:5284863/93‑1 (MQ=255)

                                              gCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGaa                                           >  1:4414997/1‑93 (MQ=255)

                                                 gCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGat                                        <  1:9148686/93‑1 (MQ=255)

                                                     aaaCTCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGattatt                                    <  2:4767650/93‑1 (MQ=255)

                                                         tCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCAc                                <  1:5520464/93‑1 (MQ=255)

                                                         tCGGCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCAc                                >  1:7213424/1‑93 (MQ=255)

                                                            gCAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGc                             <  2:4148515/93‑1 (MQ=255)

                                                             cAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCa                            <  1:6709119/93‑1 (MQ=255)

                                                             cAAGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCa                            >  2:8697195/1‑93 (MQ=255)

                                                               aGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGc                          >  1:7903159/1‑93 (MQ=255)

                                                               aGCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGc                          <  2:6504351/93‑1 (MQ=255)

                                                                gCACTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGcc                         >  1:3022784/1‑93 (MQ=255)

                                                                  ccgggCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGc                       <  1:418936/90‑1 (MQ=255)

                                                                  aCTGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGc                       >  1:7018899/1‑93 (MQ=255)

                                                                    tGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTc                     <  2:5722250/93‑1 (MQ=255)

                                                                    tGGCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTc                     <  2:1540518/93‑1 (MQ=255)

                                                                     ggCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCa                    >  1:3555559/1‑93 (MQ=255)

                                                                      gCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCAt                   <  2:8609111/93‑1 (MQ=255)

                                                                      gCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCAt                   <  1:10043852/93‑1 (MQ=255)

                                                                      gCAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCAt                   <  1:133821/93‑1 (MQ=255)

                                                                      gCAACAGTTGCATAATTACGCCTATCTTGTCGCGATTCTGGCCACAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCAt                   <  2:8977784/93‑1 (MQ=255)

                                                                       cAACAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATc                  >  2:1168077/1‑93 (MQ=255)

                                                                        aaCAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCt                 >  2:3944019/1‑93 (MQ=255)

                                                                        aaCAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCt                 >  1:8567374/1‑93 (MQ=255)

                                                                         aCAGTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTa                >  2:2158839/1‑93 (MQ=255)

                                                                            gTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgc             >  1:8518193/1‑93 (MQ=255)

                                                                            gTTGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgc             >  2:3014630/1‑93 (MQ=255)

                                                                             ttGCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgct            <  2:8894955/93‑1 (MQ=255)

                                                                               ggaTAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgg          <  2:1405606/91‑1 (MQ=255)

                                                                               gCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgg          <  2:4264320/93‑1 (MQ=255)

                                                                               gCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgg          >  1:3777895/1‑93 (MQ=255)

                                                                               gCATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgg          <  1:8798858/93‑1 (MQ=255)

                                                                               gCATAAGTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgg          <  1:4745790/93‑1 (MQ=255)

                                                                                cATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctggg         <  1:2275936/93‑1 (MQ=255)

                                                                                cATAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctggg         <  2:9370961/93‑1 (MQ=255)

                                                                                  tAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgggct       <  2:7990072/93‑1 (MQ=255)

                                                                                  tAATTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgggct       >  1:1841763/1‑93 (MQ=255)

                                                                                   aaTTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgggctg      <  2:680264/93‑1 (MQ=255)

                                                                                    aTTTCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACgctgggctgg     >  1:9907745/1‑93 (MQ=255)

                                                                                     tttCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGggctggc    >  1:3642268/1‑93 (MQ=255)

                                                                                     tttCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGggctggc    >  2:9719202/1‑93 (MQ=255)

                                                                                      ttCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGggctggct   >  1:2762277/1‑93 (MQ=255)

                                                                                      ttCGCCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGggctggct   <  1:3087895/93‑1 (MQ=255)

GCTGGGTAATTTTGCTTGCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTG  >  NC_000913/2038875‑2039054

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