BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,900,198	G→A	G94D (GGC→GAC)	ygcE →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,900,198	0	G	A	99.2%	1386.8 / ‑3.5	396	G94D (GGC→GAC)	ygcE	putative kinase
Reads supporting (aligned to +/- strand): ref base G (3/0); new base A (204/188); total (207/189)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 426 total aligned reads displayed.

TGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCT  >  NC_000913/2900108‑2900287

tgatgaGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACa                                                                                         <  2:6889470/93‑1 (MQ=255)

tgatgaGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACa                                                                                         >  1:962961/1‑93 (MQ=255)

 gatgaGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGacac                                                                                        <  2:8027300/93‑1 (MQ=255)

   tgaGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGc                                                                                      <  2:2346326/93‑1 (MQ=255)

     aGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCct                                                                                    >  2:2671214/1‑93 (MQ=255)

     aGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCct                                                                                    <  1:590255/93‑1 (MQ=255)

       tCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCctgc                                                                                  <  1:1996414/93‑1 (MQ=255)

       tCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCctgc                                                                                  >  2:8562282/1‑93 (MQ=255)

         aGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACCActgctc                                                                                >  1:1501195/1‑92 (MQ=255)

          gTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCtgctgc                                                                               >  1:1163638/1‑93 (MQ=255)

              gCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCg                                                                           >  1:6495438/1‑93 (MQ=255)

              gCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTAAAAGCCGATGACACGCCTGCTGCGCCg                                                                           >  2:4661020/1‑93 (MQ=255)

                  ggAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTgt                                                                       >  1:2669577/1‑92 (MQ=255)

                  ggAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGa                                                                       >  1:6852663/1‑93 (MQ=255)

                  ggAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGa                                                                       >  1:3752432/1‑93 (MQ=255)

                    aaTAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGAtt                                                                     <  2:828485/93‑1 (MQ=255)

                     aTAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTa                                                                    >  2:1048841/1‑93 (MQ=255)

                       aaaGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGc                                                                  >  1:7830498/1‑93 (MQ=255)

                        aagaagATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCt                                                                 >  1:6639573/1‑93 (MQ=255)

                         agaagaTATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTg                                                                <  2:1945733/93‑1 (MQ=255)

                         agaagaTATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTg                                                                >  2:1833570/1‑93 (MQ=255)

                         agaagaTATTGTCGGTATTGGTCTGGGATCCATACGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTACTGCGCCGTTGATTAGCTg                                                                <  1:8967628/93‑1 (MQ=255)

                          gaagaTATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTgg                                                               >  2:736411/1‑93 (MQ=255)

                          gaagaTATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACACCAGCTGGGCCGTTGATTAGCTgg                                                               >  1:1663007/1‑93 (MQ=255)

                            agaTATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCa                                                             >  2:3302158/1‑93 (MQ=255)

                             gaTATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAg                                                            <  1:1225165/93‑1 (MQ=255)

                             gaTATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAg                                                            >  2:7109788/1‑93 (MQ=255)

                               tatTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGa                                                          <  2:3182885/93‑1 (MQ=255)

                               tatTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGa                                                          >  2:4882064/1‑93 (MQ=255)

                               tatTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGa                                                          <  2:6621235/93‑1 (MQ=255)

                                 ttGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATg                                                        <  1:7579084/93‑1 (MQ=255)

                                 ttGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATg                                                        >  2:2992358/1‑93 (MQ=255)

                                 ttGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATg                                                        <  1:4888622/93‑1 (MQ=255)

                                   gTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCa                                                      <  1:6183475/93‑1 (MQ=255)

                                   gTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCa                                                      <  1:4040048/93‑1 (MQ=255)

                                   gTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCa                                                      >  2:8759001/1‑93 (MQ=255)

                                   gTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCa                                                      <  2:2052686/93‑1 (MQ=255)

                                      ggTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCAcgc                                                   >  1:4236337/1‑93 (MQ=255)

                                       gTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCAcgcg                                                  <  2:1718392/93‑1 (MQ=255)

                                          ttGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTa                                               >  1:7149644/1‑93 (MQ=255)

                                          ttGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTa                                               <  2:786395/93‑1 (MQ=255)

                                            ggTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTaca                                             >  2:5254772/1‑93 (MQ=255)

                                                tGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGcc                                         >  1:2894460/1‑93 (MQ=255)

                                                  ggATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCt                                       >  1:5713980/1‑93 (MQ=255)

                                                   gATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCtt                                      >  1:9481283/1‑93 (MQ=255)

                                                    atccatccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTa                                     <  1:8394469/93‑1 (MQ=255)

                                                    atccatccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTa                                     >  2:2501961/1‑93 (MQ=255)

                                                    atccatccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTa                                     <  2:7700716/93‑1 (MQ=255)

                                                     tccatccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTAc                                    <  1:1398807/93‑1 (MQ=255)

                                                      ccatccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACg                                   >  1:7696723/1‑93 (MQ=255)

                                                        atccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTCCACGCCCTTACGaa                                 <  2:821997/93‑1 (MQ=255)

                                                        atccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGaa                                 <  1:7494877/93‑1 (MQ=255)

                                                         tccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCCCGCGTTACACGCCCTTACGAac                                >  2:1754456/1‑93 (MQ=255)

                                                         tccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAac                                >  1:9135850/1‑93 (MQ=255)

                                                         tccGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAac                                >  2:561233/1‑93 (MQ=255)

                                                           cGTTGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGCacac                              >  2:3572256/1‑93 (MQ=255)

                                                             ttGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAacacac                            <  2:3814899/93‑1 (MQ=255)

                                                             ttGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAacacac                            <  1:8461424/93‑1 (MQ=255)

                                                             ttGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAacacac                            >  1:1713546/1‑93 (MQ=255)

                                                             ttGCTGCCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAacacac                            >  1:4578480/1‑93 (MQ=255)

                                                              tgctgcCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACt                           <  1:9947408/93‑1 (MQ=255)

                                                              tgctgcCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACt                           >  2:1378633/1‑93 (MQ=255)

                                                              tgctgcCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACt                           <  2:4040670/93‑1 (MQ=255)

                                                              tgctgcCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACt                           <  1:4452016/93‑1 (MQ=255)

                                                                ctgcCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTaa                         <  1:3367693/93‑1 (MQ=255)

                                                                 tgcCGTGCGTTATTGAAAGCCGATGACACGCCTGCTTCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAAc                        <  2:7446614/93‑1 (MQ=255)

                                                                  gcCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAAcc                       <  1:10014694/93‑1 (MQ=255)

                                                                  gcCGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGCTACACGCCCTTACGAACACACTAAcc                       >  1:5374833/1‑93 (MQ=255)

                                                                   ccGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAAccc                      <  1:8585474/93‑1 (MQ=255)

                                                                   ccGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAAccc                      <  1:3107547/93‑1 (MQ=255)

                                                                   ccGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAAccc                      <  2:47235/93‑1 (MQ=255)

                                                                    cGTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCt                     <  1:9762902/93‑1 (MQ=255)

                                                                     gTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTg                    <  2:5796838/93‑1 (MQ=255)

                                                                     gTGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTg                    >  1:390480/1‑93 (MQ=255)

                                                                      tccGTGATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTAACACGCCCTTACGAAGACACGATCCCTGa                   <  2:6219601/91‑1 (MQ=255)

                                                                      tGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGa                   >  1:5279593/1‑93 (MQ=255)

                                                                      tGCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGa                   <  1:7102024/93‑1 (MQ=255)

                                                                      tGCGTTATTGAAAGCCGATGACACGCCTGCTGCGACGTTGATTAGCTGGCGGGATGCACGCGTTACGGGCCCTTACGAACACACTACCCCttt                   >  1:3537857/1‑91 (MQ=255)

                                                                       gCGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGAc                  <  1:5783005/93‑1 (MQ=255)

                                                                        cGTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACg                 <  2:4238617/93‑1 (MQ=255)

                                                                        agTTATAGAACGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACTCACTAACCCAGACg                 <  2:9210965/92‑1 (MQ=255)

                                                                         gTTATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGt                <  1:2471002/93‑1 (MQ=255)

                                                                          ttATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTg               <  2:1221744/93‑1 (MQ=255)

                                                                          ttATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTg               <  1:8796731/93‑1 (MQ=255)

                                                                          ttATTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTg               >  2:6312245/1‑93 (MQ=255)

                                                                            aTTGAAAGCCGATGACACGCCTGCTGTGCAGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGc             <  2:9707531/93‑1 (MQ=255)

                                                                            aTTGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATCCACGAGTTACACGCCCTTACGCACACACTAACGCTGATGTGGc             >  1:5009522/1‑93 (MQ=255)

                                                                             ttGAAAGCCGATGACACGCCTGCTGCGCCGTTTATTAGCTGGCAGGGTGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCa            >  2:2839930/1‑93 (MQ=255)

                                                                             ttGAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCa            >  1:6343400/1‑93 (MQ=255)

                                                                               gAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCata          >  2:6007252/1‑93 (MQ=255)

                                                                               gAAAGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCata          <  2:7593397/93‑1 (MQ=255)

                                                                                aaaGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCatat         >  1:7004636/1‑93 (MQ=255)

                                                                                aaaGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCatat         <  1:5526039/93‑1 (MQ=255)

                                                                                 aaGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATg        <  1:8048031/93‑1 (MQ=255)

                                                                                  aGCCGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGt       <  1:2839013/93‑1 (MQ=255)

                                                                                    ccGATGACGCGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCa     <  1:2493194/93‑1 (MQ=255)

                                                                                    ccGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCa     <  1:9383295/93‑1 (MQ=255)

                                                                                     cGATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCAc    >  2:9087894/1‑93 (MQ=255)

                                                                                      gATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCAcc   >  2:4461823/1‑93 (MQ=255)

                                                                                      gATGACACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCAcc   >  1:5096298/1‑93 (MQ=255)

TGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACACACTAACCCTGACGTGGCATATGTCACCT  >  NC_000913/2900108‑2900287

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