BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	542,052	A→G	Q314R (CAG→CGG)	ybbZ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	542,052	0	A	G	82.7%	452.6 / 97.0	214	Q314R (CAG→CGG)	ybbZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (19/18); new base G (100/77); total (119/95)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

Only 100 of 215 total aligned reads displayed.

GGGCGCAGCACTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCA  >  NC_000913/541960‑542141

gggCGCAGCACTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCa                                                                                           >  1:263249/1‑93 (MQ=255)

   cgcAGCACTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGtt                                                                                        <  2:4756530/93‑1 (MQ=255)

    gcagcaCTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGttt                                                                                       <  2:5586283/93‑1 (MQ=255)

      agcaCTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTaa                                                                                     >  2:2465008/1‑93 (MQ=255)

       gcaCTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAAt                                                                                    <  1:6535503/93‑1 (MQ=255)

         aCTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGt                                                                                  >  1:1499165/1‑93 (MQ=255)

          cTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTa                                                                                 >  2:4963781/1‑93 (MQ=255)

          cTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTa                                                                                 >  2:342676/1‑93 (MQ=255)

           tGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTAc                                                                                >  2:2791850/1‑93 (MQ=255)

           tGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTAc                                                                                <  1:4643767/93‑1 (MQ=255)

           tGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTAc                                                                                >  1:2586237/1‑93 (MQ=255)

            ggTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTAcc                                                                               >  2:5641162/1‑93 (MQ=255)

            ggTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTAcc                                                                               <  2:7345305/93‑1 (MQ=255)

            ggTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTAcc                                                                               <  1:2292811/93‑1 (MQ=255)

             gTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCg                                                                              <  2:5568009/93‑1 (MQ=255)

              tGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCgg                                                                             >  2:6735147/1‑93 (MQ=255)

               gATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGt                                                                            <  2:4791042/93‑1 (MQ=255)

               gATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGt                                                                            <  1:2966618/93‑1 (MQ=255)

                 ttACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGa                                                                          <  1:315164/93‑1 (MQ=255)

                   aCCGGGGAAGGGCGCGTCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGAtt                                                                        >  1:4501232/1‑93 (MQ=255)

                   aCCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGAtt                                                                        <  2:4187435/93‑1 (MQ=255)

                   aCCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGAtt                                                                        <  1:6596769/93‑1 (MQ=255)

                   aCCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGAtt                                                                        >  1:570931/1‑93 (MQ=255)

                     cGGGGAAGGGCGCATCGACTCGCAAACGGGAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGCGATTgg                                                                      >  1:2376918/1‑93 (MQ=255)

                     cGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTgg                                                                      <  2:4533380/93‑1 (MQ=255)

                      ggggAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTggg                                                                     >  1:3243345/1‑93 (MQ=255)

                      ggggAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGGGTGGCGTCGGGGGCGAAGCGGGTTAATGTACCGGTGATTggg                                                                     >  1:4251577/1‑93 (MQ=255)

                       gggAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGa                                                                    <  2:6383608/93‑1 (MQ=255)

                        ggAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGAt                                                                   >  1:458216/1‑93 (MQ=255)

                        ggAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGAt                                                                   >  1:7155783/1‑93 (MQ=255)

                          aaGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTg                                                                 >  1:4213181/1‑93 (MQ=255)

                          aaGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTg                                                                 >  2:1516726/1‑93 (MQ=255)

                             ggCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTg                                                              >  2:4361734/1‑93 (MQ=255)

                                gcATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCg                                                           >  2:5089347/1‑93 (MQ=255)

                                gcATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCg                                                           >  2:1906282/1‑93 (MQ=255)

                                   tCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTAt                                                        >  1:3986395/1‑93 (MQ=255)

                                   tCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTAt                                                        >  2:4464565/1‑93 (MQ=255)

                                   tCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTAt                                                        <  2:6579100/93‑1 (MQ=255)

                                   tCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTAt                                                        >  1:6114517/1‑93 (MQ=255)

                                    ggACTCGCAAACGGCAGGCGGTAAAGCGCTGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTAtt                                                       <  2:1204157/92‑1 (MQ=255)

                                     gACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTg                                                      >  2:5500571/1‑93 (MQ=255)

                                     gACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTg                                                      <  1:3898644/93‑1 (MQ=255)

                                      aCTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGTAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTgg                                                     <  2:1069531/93‑1 (MQ=255)

                                      aCTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTgg                                                     <  2:4613660/93‑1 (MQ=255)

                                      aCTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTgg                                                     >  1:2828505/1‑93 (MQ=255)

                                      aCTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTgg                                                     >  2:7113190/1‑93 (MQ=255)

                                       cTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTggg                                                    >  2:3054960/1‑93 (MQ=255)

                                       cTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTggg                                                    >  2:7025851/1‑93 (MQ=255)

                                        tCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGt                                                   >  1:1086501/1‑93 (MQ=255)

                                        tCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGt                                                   <  2:3912034/93‑1 (MQ=255)

                                          gCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGa                                                 <  1:6964290/93‑1 (MQ=255)

                                            aaaCGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATg                                               <  2:165626/93‑1 (MQ=255)

                                             aaCGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATgg                                              <  1:4858256/93‑1 (MQ=255)

                                              aCGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGc                                             >  1:1876285/1‑93 (MQ=255)

                                                 gCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTg                                          <  2:6771706/93‑1 (MQ=255)

                                                 gCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTg                                          >  2:3255761/1‑93 (MQ=255)

                                                   aGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGa                                        >  2:317522/1‑93 (MQ=255)

                                                   aGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGa                                        <  2:7144605/93‑1 (MQ=255)

                                                      cGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAgt                                     <  2:1778198/93‑1 (MQ=255)

                                                       ggTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAgtg                                    >  1:682914/1‑93 (MQ=255)

                                                       ggTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAgnn                                    >  2:2757335/1‑91 (MQ=255)

                                                          aaaGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAgtggtg                                 <  1:3028320/93‑1 (MQ=255)

                                                          aaaGCGCCGCTGGGTGTGGCGTCGGCGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAgtggtg                                 >  1:6795390/1‑93 (MQ=255)

                                                           aaGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGc                                >  1:7311473/1‑93 (MQ=255)

                                                              cgccgcTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCAcc                             <  2:584711/93‑1 (MQ=255)

                                                               gccgcTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCa                            >  1:1414753/1‑93 (MQ=255)

                                                                 cgcTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGt                          >  2:5772383/1‑93 (MQ=255)

                                                                 cgcTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGt                          >  2:5228043/1‑93 (MQ=255)

                                                                  gcTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTa                         >  2:75293/1‑93 (MQ=255)

                                                                  gcTGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTa                         <  1:5586/93‑1 (MQ=255)

                                                                    tGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACg                       >  2:494558/1‑93 (MQ=255)

                                                                    tGGGTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACg                       >  2:511301/1‑93 (MQ=255)

                                                                    tGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACg                       >  2:882074/1‑93 (MQ=255)

                                                                     gggTGTGGCGTCGGTTGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACgg                      >  2:3712154/1‑93 (MQ=255)

                                                                      ggTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGc                     >  1:6320051/1‑93 (MQ=255)

                                                                      ggTGTGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGc                     <  2:5416002/93‑1 (MQ=255)

                                                                      ggTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGc                     <  1:1333134/93‑1 (MQ=255)

                                                                       gtgtGGCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCa                    <  1:2259298/93‑1 (MQ=255)

                                                                         gtgGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCAtt                  >  1:3482960/1‑93 (MQ=255)

                                                                           ggCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGa                <  1:6365010/93‑1 (MQ=255)

                                                                           ggCGTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGa                >  2:6195402/1‑93 (MQ=255)

                                                                           ggCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGa                <  1:366073/93‑1 (MQ=255)

                                                                            gCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGAc               <  2:5972930/93‑1 (MQ=255)

                                                                              gTCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGAcgc             >  1:2484106/1‑93 (MQ=255)

                                                                               tCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGAcgcg            >  1:6907208/1‑93 (MQ=255)

                                                                               tCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGAcgcg            >  1:6489824/1‑93 (MQ=255)

                                                                               tCGGTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGAcgcg            >  2:2338932/1‑93 (MQ=255)

                                                                                  gTGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTa         <  1:4682971/93‑1 (MQ=255)

                                                                                   tGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTAt        <  1:982573/93‑1 (MQ=255)

                                                                                   tGGCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTAt        >  1:1606091/1‑93 (MQ=255)

                                                                                   tGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTAt        <  1:5528511/93‑1 (MQ=255)

                                                                                    ggggAAGCGGTGTAATGTACCGATGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTAtt       <  2:5186923/90‑1 (MQ=255)

                                                                                     gCGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTc      <  2:2646990/93‑1 (MQ=255)

                                                                                     gCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTc      <  1:74818/93‑1 (MQ=255)

                                                                                      cGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCa     >  2:5516517/1‑93 (MQ=255)

                                                                                      cGAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCa     <  2:4571419/93‑1 (MQ=255)

                                                                                       gAAGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAg    >  2:3989327/1‑93 (MQ=255)

                                                                                       gAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAg    >  1:1276887/1‑93 (MQ=255)

                                                                                        aaGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGc   <  1:4299118/93‑1 (MQ=255)

                                                                                         aGCGGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCa  >  1:2323213/1‑93 (MQ=255)

GGGCGCAGCACTGGTGATTACCGGGGAAGGGCGCATCGACTCGCAAACGGCAGGCGGTAAAGCGCCGCTGGGTGTGGCGTCGGTGGCGAAGCAGTTTAATGTACCGGTGATTGGGATTGCTGGCGTATTGGGTGATGGCGTGGAAGTGGTGCACCAGTACGGCATTGACGCGGTATTCAGCA  >  NC_000913/541960‑542141

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