BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,289,431	T→C	G150G (GGT→GGC)	yraK →	putative fimbrial protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,289,431	0	T	C	100.0%	625.9 / NA	200	G150G (GGT→GGC)	yraK	putative fimbrial protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (101/99); total (101/99)

Only 100 of 210 total aligned reads displayed.

CGACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGTGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAT  >  NC_000913/3289341‑3289520

cGACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGa                                                                                         >  1:3730234/1‑93 (MQ=255)

cGACCAGGTTGCTCCTGTTACCCCCCGATCCACACAGATAAAAAACACCCGGGTATATCCTCCTGCAACCCAAACGGGGTTTTGCCAGGGCGa                                                                                         <  2:418945/93‑1 (MQ=255)

  aCCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATa                                                                                       >  2:3566381/1‑93 (MQ=255)

    cAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACg                                                                                     >  2:1500030/1‑93 (MQ=255)

     aGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGt                                                                                    <  2:2967558/93‑1 (MQ=255)

      ggTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtc                                                                                   <  2:882125/93‑1 (MQ=255)

        ttGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtctc                                                                                 >  1:3199688/1‑93 (MQ=255)

         tGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCg                                                                                >  1:3234893/1‑93 (MQ=255)

           tctcTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGtt                                                                              >  2:3788036/1‑93 (MQ=255)

            ctctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTa                                                                             <  1:4368441/93‑1 (MQ=255)

             tctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTAt                                                                            >  1:1444183/1‑93 (MQ=255)

             tctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTAt                                                                            >  1:202979/1‑93 (MQ=255)

             tctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGGCTCGTTAt                                                                            <  2:4580658/93‑1 (MQ=255)

              ctGTTATCACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATg                                                                           <  1:2139644/93‑1 (MQ=255)

              ctGTTAACACCCGATTCACACAGATAGAAAACAACGCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATg                                                                           <  2:1727333/93‑1 (MQ=255)

              ctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATg                                                                           <  1:629504/93‑1 (MQ=255)

               tGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATgg                                                                          <  1:4541216/93‑1 (MQ=255)

                  tAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGAcc                                                                       <  2:4139280/93‑1 (MQ=255)

                  tAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGAcc                                                                       <  1:3469366/93‑1 (MQ=255)

                  tAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGAAcc                                                                       <  2:732709/93‑1 (MQ=255)

                   aaCACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGAccc                                                                      <  1:3159983/93‑1 (MQ=255)

                    acacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCg                                                                     <  1:4114917/93‑1 (MQ=255)

                     cacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGt                                                                    >  1:3929210/1‑93 (MQ=255)

                     cacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGCTATGGACCCGt                                                                    >  2:623497/1‑93 (MQ=255)

                      acCCGATTCACACAGATAAAAAACAACACGGTATATTATGCTGCAACCCAAACGGGTGATTGCCAGGGAGATACGTCTCGTTATGGACCCGtt                                                                   >  2:3661811/1‑93 (MQ=255)

                        ccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTaa                                                                 <  1:1685785/93‑1 (MQ=255)

                        ccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTaa                                                                 >  1:4357265/1‑93 (MQ=255)

                         cGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAat                                                                >  2:811331/1‑93 (MQ=255)

                         cGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAat                                                                >  1:1254174/1‑93 (MQ=255)

                          gATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAata                                                               <  2:924410/93‑1 (MQ=255)

                          gATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAata                                                               <  1:2338238/93‑1 (MQ=255)

                           aTTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAatat                                                              <  1:4203798/93‑1 (MQ=255)

                            ttCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATAtt                                                             <  1:3508246/93‑1 (MQ=255)

                            ttCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATAtt                                                             <  2:673462/93‑1 (MQ=255)

                              cacacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTgg                                                           <  1:1394013/93‑1 (MQ=255)

                              cacacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTgg                                                           <  1:1534673/93‑1 (MQ=255)

                               acacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGt                                                          >  1:3083862/1‑93 (MQ=255)

                                cacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTg                                                         >  1:1405968/1‑93 (MQ=255)

                                 acaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGc                                                        <  1:1746847/93‑1 (MQ=255)

                                 acaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGc                                                        <  2:2356259/93‑1 (MQ=255)

                                 acaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGc                                                        >  2:3240549/1‑93 (MQ=255)

                                 acaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGc                                                        >  1:2759214/1‑93 (MQ=255)

                                  caGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCg                                                       >  2:3588234/1‑93 (MQ=255)

                                   aGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGa                                                      >  2:4012079/1‑93 (MQ=255)

                                    gATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGaa                                                     <  1:4231689/93‑1 (MQ=255)

                                    gATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGaa                                                     <  1:2173395/93‑1 (MQ=255)

                                    gATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGaa                                                     <  2:2399694/93‑1 (MQ=255)

                                      tAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAaca                                                   >  1:366640/1‑93 (MQ=255)

                                      tAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAaca                                                   >  1:4644229/1‑93 (MQ=255)

                                         aaaaCAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAAcacca                                                <  2:1565426/93‑1 (MQ=255)

                                           aacaacACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACaccacc                                              <  2:2259006/93‑1 (MQ=255)

                                           aacaacACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACaccacc                                              <  2:4323223/93‑1 (MQ=255)

                                            acaacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCt                                             <  1:434917/93‑1 (MQ=255)

                                            acaacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCt                                             >  2:2701021/1‑93 (MQ=255)

                                            acaacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCt                                             >  2:1785366/1‑93 (MQ=255)

                                            acaacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCt                                             >  2:4041689/1‑93 (MQ=255)

                                            acaacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCt                                             >  1:2151019/1‑93 (MQ=255)

                                               acacGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTa                                          >  2:2015489/1‑93 (MQ=255)

                                               acacGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTa                                          <  1:457278/93‑1 (MQ=255)

                                                cacGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAc                                         <  2:53668/93‑1 (MQ=255)

                                                 acGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAcc                                        >  1:1343706/1‑93 (MQ=255)

                                                  cGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAccc                                       <  1:65887/93‑1 (MQ=255)

                                                  cGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAccc                                       <  1:325242/93‑1 (MQ=255)

                                                   ggTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCt                                      >  1:2841912/1‑93 (MQ=255)

                                                   ggTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCt                                      >  1:3858713/1‑93 (MQ=255)

                                                    gTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTg                                     >  1:3346685/1‑93 (MQ=255)

                                                    gTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTg                                     >  1:2627602/1‑93 (MQ=255)

                                                     tataTTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGt                                    <  2:3842326/93‑1 (MQ=255)

                                                      atatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTa                                   >  2:3078043/1‑93 (MQ=255)

                                                      atatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTa                                   >  1:1111495/1‑93 (MQ=255)

                                                       tatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACATCACCTTTACCCTGTAt                                  >  2:3988547/1‑93 (MQ=255)

                                                       tatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTAt                                  >  2:3217825/1‑93 (MQ=255)

                                                        atTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATg                                 <  2:2795344/93‑1 (MQ=255)

                                                        atTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATg                                 >  2:2293921/1‑93 (MQ=255)

                                                        atTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATg                                 >  1:1194060/1‑93 (MQ=255)

                                                        atTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATg                                 >  2:3907121/1‑93 (MQ=255)

                                                          tCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTc                               >  1:4599305/1‑93 (MQ=255)

                                                           ctgctgCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCa                              >  2:4627401/1‑93 (MQ=255)

                                                             gctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCAcc                            >  2:3287065/1‑93 (MQ=255)

                                                             gctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCAcc                            >  1:4559402/1‑93 (MQ=255)

                                                                gcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAg                         >  1:2997422/1‑93 (MQ=255)

                                                                    cccAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAt                     >  2:3685102/1‑93 (MQ=255)

                                                                    cccAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAt                     <  2:1209242/93‑1 (MQ=255)

                                                                      cAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAttt                   <  2:4278733/93‑1 (MQ=255)

                                                                       aaaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtc                  <  2:1080783/93‑1 (MQ=255)

                                                                       aaaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtc                  >  2:3442353/1‑93 (MQ=255)

                                                                       aaaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATGGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtc                  <  1:1494383/93‑1 (MQ=255)

                                                                         aCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtctc                <  2:2676344/93‑1 (MQ=255)

                                                                         aCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtctc                <  1:4323652/93‑1 (MQ=255)

                                                                          cGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCg               <  1:1669789/93‑1 (MQ=255)

                                                                           gggTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCgg              <  2:985389/93‑1 (MQ=255)

                                                                            ggTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGc             <  2:2208944/93‑1 (MQ=255)

                                                                             gtgtTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCt            >  1:790064/1‑93 (MQ=255)

                                                                             gtgtTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCt            >  1:1060987/1‑93 (MQ=255)

                                                                              tgtTTGCCAGGGCGATACGTCTCGTAATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTc           <  2:1979915/93‑1 (MQ=255)

                                                                                      aGGGGGATACGTCTTGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   <  1:1709233/93‑1 (MQ=255)

                                                                                      aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   >  1:4000626/1‑93 (MQ=255)

                                                                                      aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   >  1:3578366/1‑93 (MQ=255)

                                                                                       gggCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAt  <  1:3224288/93‑1 (MQ=255)

                                                                                       gggCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAt  <  2:1308406/93‑1 (MQ=255)

CGACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGTGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAT  >  NC_000913/3289341‑3289520

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