BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	708,421	G→A	G289S (GGC→AGC)	ybfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	708,421	0	G	A	84.7%	672.1 / 79.3	248	G289S (GGC→AGC)	ybfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (20/18); new base A (113/97); total (133/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 250 total aligned reads displayed.

TGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGAC  >  NC_000913/708330‑708510

 gATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGc                                                                                         <  2:3614776/93‑1 (MQ=255)

  atatCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCa                                                                                        <  1:4571049/93‑1 (MQ=255)

   tatCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCAc                                                                                       >  1:4161704/1‑93 (MQ=255)

   tatCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCAc                                                                                       <  2:5435976/93‑1 (MQ=255)

   tatCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCAc                                                                                       >  2:1923187/1‑93 (MQ=255)

   tatCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCAc                                                                                       >  2:4001905/1‑93 (MQ=255)

     tCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCg                                                                                     >  2:6962273/1‑93 (MQ=255)

       gCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGcc                                                                                   <  1:6941929/93‑1 (MQ=255)

         cGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTg                                                                                 >  1:2544188/1‑93 (MQ=255)

         cGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTg                                                                                 <  1:3907177/93‑1 (MQ=255)

         cGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTg                                                                                 >  1:6757547/1‑93 (MQ=255)

         cGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTg                                                                                 >  2:715706/1‑93 (MQ=255)

         cGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCAAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTg                                                                                 <  2:5875917/93‑1 (MQ=255)

          ggTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTgg                                                                                >  2:5972552/1‑93 (MQ=255)

           gTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGc                                                                               >  2:390220/1‑93 (MQ=255)

           gTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGc                                                                               >  1:5917384/1‑93 (MQ=255)

             aGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTg                                                                             <  2:3398793/93‑1 (MQ=255)

              gCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGc                                                                            >  2:7212958/1‑93 (MQ=255)

                 cGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAg                                                                         >  1:4275642/1‑93 (MQ=255)

                 cGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAg                                                                         >  2:4936928/1‑93 (MQ=255)

                 cGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAg                                                                         >  1:2036532/1‑93 (MQ=255)

                  gTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGc                                                                        >  2:6352757/1‑93 (MQ=255)

                   ttAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCg                                                                       >  2:4139033/1‑93 (MQ=255)

                   ttAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCg                                                                       >  1:1217022/1‑93 (MQ=255)

                   ttAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCg                                                                       <  2:162165/93‑1 (MQ=255)

                    tAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGt                                                                      >  1:4377377/1‑93 (MQ=255)

                      accaccAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTg                                                                    <  2:7483621/93‑1 (MQ=255)

                      accaccAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTg                                                                    >  2:5218472/1‑93 (MQ=255)

                       ccaccaGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGa                                                                   >  2:4298032/1‑93 (MQ=255)

                         accaGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGAcc                                                                 >  1:4570937/1‑93 (MQ=255)

                         accaGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGAcc                                                                 <  2:4381235/93‑1 (MQ=255)

                          ccaGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCt                                                                <  2:861392/93‑1 (MQ=255)

                             gCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGTCTGCAAGCGTTGACCTTTg                                                             >  2:958828/1‑93 (MQ=255)

                              cTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTgg                                                            <  2:3459568/93‑1 (MQ=255)

                                aCCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGtt                                                          <  2:263501/93‑1 (MQ=255)

                                 ccAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTa                                                         <  2:7020797/93‑1 (MQ=255)

                                 ccAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTa                                                         <  2:2211401/93‑1 (MQ=255)

                                 ccAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTa                                                         >  1:7191944/1‑93 (MQ=255)

                                  cAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTAc                                                        <  1:1627645/93‑1 (MQ=255)

                                  cAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTAc                                                        <  1:6276927/93‑1 (MQ=255)

                                  cAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTAc                                                        >  2:1286725/1‑93 (MQ=255)

                                   aGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTAcc                                                       >  2:5649903/1‑93 (MQ=255)

                                    gTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCg                                                      >  2:3007565/1‑93 (MQ=255)

                                    gTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCg                                                      <  2:4237486/93‑1 (MQ=255)

                                     ttCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCgg                                                     >  2:5342200/1‑93 (MQ=255)

                                     ttCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCgg                                                     <  2:781880/93‑1 (MQ=255)

                                       cTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGggc                                                   <  2:7397731/93‑1 (MQ=255)

                                       cTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGggc                                                   >  1:5151108/1‑93 (MQ=255)

                                        tACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGggcg                                                  >  2:6804817/1‑93 (MQ=255)

                                           ggTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCt                                               >  1:869476/1‑93 (MQ=255)

                                              aCCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGAc                                            <  1:7084263/93‑1 (MQ=255)

                                               cccGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACg                                           <  2:2898351/93‑1 (MQ=255)

                                                 cGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTa                                         <  1:57918/93‑1 (MQ=255)

                                                 cGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTa                                         <  1:2920853/93‑1 (MQ=255)

                                                  gagaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAg                                        >  1:1271314/1‑93 (MQ=255)

                                                  gagaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAg                                        >  1:1717499/1‑93 (MQ=255)

                                                  gagaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAg                                        >  2:3520266/1‑93 (MQ=255)

                                                   agaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGt                                       >  1:4973748/1‑93 (MQ=255)

                                                   agaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGt                                       >  1:2135398/1‑93 (MQ=255)

                                                    gaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTg                                      <  1:5193581/93‑1 (MQ=255)

                                                      tAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGa                                    <  1:4825153/93‑1 (MQ=255)

                                                      tAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGa                                    >  1:3022938/1‑93 (MQ=255)

                                                        aaGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGAtt                                  <  2:2446256/93‑1 (MQ=255)

                                                         aGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGAttt                                 >  2:4830446/1‑93 (MQ=255)

                                                           ttGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTgc                               >  2:3818425/1‑93 (MQ=255)

                                                            tGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTgcg                              >  2:407427/1‑93 (MQ=255)

                                                             gACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTgcgc                             <  1:237239/93‑1 (MQ=255)

                                                              aCGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGTCCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGcc                            <  1:7510183/93‑1 (MQ=255)

                                                                gATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTc                          <  1:638454/93‑1 (MQ=255)

                                                                 aTCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCg                         >  1:376879/1‑93 (MQ=255)

                                                                 aTCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCg                         >  2:1104187/1‑93 (MQ=255)

                                                                  tCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGa                        <  1:6566438/93‑1 (MQ=255)

                                                                   cGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGaa                       <  2:3703507/93‑1 (MQ=255)

                                                                   cGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGaa                       <  1:346405/93‑1 (MQ=255)

                                                                    gCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAg                      <  1:647899/93‑1 (MQ=255)

                                                                     cAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAgg                     <  2:5139449/93‑1 (MQ=255)

                                                                     cAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAgg                     <  2:123112/93‑1 (MQ=255)

                                                                      aGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGc                    <  1:5971750/93‑1 (MQ=255)

                                                                      aGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGc                    >  2:1974452/1‑93 (MQ=255)

                                                                       gCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCa                   >  1:2337427/1‑93 (MQ=255)

                                                                         gTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCAcc                 >  2:2508722/1‑93 (MQ=255)

                                                                          tCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCt                >  1:2963391/1‑93 (MQ=255)

                                                                           cAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTg               <  1:1918980/93‑1 (MQ=255)

                                                                           cAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTg               <  2:4471178/93‑1 (MQ=255)

                                                                            aaCGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTgg              <  2:2797935/93‑1 (MQ=255)

                                                                             aCGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTggg             <  2:5557734/93‑1 (MQ=255)

                                                                               gACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGtt           >  1:5427110/1‑93 (MQ=255)

                                                                                aCCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTa          >  2:943427/1‑93 (MQ=255)

                                                                                aCCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTa          >  2:2254480/1‑93 (MQ=255)

                                                                                 ccTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTaa         <  1:3841385/93‑1 (MQ=255)

                                                                                 ccTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTaa         >  1:1187323/1‑93 (MQ=255)

                                                                                 ccTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTaa         <  1:2607145/93‑1 (MQ=255)

                                                                                  cTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAg        >  2:6714368/1‑93 (MQ=255)

                                                                                  cTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAg        >  1:1486431/1‑93 (MQ=255)

                                                                                   tttATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAgg       >  1:7375168/1‑93 (MQ=255)

                                                                                     tATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCt     >  1:1968712/1‑93 (MQ=255)

                                                                                       tGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGa   >  2:1550745/1‑93 (MQ=255)

                                                                                        gACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGAc  <  2:4871000/93‑1 (MQ=255)

                                                                                        gACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGAc  >  2:6276861/1‑93 (MQ=255)

                                                                                        gACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGAc  <  2:4747431/93‑1 (MQ=255)

TGATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGAC  >  NC_000913/708330‑708510

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