BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,020,045	0	C	T	65.0%	475.8 / inf	411	A54V (GCG→GTG)	b3838	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (73/71); new base T (133/134); total (206/205)
Fisher's exact test for biased strand distribution p-value = 9.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-03
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 444 total aligned reads displayed.

ACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGCACACCA  >  NC_000913/4019953‑4020137

acTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGc                                                                                              <  2:3110216/93‑1 (MQ=255)

  tGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAgtgt                                                                                            >  2:5992413/1‑93 (MQ=255)

  tGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAgtgt                                                                                            >  2:9567884/1‑93 (MQ=255)

  tGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAgtgt                                                                                            <  1:2371950/93‑1 (MQ=255)

  tGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGATGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAgtgt                                                                                            >  2:893978/1‑93 (MQ=255)

     cccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCga                                                                                         <  1:9750801/93‑1 (MQ=255)

       cAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgatg                                                                                       >  1:1879154/1‑93 (MQ=255)

        aGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCgatgg                                                                                      >  1:8788287/1‑93 (MQ=255)

         ggAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgatgga                                                                                     >  1:4355773/1‑93 (MQ=255)

           aGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgatggatg                                                                                   >  1:5515242/1‑93 (MQ=255)

           aGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgatggatg                                                                                   <  1:3779517/93‑1 (MQ=255)

            gTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGa                                                                                  >  2:8241854/1‑93 (MQ=255)

                    tCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGcc                                                                          <  1:6409403/93‑1 (MQ=255)

                    tCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGcc                                                                          >  2:1309706/1‑93 (MQ=255)

                     ccAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCa                                                                         <  1:7167079/93‑1 (MQ=255)

                     ccAGGAGTTTAAGGACAGTCGGAAAAAGGATGAAAAGGCAAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCAATGGATGAACTACGCCa                                                                         <  2:4882875/93‑1 (MQ=255)

                      cAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAg                                                                        <  2:8091767/93‑1 (MQ=255)

                      cAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAg                                                                        <  2:9861121/93‑1 (MQ=255)

                      cAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAg                                                                        >  2:6810030/1‑93 (MQ=255)

                      cAGGAGTTTCAGGACAGTCTGAAAAAGGGTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGAGGGATGAAGTACGCCAg                                                                        >  1:5648037/1‑93 (MQ=255)

                         gAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGcc                                                                     >  2:9773809/1‑93 (MQ=255)

                            tttCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCcgcg                                                                  >  2:7132007/1‑93 (MQ=255)

                             ttCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCgg                                                                 >  1:2233280/1‑93 (MQ=255)

                             ttCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCgg                                                                 >  2:8486832/1‑93 (MQ=255)

                              tCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGa                                                                >  2:2956928/1‑93 (MQ=255)

                               cAGGACAGTCTGAAAAAGGTTGAAAAGGCGGGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAg                                                               >  1:7426873/1‑93 (MQ=255)

                               cAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGANGAACTACGCCAGGCCGCGGAg                                                               >  1:692898/1‑93 (MQ=255)

                                aGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAgc                                                              >  2:5101988/1‑92 (MQ=255)

                                 ggACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTc                                                             >  2:6544086/1‑93 (MQ=255)

                                 ggACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTc                                                             >  2:8586811/1‑93 (MQ=255)

                                 ggACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTc                                                             >  1:4255688/1‑93 (MQ=255)

                                 ggACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgatgga                                                                                     <  1:9797603/69‑1 (MQ=255)

                                 ggACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTc                                                             >  2:8781155/1‑93 (MQ=255)

                                 agACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTc                                                             <  1:8562468/92‑1 (MQ=255)

                                   aCAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGa                                                           >  2:7304695/1‑93 (MQ=255)

                                   aCAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGa                                                           >  1:6700106/1‑93 (MQ=255)

                                   aCAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGa                                                           >  2:5343601/1‑93 (MQ=255)

                                    cAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGAt                                                          >  2:3671311/1‑93 (MQ=255)

                                    cAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGAt                                                          <  1:9120255/93‑1 (MQ=255)

                                    cAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGAt                                                          >  2:3302610/1‑93 (MQ=255)

                                      gTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGa                                                        <  2:3930035/93‑1 (MQ=255)

                                      gTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGa                                                        <  2:46011/93‑1 (MQ=255)

                                        cTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAg                                                      <  2:4017482/93‑1 (MQ=255)

                                         tGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGc                                                     <  1:6944989/93‑1 (MQ=255)

                                         tGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGc                                                     >  1:5784632/1‑93 (MQ=255)

                                           aaaaaGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGt                                                   >  2:1948009/1‑93 (MQ=255)

                                            aaaaGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGtt                                                  <  2:835562/93‑1 (MQ=255)

                                            aaaaGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGtt                                                  >  2:5755630/1‑93 (MQ=255)

                                            aaaaGGTTAAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGtt                                                  <  2:7486674/93‑1 (MQ=255)

                                             aaaGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCATGCCGCGGAGTCGATGAAGCGTTc                                                 >  2:5595397/1‑93 (MQ=255)

                                              aaGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTcc                                                <  1:3632703/93‑1 (MQ=255)

                                               aGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCt                                               <  2:2701085/93‑1 (MQ=255)

                                                ggTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTa                                              <  1:100192/93‑1 (MQ=255)

                                                 gTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTAc                                             >  2:2408818/1‑93 (MQ=255)

                                                 gTTGAAAAGCGGAGCCTCACAAACCTGGCGCCCGAACTG‑AAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACg                                            <  1:1324932/93‑1 (MQ=255)

                                                   tGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGt                                           >  1:6759959/1‑93 (MQ=255)

                                                    gAAAAGCCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGtt                                          <  1:49704/93‑1 (MQ=255)

                                                      aaaGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGc                                        >  1:1009532/1‑93 (MQ=255)

                                                       aaGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCa                                       >  1:1672508/1‑93 (MQ=255)

                                                         ggCGAGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCaaa                                     >  1:6842666/1‑93 (MQ=255)

                                                             aGCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGAt                                 >  1:6071076/1‑93 (MQ=255)

                                                             aGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGAt                                 <  2:1659269/93‑1 (MQ=255)

                                                              gCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATc                                <  2:6313181/93‑1 (MQ=255)

                                                              gCCTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGCCGATGAAGCGTTCCTACGTTGCAAACGATc                                >  2:6219029/1‑93 (MQ=255)

                                                               ccTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATcc                               <  2:1749873/93‑1 (MQ=255)

                                                                cTCACTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCt                              <  2:3024573/93‑1 (MQ=255)

                                                                   aCTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaa                           <  1:2040989/93‑1 (MQ=255)

                                                                    cTAACCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaaa                          <  1:4540118/93‑1 (MQ=255)

                                                                    cTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTCCCAACGCTCCTGaaa                          >  1:7820917/1‑93 (MQ=255)

                                                                     tAACCTGACGCCCCAACGGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaaaa                         <  1:1426338/93‑1 (MQ=255)

                                                                      aaCCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAg                        <  2:2236071/93‑1 (MQ=255)

                                                                      aaCCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAg                        >  1:4748907/1‑93 (MQ=255)

                                                                       aCCTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAgg                       >  2:1534280/1‑93 (MQ=255)

                                                                        ccTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgc                      >  1:5165388/1‑93 (MQ=255)

                                                                        ccTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgc                      <  1:4172240/93‑1 (MQ=255)

                                                                         cTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcg                     >  1:2663566/1‑93 (MQ=255)

                                                                         cTGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcg                     <  2:350383/93‑1 (MQ=255)

                                                                          tGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcga                    <  2:5537610/93‑1 (MQ=255)

                                                                          tGACGCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcga                    <  2:7065463/93‑1 (MQ=255)

                                                                            ccGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgagc                  <  2:2092342/92‑1 (MQ=255)

                                                                            aCGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgagc                  <  1:9527686/93‑1 (MQ=255)

                                                                              gCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgagcga                <  2:1126146/93‑1 (MQ=255)

                                                                              gCCCGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGgcgagcga                <  1:2995631/93‑1 (MQ=255)

                                                                               cccGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGAt               <  2:5215644/93‑1 (MQ=255)

                                                                               cccGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGAt               >  2:8167992/1‑93 (MQ=255)

                                                                                ccGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATg              >  2:9187259/1‑93 (MQ=255)

                                                                                ccGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATg              <  2:3347840/93‑1 (MQ=255)

                                                                                 cGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGa             <  1:8817403/93‑1 (MQ=255)

                                                                                 cGAACTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGa             <  1:3330412/93‑1 (MQ=255)

                                                                                 cGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGa             >  1:5943964/1‑93 (MQ=255)

                                                                                   aaCTGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAg           <  2:4434835/93‑1 (MQ=255)

                                                                                      tGAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAgcgc        <  1:9322333/93‑1 (MQ=255)

                                                                                       gAAAGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGca       >  2:675765/1‑93 (MQ=255)

                                                                                        aaaGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcac      <  1:6635466/93‑1 (MQ=255)

                                                                                        aaaGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcac      >  2:4683723/1‑93 (MQ=255)

                                                                                        aaaGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcac      <  1:9906838/93‑1 (MQ=255)

                                                                                         aaGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcaca     <  1:4052136/93‑1 (MQ=255)

                                                                                         aaGTGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcaca     <  2:5890705/93‑1 (MQ=255)

                                                                                         aaGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGcaca     >  2:9658911/1‑93 (MQ=255)

                                                                                          aGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGCacac    <  1:3675487/93‑1 (MQ=255)

ACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAAGGCGAGCGATGAAGCGCACACCA  >  NC_000913/4019953‑4020137

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